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COMPUTATION  METHODS 
BASIC  FUNCTIONS 


^INTRODUCTION 


The  computation  methods  described  in  this  section  are  for  use  in  fatigue 
analyses.  The  functions  to  be  computed  are  the  Gamma  Function,  the  two 
Incomplete  Gamma  Functions,  an  Error  Function  and  its  Inverse,  a 
Probability  of  Failure  Function,  and  a  transcendental  accelerated  test 
level  function. 


These  methods  are  intended  to  be  user  oriented.  The  user  is  given  a 
choice  of  methods.  Tables  and  interpolation  methods  are  included  for 
those  who  prefer  to  use  tables.  Techniques  are  also  included  for  use 
with  calculators,  with  the  modern  TI-59,  HP-67,  HP-34C  and  HP-4 1C 
Programmable  Calculators  and  with  large  computers  (i.e.  Basic  Language). 
Both  numerical  integration  and  closed  form  equation  methods  are  given. 
Examples  using  each  method  are  worked  out.  Computer  program  listings 
are  also  shown. 


It  should  be  noted  that  the  same  programs  can  be  used  for  both  the 
HP-67  and  HP-4 1C  programmable  calculators.  Only  the  HP-67  is  referred 
to  in  the  listings  for  simplicity. 
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A.  GAMMA  FUNCTION  COMPUTATION 


Definition: 

The  Gamma  Function  is  defined  as  follows: 


Cw3 


T(a)  is  undefined  for  ct=  0  and  for  negative  integer  values  of  a. 
For  fatigue  analyses  a  >  1.  Figure  1  shows  T( a)  versus  a  for 
0  <  a  <  4,  The  curve  increases  monotonically  for  a  >  4. 


Figure  1  Plot  of  Gamma  Function  Versus  a 
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TABULAR  METHOD 


Table  I  (Table  6.1  pages  267-270  of  Abramowitz  j_l  _  )  gives  values  to  ten 
places  of  T  (a)  versus  a  for  1  <_  a  <_  2.  Table  II  (Table  6.3  pages  272, 
273  of  Abramowitz  £l  j  )  gives  values  to  eleven  places  of  T  (a)  for 
integer  and  half-integer  values  of  a  for  1  <  a  _<  101. 

EXAMPLE 

FIND:  r  (1.225) 

SOLUTION:  In  Table  I  (Table  6.1  page  267  [_lj  ) 

T  (1.225)  =  0.9119156071 

EXAMPLE: 

FIND:  I  (8.5) 

SOLUTION:  In  Table  II  (Table  6.3  page  272  [l]  ) 

r  (8.5)  =  1.4034407x10^ 

The  following  expression  can  be  used  recursively: 
r  (a)  =  (a-1)  r  (a-T) 

EXAMPLE : 

FIND:  T  (5.64) 

SOLUTION:  T  (5.64)  =  4.64  x  3.64  x  2.64  x  1.64  x  F(1.64) 

In  Table  6.1  page  269  1  i 

T(1.64)  =  0.8986420302 
1(5.64)  =  65.71338911 
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Table  I  GAMMA  FUNCTION  (1  £  a  <  2) 


GAMMA  FUNCTION  AND  RELATED  FUNCTIONS 
GAMMA,  DIGAMMA  AND  TRIGAMMA  FUNCTIONS  Table  6.1 


M 

r« 

In  r(e) 

*(z) 

*■(') 

L000 

1.00000  00000 

0.00000 

00000 

-0.57721 

56649 

1.64493 

40668 

0.000 

1.005 

0.99713 

85354 

-0.00286 

55666 

-0.56902  09113 

1.63299 

41567 

0.005 

1.010 

0.99432 

58512 

-0.00569 

03079 

-0.56088 

54579 

1.62121 

35283 

0.010 

1.015 

0.99156 

12888 

-0.00847 

45187 

-0.55280 

85156 

1.60958  91824 

0.015 

1.070 

0.98884 

42033 

-0.01121 

84893 

-0.54478  93105 

1.59811 

81919 

0.020 

1.025 

0.98617 

39633 

-0.01392 

25067 

-0.53682 

70828 

1.58679 

76993 

0.025 

1.030 

0.98354 

99506 

-0.01658 

68539 

-0.52892 

10873 

1.57562 

49154 

0.030 

1.035 

0.98097 

15606 

-0.01921 

18101 

-0.52107 

05921 

1.56459 

71163 

0.035 

1.040 

0.97843 

82009 

-0.02179 

76511 

-0.51327 

48789 

1.55371 

16426 

0.040 

1.045 

0.97594 

92919 

-0.02434 

46490 

-0.50553 

32428 

1.54296 

58968 

0.045 

1.050 

0.97350 

42656 

-0.02685 

30725 

-0.49784 

49913 

1.53235 

73421 

0.050 

1.055 

0.97110  25663 

-0.02932 

31868 

-0.49020 

94448 

1.52188 

35001 

0.055 

1.060 

0.96874 

36495 

-0.03175 

52537 

-0.48262 

59358 

1.51154 

19500 

0.060 

1.065 

0.96642 

69823 

-0.03414 

95318 

-0.47509 

38088 

1.50133 

03259 

0.065 

1.070 

0.96415 

20425 

-0.03650 

62763 

-0.46761 

24199 

1.49124 

63164 

0.070 

1.075 

0.96191 

83189 

-0.03882 

57395 

-0.46018 

11367 

1.48128 

76622 

0.075 

1.080 

0.95972 

53107 

-0.04110 

81702 

-0.45279 

93380 

1.47145 

21556 

0.080 

1.085 

0.95757  25273 

-0.04335 

38143 

-0.44546  64135 

1.46173 

76377 

0.085 

1.090 

0.95545 

94882 

-0.04556  29148 

-0.43818 

17635 

1.45214 

19988 

0.090 

1.095 

0.95338 

57227 

-0.04773 

57114 

-0.43094  47988 

1.44266 

31755 

0.095 

1.100 

0.95135 

07699 

-0.04987 

24413 

-0.42375 

49404 

1.43329 

91508 

0.100 

1.105 

0.94935 

41778 

-0.05197 

33384 

-0.41661 

16193 

1.42404 

79514 

0.105 

1.110 

0.94739 

55040 

-0.05403 

86341 

-0.40951 

42761 

1.41490 

76482 

0.110 

1.115 

0.94547 

43149 

-0.05606 

85568 

-0.40246  23611 

1.40587 

63535 

0.115 

1.120 

0.94359 

01856 

-0.05806 

33325 

-0.39545 

53339 

1.39695 

22213 

0.120 

1.125 

0.94174 

26997 

-0.06002 

31841 

-0.38849 

26633 

1.38813 

34449 

0.125 

1.130 

0.93993 

14497 

-0.06194 

83322 

-0.38157 

38268 

1.37941 

82573 

0.130 

1.135 

0.93815 

60356 

-0.06383 

89946 

-0.37469 

83110 

1.37080 

49288 

0.135 

1.140 

0.93641 

60657 

-0.06569  53867 

-0.36786  56106 

1.36229 

17670 

0.140 

1.145 

0.93471 

11562 

-0.06751 

77212 

-0.36107 

52291 

1.35387 

71’52 

0.145 

1.150 

0.93304 

09311 

-0.06930 

620B7 

-0.35432 

66780 

1.34555 

93520 

0.150 

1.155 

0.93140 

50217 

-0.07106 

10569 

-0.34761 

94768 

1.33733 

68900 

0.155 

1.160 

0.92980 

30666 

-0.07278 

24716 

-0.34095 

31528 

1.32920 

81752 

0.160 

1.165 

0.92823 

47120 

-0.07447 

06558 

-0.33432 

72413 

1.32117 

16859 

0.165 

1.170 

0.92669  96106 

-0.07612 

58106 

-0.32774 

12847 

1.31322 

59322 

0.170 

1,175 

0.92519 

74225 

-0.07774 

81345 

-0.32119 

48332 

1.30536 

94548 

0.175 

1.180 

0.92372 

78143 

-0.07933 

78240 

-0.31468 

74438 

1.29760 

08248 

0.180 

1.185 

0.92229 

04591 

-0.08089 

50733 

-0.30821 

86809 

1.28991 

86421 

0.185 

1.190 

0.92088 

50371 

-0.08242 

00745 

-0.30178 

81156 

1.28232 

15358 

0.190 

1.195 

0.91951 

12341 

-0.08391 

30174 

-0.29539 

53259 

1.27480 

81622 

0.195 

1.200 

0.91816 

87424 

-0.08537 

40900 

-0.28903  98966 

1.26737 

72054 

0.200 

1.205 

0.91685 

72606 

-0.08680 

34780 

-0.28272 

14187 

1.26002 

73755 

0.205 

1.210 

0.91557  64930 

-0.08820 

13651 

-0.27643  94897 

1.25275 

74090 

0.210 

1.215 

0.91432 

61500 

-0.08956 

79331 

-0.27019 

37135 

1.24556 

60671 

0.215 

1.22Q 

0.91310  59475 

-0.09090 

33619 

-0.26398 

37000 

1.23845 

21360 

0.220 

1.225 

1.230 

1.235 

0.91191 

56071 

-0.09220 

78291 

-0.25780 

90652 

1.23141 

44258 

0.225 

0.91075 

48564 

-0.09348 

15108 

-0.25166  94307 

1.22445 

17702 

0.230 

0.90962 

34274 

-0.09472  45811 

-0.24556  44243 

1.21756 

30254 

0.235 

1.240 

1.245 

0.90852 

10583 

-0.09593 

72122 

-0.23949 

36791 

1.21074 

70707 

0.240 

0.90744 

74922 

-0.09711 

95744 

-0.23345 

68341 

1.20400 

28063 

0.245 

1.250 

0.90640 

24771 

-0.09827 

18364 

-0.22745 

35334 

1.19732 

91545 

0.250 

y 

I 

In  y! 

ih*' 

u 

r('6)61 

f(-6)51 

n-6)7i 

f(-5)21 

l  & 

J 

l  5 

•  J 

l  5 

J 

l  5 

J 

Fw  '>2  see  Example*  1-1.  logl°  ""0  '13429  44819 


^"piled  from  H.  T.  Davis,  Tables  of  the  higher  mathematical  functions,  2  vols.  (Principia 
*^«ss,  Bloomington,  Ind.,  1933, 1935)  (with  permission).  Known  error  lias  been  corrected. 
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Table  6.1 


Table  I  (Cont'd) 

GAMMA  FUNCTION  AND  It)  LATKD  FUNCTIONS 

<;amma,  dicamma  \m»  tkh.\mma  u  nctions 


I 

1.250 

1.255 

1.260 

1.265 

1.270 

r(-r) 

0.90640  24771 
0.90538  57663 
0.90439  71178 
0.90343  62946 
0.90250  30645 

In  r(j-) 

-0.09827  18364 
-0.09939  41651 
-0.10048  67254 
-0.10154  96809 
-0.10258  31932 

i(j) 

-0.22745  35334 
-0.22148  34266 
-0.21554  61686 
-0.20964  14193 
-0.20376  88437 

1.19732  91545 
1.19072  50579 
1.18418  94799 
1.17772  14030 
1.17131  98301 

0.250 

0.255 

0.260 

0.265 

0.270 

1.275 

1.280 

1.285 

1.290 

1.295 

0.90159  71994 
0.90071  84765 
0.89986  66769 
0.89904  15863 
0.89824  29947 

-0.10358  74224 
-0.10456  25269 
-0.10550  86634 
-0.10642  59872 
-0.10731  46519 

-0.19792  81118 
-0.19211  88983 
-0.18634  08828 
-0.18059  37494 
-0.17487  71870 

1.16498  37821 
1.15871  22990 
1.15250  44385 
1.14635  92764 
1.14027  59053 

0.275 

0.280 

0.285 

0.290 

0.295 

1.300 

1.305 

1.310 

1.315 

1.320 

0.89747-  06963 
0.89672  44895 
0.89600  41767 
0.89530  95644 
0.89464  04630 

-0.10817  48095 
-0.10900  66107 
-0.10981  02045 
-0.11058  57384 
-0.11133  33587 

-0.16919  08889 
-0.16353  45526 
-0.15790  78803 
-0.15231  05782 
-0.14674  23568 

1.13425  34350 
1.12829  09915 
1.12238  77175 
1.11654  27706 
1.11075  53246 

0.300 

0.305 

0.310 

0.315 

0.320 

1.325 

1.330 

1.335 

1.340 

1.345 

0.89399  66866 
0.89337  80535 
0.89278  43850 
0.89221  55072 
0.89167  12485 

-0.11205  32100 
-0.11274  54356 
-0.11341  01772 
-0.11404  75756 
-0.11465  77697 

-0.14120  29305 
-0.13569  20180 
-0.13020  93416 
-0.12475  46279 
-0.11932  76069 

1.10502  45678 
1.09934  97037 
1.09372  99497 
1.08816  45379 
1.08265  27136 

0.325 

0.330 

0.335 

0.340 

0.345 

1.350 

1.355 

1.360 

1.365 

1.370 

0.89115  14420 
0.89065  59235 
0.89018  45324 
0.88973  71116 
0.88931  35074 

-0.11524  08974 
-0.11579  70951 
-0.11632  64980 
-0.11682  92401 
-0.11730  54539 

-0.11392  80127 
-0.10855  55827 
-0.10321  00582 
-0.09789  11840 
-0.09259  87082 

1.07719  37361 
1.07178  68773 
1.06643  14226 
1.06112  66696 
1.05587  19286 

0.350 

0.355 

0.360 

0.365 

0.370 

1.375 

1.380 

1.385 

1.390 

1.395 

0.88891  35692 
0.88853  71494 
0.88818  41041 
0.88785  42918 
0.88754  75748 

-0.11775  52707 
-0.11817  88209 
-0.11857  62331 
-0.11894  76353 
-0.11929  31538 

-0.08733  23825 
-0.08209  19619 
-0.07687  72046 
-0.07168  78723 
-0.06652  37297 

1.05066  65216 
1.04550  97829 
1.04040  10578 
1.03533  97036 
1.03032  50881 

0.375 

0.380 

0.385 

0.390 

0.395 

1.400 

1.405 

1.410 

1.415 

1.420 

0.88726  38175 
0.88700  28884 
0.88676  46576 
0.88654  89993 
0.88635  57896 

-0.11961  29142 
-0.11990  70405 
-0.12017  56559 
-0.12041  88823 
-0.12063  68406 

-0.06138  45446 
-0.05627  00879 
-0.05118  01337 
-0.04611  44589 
-0.04107  28433 

1.02535  65905 
1.02043  36002 
1.01555  55173 
1.01072  17518 
1.00593  17241 

0.400 

0.405 

0.410 

0.415 

0.420 

1.425 

1.430 

1.435 

1.440 

1.445 

0.88618  49081 
0.88603  62361 
0.88590  96587 
0.88580  50635 
0.88572  23397 

-0.12082  96505 
-0.12099  74307 
-0.12114  02987 
-0.12125  83713 
-0.12135  17638 

-0.03605  50697 
-0.03106  09237 
-0.02609  01935 
-0.02114  26703 
-0.01621  81479 

1.00118  48640 
0.99648  06113 
0.99181  84147 
0.98719  77326 
0.98261  80318 

0.425 

0.430 

0.435 

0.440 

0.445 

1.450 

1.455 

1.460 

1.465 

1.470 

0.88566  13803 
0.88562  20800 
0.88560  43364 
0.88560  80495 
0.88563  31217 

-0.12142  05907 
-0.12146  49657 
-0.12148  50010 
-0.12148  08083 
-0.12145  24980 

-0.01131  64226 
-0.00643  72934 
-0.00158  05620 
♦  0.00325  39677 
0.00806  64890 

0.97807  87886 
0.97357  94874 
0.96911  96215 
0.96469  86921 
0.96031  62091 

0.450 

0.455 

0.460 

0.465 

0.470 

1.475 

1.480 

1.485 

1.490 

1.495 

0.88567  94575 
0.88574  69646 
0.88583  55520 
0.88594  51316 
0.88607  56174 

-0.12140  01797 
-0.12132  39621 
-0.12122  39528 
-0.12110  02585 
-0.12095  29852 

0.01285  71930 
0.01762  62684 
0.02237  39013 
0.02710  02758 
0.03180  55736 

0.95597  16896 
0.95166  46592 
0.94739  46509 
0.94316  12052 
0.93896  38700 

0.475 

0.480 

0.485 

0.490 

0.495 

1.500 

0.88622  69255 

A 

[‘I’4] 

-0.12078  22376 

In  y\ 

m 

0.03648  99740 

[T] 

0.93480  22005 

m 

0.500 

tl 

log,o  <’=0.43429  44819 

•S «  p «*«  n. 
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GAMMA,  DIGAMMA  AND  TWIG  VMM  A  U  NCTIONS  (  Table  6.1 


I 

1.500 

1.505 

1.510 

1.515 

1.520 

r(j) 

0.88622  69255 
0.88639  89744 
0.88659  16850 
0.88680  49797 
0.88703  87833 

In  r(i) 

-0.12078  22376 
-0.12058  81200 
-0.12037  07353 
-0.12013  01860 
-0.11986  65735 

<H  r) 

0.03648  99740 
0.04115  36543 
0.04579  67896 
0.05041  95527 
0.05502  21146 

*'U) 

0.93480  22005 
0.93067  57588 
0.92658  41142 
0.92252  68425 
0.91850  35265 

0.500 

0.505 

0.510 

0.515 

0.520 

1.525 

1.550 

1.535 

1.540 

1.545 

0.88729  30231 
0.88756  76278 
0.88786  25287 
0.88817  76586 
0.88851  29527 

-0.11957  99983 
-0.11927  05601 
-0.11893  83580 
-0.11858  34900 
-0.11820  60534 

0.05960  46439 
0.06416  73074 
0.06871  02697 
0.07323  36936 
0.07773  77400 

0.91451  37552 
0.91055  71245 
0.90663  32361 
0.90274  16984 
0.89888  21253 

0.525 

0.530 

0.535 

0.540 

0.545 

1.550 

1.555 

1.560 

1.565 

1.570 

0.88886  83478 
0.88924  37830 
0.88963  91990 
0.69005  45387 
0.89048  97463 

-0.11780  61446 
-0.11738  38595 
-0.11693  92928 
-0.11647  25388 
-0.11598  36908 

0.08222  25675 
0.08668  83334 
0.09113  51925 
0.09556  32984 
0.09997  28024 

0.89505  41371 
0.89125  73596 
0.88749  14249 
0.88375  59699 
0.88005  06378 

0.550 

0.555 

0.560 

0.565 

0.570 

1.575 

1.580 

1.585 

1.590 

1.596 

0.89094  47686 
0.89141  95537 
0.89191  40515 
0.89242  82141 
0.89296  19949 

-0.11547  28415 
-0.11494  00828 
-0.11438  55058 
-0.11380  92009 
-0.11321  12579 

0.10436  38544 
0.10873  66023 
0.11309  11923 
0.11742  77690 
0.12174  64754 

0.87637  50766 
0.87272  89402 
0.86911  18871 
0.86552  35815 
0.86196  36921 

0.575 

0.580 

0.585 

0.590 

0.595 

1.600 

1.606 

1.610 

1.615 

1.620 

0.89351  53493 
0.89408  82342 
0.89468  06085 
0.89529  24327 
0.89592  36685 

-0.11259  17657 
-0.11195  08127 
-0.11128  84864 
-0.11060  48737 
-0.10990  00610 

0.12604  74528 
0.13033  08407 
0.13459  67772 
0.13884  53988 
0.14307  68404 

0.85843  18931 
0.85492  78630 
0.85145  12856 
0.84800  18488 
0.84457  92455 

0.600 

0.605 

0.610 

0.615 

0.620 

1.625 

1.630 

1.635 

1.640 

1.645 

0.89657  42800 
0.89724  42326 
0.89793  34930 
0.89864  20302 
0.89936  98138 

-0.10917  41338 
-0.10842  71769 
-0.10765  92746 
-0.10687  05105 
-0.10606  09676 

0.14729  12354 
0.15148  87158 
0.15566  94120 
0.15983  34529 
0.16398  09660 

0.84118  31730 
0.83781  33330 
0.83446  94315 
0.83115  11790 
0.82785  82897 

0.625 

0.630 

0.635 

0.640 

0.645 

1.650 

1.656 

1.660 

1.665 

1.670 

0.90011  68163 
0.90088  30104 
0.90166  83712 
0.90247  28748 
0.90329  64995 

-0.10523  07282 
-0.10437  98739 
-0.10350  84860 
-0.10261  66447 
-0.10170  44301 

0.16811  20776 
0.17222  69122 
0.17632  55933 
0.18040  82427 
0.18447  49813 

0.82459  04826 
0.82134  74802 
0.81812  90092 
0.81493  48001 
0.81176  45875 

0.650 

0.655 

0.660 

0.665 

0.670 

1.675 

1.680 

1.685 

1.690 

1.695 

0.90413  92243 
0.90500  10302 
0.90588  18996 
0.90678  18160 
0.90770  07650 

-0.10077  19212 
-0.09981  91969 
-0.09884  63351 
-0.09785  34135 
-0.09684  05088 

0.18852  59282 
0.19256  12015 
0.19658  09180 
0.20058  51931 
0.20457  41410 

0.80861  81094 
0.80549  51079 
0.80239  53282 
0.79931  85198 
0.79626  44350 

0.675 

0.680 

0.685 

0.690 

0.695 

1.700 

1.705 

1.710 

1.715 

1.720 

0.9086  '  87329 
0.9095 3  57079 
0.91057  16796 
0.91196  6b390 
0.91258  05779 

-0.09580  76974 
-0.09475  50552 
-0.09368  26573 
-0.09259  05785 
-0.09147  88929 

0.20854  78749 
0.21250  65064 
0.21645  01462 
0.22037  89037 
0.22429  28871 

0.79323  28302 
0.79022  34645 
0.78723  61012 
0.78427  05060 
0.78132  64486 

0.700 

0.705 

0.710 

0.715 

0.720 

1.725 

1.730 

1.735 

1.740 

1.745 

0.913ol  34904 
0.91466  53712 
0.91573  62171 
0.91682  60252 
0.91793  47950 

-0.09034  76741 
-0.08919  69951 
-0.08802  69286 
-0.08683  75466 
-0.08562  89203 

0.22819  22037 
0.23207  69593 
0.23594  72589 
0.23980  32061 
0.24364  49038 

0.77840  37011 
0.77550  20396 
0.77262  12424 
0.76976  10915 
0.76692  13714 

0.725 

0.730 

0.735 

0.740 

0.745 

1.750 

0.91906  25268 

-0.08440  11210 

0.24747  24535 

0.76410  18699 

7“  iy  In  >/. 

0.750 

v! 

In  i/l 

In  ft- 

V 

r  n 

[•  ri  r  n 

loRu,  •  0.43429  41819 
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I. ,1.1.  (.1  (.  VMM  It|(.  VMM  \  \\l>  I  HM.  \  M  'I  \  11  M  IIUNS 


n  m  In  r( i )  J < < )  i  (r\ 


1.750 

0.91906 

25268 

-0.08440 

11210 

0.24747 

24535 

0.76410 

18699 

0.750 

1.755 

0.92020 

92224 

-0.08315 

42192 

0.25128 

59559 

0.76130 

23773 

0.755 

1.760 

0.92137 

48846 

-0.08188 

82847 

0.25508 

55103 

0.75852 

26870 

0.760 

1.765 

0.92255 

95178 

-0.08060 

33871 

0.25887 

12154 

0.75576 

25950 

0.765 

1.770 

0.92376 

31277 

-0.07929 

95955 

0.26264 

31686 

0.75302 

19003 

0.770 

1.775 

0.92498 

57211 

-0.07797 

69782 

0.26640 

14664 

0.75030 

04040 

0.775 

1.780 

0.92622 

73062 

-0.07663 

56034 

0.27014 

62043 

0.74759 

79107 

0.780 

1.785 

0.92748 

78926 

-0.07527 

55386 

0.27387 

74769 

0.74491 

422d8 

0.785 

1.790 

0.92876 

74904 

-0.07389 

68509 

0.27759 

53776 

0.74224 

91617 

0.790 

1.795 

0.93006 

61123 

-0.07249 

96070 

0.28129 

99992 

0.73960 

25271 

0.795 

1.800 

0.93138 

37710 

-0.07108 

38729 

0.28499 

14333 

0.73697 

41375 

0.800 

1.805 

0.93272 

04811 

-0.06964 

97145 

0.28866 

97707 

0.73436 

38093 

0.805 

1.810 

0.93407 

62585 

-0.06819 

71969 

0.29233 

51012 

0.73177 

13620 

0.810 

1.815 

0.93545 

11198 

-0.06672 

63850 

0.29598 

75138 

0.72919 

661 66 

0.815 

1.820 

0.93684 

50832 

-0.06523 

73431 

0.29962 

70966 

0.72663 

93972 

0.820 

1.825 

0.93825 

81682 

-0.06373 

01353 

0.30325 

39367 

0.72409 

95297 

0.825 

1.830 

0.93969 

03951 

-0.06220 

48248 

0.30686 

81205 

0.72157 

68426 

0.830 

1.835 

0.94114 

17859 

-0.06066 

14750 

0.31046 

97335 

0.71907 

11662 

0.835 

1.840 

0.94261 

23634 

-0.05910 

01483 

0.31405 

88602 

0.71658 

23333 

0.840 

1.845 

0.94410 

21519 

-0.05752 

09071 

0.31763 

55846 

0.71411 

01788 

0.845 

1.850 

0.94561 

11764 

-0.05592 

38130 

0.32119 

99895 

0.71165 

45396 

0.850 

1.855 

0.94713 

94637 

-0.05430 

89276 

0.32475 

21572 

0.70921 

52546 

0.855 

1.860 

0.94868 

70417 

-0.05267 

63117 

0.32829 

21691 

0.70679 

21650 

0.860 

1.865 

0.95025 

39389 

-0.05102 

60260 

0.33182 

01056 

0.70438 

51138 

0.865 

1.870 

0.95184 

01855 

-0.04935 

81307 

0.33533 

60467 

0.70199 

39461 

0.870 

1.875 

0.95344 

58127 

-0.04767 

26854 

0.33884 

00713 

0.69961 

85089 

0.875 

1.880 

0.95507 

08530 

-0.04596 

97497 

0.34233 

22577 

0.69725 

86512 

0.880 

1.885 

0.95671 

53398 

-0.04424 

93824 

0.34581 

26835 

0.69491 

42236 

0.885 

1.890 

0.95837 

93077 

-0.04251 

16423 

0.34928 

14255 

0.69258 

50790 

0.890 

1.895 

0.96006 

27927 

-0.04075 

65875 

0.35273 

85596 

0.69027 

10717 

0.895 

1.900 

0.96176 

58319 

-0.03898 

42759 

0.35618 

41612 

0.6B797 

20582 

0.900 

1.905 

0.96348 

84632 

-0.03719 

47650 

0.35961 

83049 

0.68568 

78965 

0.905 

1.910 

0.96523 

07261 

-0.03538 

81118 

0.36304 

10646 

0.68341 

84465 

0.910 

1.915 

0.96699 

26608 

-0.03356 

43732 

0.36645 

25136 

0.68116 

35696 

0.915 

1.920 

0.96877 

43090 

-0.03172 

36054 

0.36985 

27244 

0.67892 

31293 

0.920 

1.925 

0.97057 

57134 

-0.02986 

58646 

0.37324 

17688 

0.67669 

69903 

0.925 

1.930 

0.97239 

69178 

-0.02799 

12062 

0.37661 

97179 

0.67448 

50194 

0.930 

1.935 

0.97423 

79672 

-0.02609 

96858 

0.37998 

66424 

0.67228 

70846 

0.935 

1.940 

0.97609 

89075 

-0.02419 

13581 

0.38334 

26119 

0.67010 

30559 

0.940 

1.945 

0.97797 

97861 

-0.02226 

62778 

0.38668 

76959 

0.66793 

28044 

0.945 

1.950 

0.97988 

06513 

-0.02032 

44991 

0.39002 

19627 

0.66577 

62034 

0.950 

1.955 

0.98180 

15524 

-0.01836 

60761 

0.39334 

54805 

0.66363 

31270 

0.955 

1.960 

0.98374 

25404 

-0.01639 

10621 

0.39665 

83163 

0.66150 

34514 

0.960 

1.965 

0.98570 

36664 

-0.01439 

95106 

0.39996 

05371 

0.65938 

70538 

0.965 

1.970 

0.98768 

49838 

-0.01239 

14744 

0.40325 

22088 

0.65728 

38134 

0.970 

1.975 

0.98968 

65462 

-0.01036 

70060 

0.40653 

33970 

0.65519 

36104 

0.975 

1.980 

0.99170 

84087 

-0.00832 

61578 

0.40980 

41664 

0.65311 

63266 

0.980 

1.985 

0.99375 

06274 

-0.00626 

89816 

0.41306 

45816 

0.65105 

18450 

0.985 

1.990 

0.99581 

32598 

-0.00419 

55291 

0.41631 

47060 

0.64900 

00505 

0,990 

1.995 

0.99789 

63643 

-0.00210 

58516 

0.41955 

46030 

0.64696 

08286 

0.995 

2.000 

1.00000 

00000 

0.00000 

00000 

0.42278 

43351 

0.64493 

40668 

1.000 

,1 

,n 

1 

l\ 

In  >/! 

t  In  y! 

dtfl 

In  v ! 

V 

r(- 

6)2] 

n 

6)2] 

n 

6)2] 

?,21 

L 

4  J 

L 

4  J 

1 

4  J 

L 

4  J 

log.o  <-0.43429  44819 
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iin'iiiii.v  imi  i,i:i.atki>  kiwc tions 


<.  \  MM  \  \M»  UK.  VMM  \  II  V  IHt\»  I  Hit  I  \  lit.  I  It  \\l*  II  \l  II  \  I  I  ill!  \  \ll  |> 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 
52 

33 

34 

35 

36 

37 

38 

39 

40 

4! 

42 

43 

44 

45 


r(») 

(  0)1.00000  00000 
(  0)1.00000  00000 
0)  2.00000  00000 
(  0)6.00000  00000 
(  1)2.40000  00000 

(  2)1.20000  00000 
(  2)7.20000  00000 
(  3)5.04000  00000 
(  4)4.03200  00000 
(  5)3.62880  00000 

1  6)3.62880  00000 
7)3.99168  00000 
8)4.79001  60000 
9)6.22702  08000 
(10)8.71782  91200 

(  12)  1.30767  43680 
(  13)  2.09227  89888 
(14)  3.55687  42810 
(  15)  6.40237  37057 
(17)1.21645  10041 

(  18)2.43290  20082 
(19)5,10909  42172 
(21  1.12400  07278 
( 22) 2. 58520  16739 
(23)6.20448  40173 

(  25)  1.55112  10043 
(26)4.03291  46113 
(28)1.08888  69450 
(29)  3.04888  34461 
(  30)8.84176  19937 

(  32)2.65252  85981 
(  33)8.22283  86542 
(35)2.63130  83693 
(36)8.68331  76188 
(  38)2.95232  79904 

(40)  1.03331  47966 
(41)  3.71993  32679 
(  43)1.37637  53091 
(44)5.23022  61747 
(46)2.03978  82081 

(47)8.15915  28325 
(  49)  3.34525  26613 
(  51)  1.40500  61178 
(52)6.04152  63063 
(  54)2.65827  15748 


1/' 


1 

1:8 


(- 

<- 


(-13)  7.6 

(-14)4.7 
)  2.8 
) 1-5 
)  8.2 


-16 

-is; 


-20)1.9 
(-22)  8.8 
(-23)  3.8 
(-24)1.6 


3 

f: 


37 


(  * 


-50)2.98931 
-52)  7.11740 
-5M  1.65521 


46  ( 

56) 

11.19622  22087 

[  -57) 

47  1 

'57 

15.50262  21598 

!-58) 

48  ( 

|  59 

12.58623  24151 

-60) 

49  | 

!  61 

11.24139  15593 

-62 

50  | 

[  62 

16.08281  86403 

-63 

000 

(  - 1)8.86226 

93 

-0.57721  56649 

1.08443 

755 

0.57721 

566 

000 

0)  1.32934 

04 

<0.42278  43351 

1.04220 

712 

0.27036 

285 

000 

(  0  3.32335 

10 

0.92278  43351 

1.02806 

452 

0.17582 

795 

667 

(  1)  1.16317 

28 

1.25611  76684 

1.02100 

830 

0.13017 

669 

667 

(  1)  5.23427 

78 

1.50611  76684 

1.01678 

399 

0.10332 

02^ 

333 

(  2)2.87885 

28 

1.70611  76684 

1.01397 

285 

0.08564 

180 

889 

(  3)1.87125 

43 

1.87278  43351 

1.01196 

776 

0.07312 

581 

698 

4)  1.40344 

07 

2.01564  14780 

1.01046 

565 

0.06380 

006 

873 

(  5)1.19292 

46 

2.14064  14780 

1.00929 

843 

0.05658 

31  n 

192 

(  6)  1.1  3327 

84 

2.25175  25891 

1.00836 

536 

0.05083  250 

192 

(  7)  1.18994 

23 

2.35175  25891 

1.00760 

243 

0.04614 

268 

084 

(  8)1.36843 

37 

2.44266  16800 

1.00696 

700 

0.04224 

497 

570 

(  9)1.71054 

21 

2.52599  50133 

1.00642 

958 

0.03895 

434 

438 

(  10)2.30923 

18 

2.60291  80902 

1.00596 

9)1 

0.03613 

924 

456 

(  1  1)3.34838 

61 

2.67434  66617 

1.00557 

019 

0.03370 

354 

373 

(12)5.18999 

85 

2.74101  33283 

1.00522 

124 

0.03157 

539 

733 

(  1  3)  8.56349 

74 

2.80351  33283 

1.00491 

343 

0.02970  002 

725 

(  15)  1.49861 

21 

2.86233  68577 

1.00463 

988 

0.02803 

490 

070 

(  16)2.77243 

23 

2.91789  24133 

1.00439 

519 

0.02654 

6S7 

525 

(  17)  5.40624 

30 

2.97052  39922 

1.00417 

501 

0.02520 

828 

762 

(19)1.10827 

98 

3.02052  39922 

1.00397 

584 

0.02399  845 

411 

(20)2.38280 

16 

3.06814  30399 

1.00379 

480 

0.02289  941 

139 

(21)5.36130 

j  ■, 

3.11359  75853 

1.00362 

953 

0.02189 

663 

017 

(  23)1.25990 

63 

3.15707  58462 

1.00347 

806 

0.02097 

798 

757 

(  24)  3.08677 

05 

3.19874  25129 

1.00333 

872 

0.02013 

331 

029 

(25)  7.87126 

49 

3.23874  25129 

1.00321 

011 

0.01935  403 

626 

(27  2.08588 

52 

3.27720  40513 

1.00309 

105 

0.01863 

281 

986 

(28)  5.73618 

43 

3.31424  10884 

1.00298 

050 

0.01796 

342 

924 

(  30)1.63481 

25 

3.34995  53741 

1.00787 

758 

0.01734 

046 

629 

(  31)4.82269 

69 

3.38443  81327 

1.00278 

154 

0.01675 

925 

763 

(  33)  1.47092 

26 

3.41777  14660 

1.00269 

170 

0.01621 

574 

504 

(  34)4.63340 

61 

3.45002  95305 

1.00760 

748 

0.01570 

637 

076 

36  1.50585 

70 

3.48127  95305 

1.00252 

837 

0.01522 

803 

356 

(  37)  5.04462 

09 

3.51158  25608 

1.00245 

392 

0.01477 

796 

754 

<  39)  1.74039 

42 

3.54099  43255 

1.00238 

372 

0.01435 

374 

296 

(40)6.17839 

94 

3.56956  57541 

1.00231 

744 

0.01395 

318 

02  7 

(  42)2.25511 

58 

3.59734  35319 

1.00225 

474 

0.01357 

438 

018 

(  43)8.45668 

42 

3.62437  05589 

1.00219 

534 

0.01321 

560 

320 

(  45)  3.25582 

34 

3.65068  63484 

1.00213 

899 

0.01287 

530 

976 

(47)  1.28605 

02 

3.67632  73740 

1.00208 

546 

0.01255 

208 

744 

(  48)5.20850 

35 

3.70132  73740 

1.00203 

455 

0.01224 

469 

083 

(50)2.16152 

90 

3.72571  76179 

1.00198 

606 

0.01195 

200 

673 

51  9.18649 

81 

3.74952  71417 

1.00193 

983 

0.01167 

297 

087 

(  53)  3.99612 

67 

3.77278  29557 

1.00189 

570 

0.01140 

668 

288 

(  55)  1.77827 

64 

3.79551  02284 

1.00185 

354 

0.01115 

226 

084 

(  56)8.09115 

74 

3.81773  24506 

1.00181 

321 

0.01090 

895 

540 

(  58)  3.76238 

82 

3.83947  15811 

1.00177 

460 

0.01067 

602 

851 

(60  1.78713 

44 

3.86074  81768 

1.00173 

759 

0.01045 

283 

607 

(61)8.66760 

18 

3.88158  15102 

1.00170 

210 

0.01023 

879 

471 

(63)4.29046 

29 

3.90198  96734 

1.00166 

803 

0.01003 

333 

942 

(65)2.16668 

38 

3.92198  96734 

1.00163 

530 

0.00983 

596 

1! 

* 

an 

2.)!»" 

*<  ".Mm) 

*{") 

--ln»  /.t(») 

(2-)*-  2.50662  82746  31001 

51  (64)3.04140  93202  ( -65)  3.2E 

(»  D!  l/( 

a !  (2 *))»"  '  -  »  Mm)  r(< 

<H»l  compiled  from  H.  T.  Davis,  Tables  of  the  higher  mathematical  functions,  2  vols.  (Principia  Press, 
Bloomington,  Ind.,  1933,  1935)  (with  permission). 
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II 

r(» 

i 

l/i (»l 

mi  *■  5) 

<(») 

l'i('i) 

51 

(  64)  3.04140 

93202 

’-  65)  3.28794 

942 

(  65)2.16668 

38 

3.92198 

96734 

1.00163 

530 

0.00983 

596 

52 

(  66)1.55111 

87533 

-  67)6.44695 

964 

(  67)  1.11584 

21 

3.94159 

75166 

1.00160 

383 

0.00964 

620 

53 

(  67)  8.06581 

75171 

!-  68)  1.23979 

993 

(  68)5.85817 

12 

3.96082 

82858 

1.00157 

355 

0.00946 

363 

54 

(  69)4.27488 

32841 

!-  70)2.33924 

515 

(  70)  3.13412 

}6 

3.97969 

62103 

1.00154 

438 

0.00928 

784 

55 

(  71)2.30843 

69734 

[-  72)  4.33193 

547 

(  72)  1.70809 

63 

3.99821 

47288 

1.00151 

628 

0.00911 

846 

56 

(  73)  1.26964 

03354 

’  -  74)  7.87624 

631 

(  73)9.47993 

44 

4.01639 

65470 

1.00148 

919 

0.00895 

514 

57 

74)  7.10998 

58780 

’-  75)  1.40647 

255 

75)5.35616 

29 

4.03425 

36899 

1.00146 

304 

0.00879 

758 

58 

(  76)4.05269 

19505 

-  77)2.46749 

571 

(  7  7)  3.079  79 

37 

4.05179 

75495 

1.00143 

780 

0.00864 

546 

59 

(  78)  2.35056 

13313 

!-  79)  4.25430 

295 

(  79)  1.80167 

93 

4.06903 

89288 

1.00141 

341 

0.00849 

852 

60 

(  80)  1.38683 

11855 

[  -  81)  7.21068 

296 

(  81)1.07199 

92 

4.08598 

80814 

1.00138 

984 

0.00835 

648 

61 

(  81)8.32098 

71127  .  (-  82)  1.201  78 

049 

(  82)6.48559 

51 

4.10265 

47481 

1.00136 

704 

0.00821 

912 

62 

(  83)5.07580 

21388  1 

’-  84)  1.97013 

196 

(  84)  3.98864 

10 

4.11904 

81907 

1.00134 

498 

0.00808 

619 

63 

(  85)  3.14699 

73260  | 

-  86)  3.17763 

219 

(  86)2.49290 

06 

4.13517 

72229 

1.00132 

?62 

0.00795 

750 

64 

(  87)1.98260 

83154  (-  88)  5.04386 

062 

(  88)1.58299 

19 

4.15105 

02388 

1.00130 

292 

0.00783 

284 

65 

(  89)  1.26886 

93219  | 

[-  90)  7.88103 

221 

(  90)  1.02102 

98 

4.16667 

52388 

1.00128 

286 

0.00771 

203 

66 

(  90)8.24765 

05921  ( 

'  -  91)  1.21246 

649 

(  91)6.68774 

50 

4.18205 

98542 

1.00126 

341 

0.00759 

489 

67 

(  92)  5.44344 

93908  ( 

1  -  93)  1.83707 

044 

(  93)  4.44735 

04 

4.19721 

13693 

1.00124 

455 

0.00748 

125 

68 

(  94)  3.6471  1 

10918  ( 

'-  95)  2.74189 

619 

95)3.00196 

15 

4.21213 

67425 

1.00122 

62  3 

0.00737 

096 

64 

(  96)2.48003 

55424  ( 

-  97)4.03220 

028 

(  97)2.05634 

36 

4.22684 

26248 

1.00120 

845 

0.00726 

388 

70 

(  98)1.71122 

45243  I 

:  -  99)  5.84376 

852 

(  99)1.42915 

88 

4.24133 

53785 

1.00119 

118 

0.00715 

986 

71 

(  100)  1.19785 

71670  ( 

'  -101)  8.34824 

074 

(  101)1.00755 

70 

4.25562 

10927 

1.00117 

439 

0.00705 

878 

72 

(  101)8.50478 

58857  ( 

i -102)  1.1  7580 

85o 

(  102)  7.20403 

24 

4.26970 

55998 

1.00115 

807 

0.00696 

052 

73 

(  103)6.12344 

58377  (-104)  1.63306 

744 

104)5.22292 

35 

4.28359 

44887 

1.00114 

220 

0.00686 

495 

74 

(  105)  4.47011 

54615  ( 

;  -106)2.23707 

868 

(  106)  3.83884 

87 

4.29729 

31188 

1.00112 

675 

0.00677 

197 

7C 

(  107)  3.30788 

54415  | 

i  -108)  3.02307 

930 

(  108)2.85994 

23 

4.31080 

66323 

1.00111 

172 

0.00668 

148 

76 

(  109)  2.48091 

40811  ( 

-110)  4.03077 

240 

(  110)2.15925 

64 

4.32413 

99657 

1.00109 

709 

0.00659 

337 

77 

(  111)  1.88549 

4701  7  ( 

‘-112)5.30364 

789 

(  112)1.65183 

12 

4.33729 

78604 

1.00108 

283 

0.00650 

756 

78 

(113)1.45183 

09203  ( 

-114  6.88785 

441 

(  114)  1.28016 

92 

4.35028 

48734 

1.00106 

894 

0.00642 

395 

79 

(  11  5)  1.13242 

81178  ( 

1  -116)8.83058 

257 

(116)1.00493 

28 

4.36310 

53862 

1.00105 

540 

0.00634 

247 

80 

(  116)8.94618 

21308  ( 

\  -117)  1.11779 

526 

(  117)  7.98921 

57 

4.37576 

36140 

1.00104 

220 

0.00626 

302 

81 

(  118)  7.15694 

^7046  ( 

-119)  1.39724 

408 

(  119)6.43131 

87 

4.38826 

36140 

1.00102 

933 

0.00618 

554 

82 

(  120)  5.79712 

60207  ( 

-121  1.72499 

269 

121)5.24152 

47 

4.40060 

92931 

1.00101 

677 

0.00610 

995 

83 

(  122)4.75364 

33370  ( 

1  -123)  2.10364 

962 

(  123)4.32425 

79 

4.41280 

44150 

1.00100 

452 

0.00603 

619 

84 

(  124)  3.94552 

39697  ( 

-125)  2.53451 

761 

125)  3.61075 

53 

4.42485 

26078 

1.00099 

255 

0.00596 

419 

85 

(  126)  3.31424 

01346  ( 

-127)  3.01728 

287 

(  127)  3.05108 

83 

4.43675 

73697 

1.00098 

087 

0.00589 

389 

86 

(128)2.81710 

41144  | 

I  -129)  3.54974 

456 

(  129)2.60868 

05 

4.44852 

20756 

1.00096 

946 

0.00582 

522 

37 

(  130)  2.42270 

95384  ( 

-131)  4.12760 

995 

(  131)2.25650 

86 

4.46014 

99825 

1.00095 

831 

0.00575 

814 

88 

(  132)2.10775 

72984  ( 

1  -133)  4.74437 

926 

(  133)  1.97444 

50 

4.47164 

42354 

1.00094 

741 

0.00569 

258 

89 

(  134)  1.85482 

64226  | 

-135)  5.39134 

006 

1  35)  1.74738 

38 

4.48300 

78718 

1.00093 

676 

0.00562 

850 

90 

(  136)  1.65079 

55161  I 

I  -137)  6.05768 

546 

(  137)  1.56390 

85 

4.49424 

38268 

1.00092 

635 

0.00556 

584 

91 

(  138)  1.48571 

59645  I 

[  -139)  6.73076 

163 

(  139)  1.41533 

72 

4.50535 

49379 

1.00091 

617 

0.00550 

457 

92 

(  140)  1.35200 

15277  I 

-141  7.39644 

134 

(  141)  1.29503 

36 

4.51634 

39489 

1.00090 

620 

0.00544 

463 

93 

( 142) 1.24384 

14055  ( 

’  -143)8.03961 

016 

(  143)  1.19790 

60 

4.52721 

35142 

1.00089 

646 

0.00538 

598 

94 

(  144)  1.15677 

25071  ( 

’  -145  8.64474 

211 

(  145)  1.12004 

22 

4.53796 

62023 

1.00088 

691 

0.00532 

858 

95 

(  146)  1.08736 

61567  1 

[  -147)9.19653 

415 

(  147)  1.05843 

98 

4.54860 

45002 

1.00087 

757 

0.00527 

239 

96 

(  148)  1.03299 

78488  1 

[ -149)9.68056 

227 

(  149)1.01081 

00 

4.55913 

08160 

1.00086 

843 

0.00521 

738 

97 

(  149)9.91677 

93487  1 

I  -150)  1.00839 

190 

150  9.75431 

69 

4.56954 

74827 

1.00085 

94  7 

0.00516 

350 

98 

( 151)9.61927 

59682  1 

[  -152)  1.03957 

928 

152)9.51045 

90 

4.57985 

67610 

1.00085 

070 

0.00511 

072 

99 

(  153)9.42689 

04489  1 

(-154)  1.06079 

519 

( 154  9.36780 

21 

4.59006 

08426 

1.00084 

210 

0.00505 

901 

100 

(  155)9.33262 

15444  1 

(-156)  1.07151 

029 

(  156)9.32096 

31 

4.60016 

18527 

1.00083 

368 

0.00500 

833 

101 

(  157)9.33262 

15444  (-158)  1.07151 

029 

(  158)9.36756 

79 

4.61016 

18527 

1.00082 

542 

0.00495 

866 
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CALCULATOR  METHOD 


The  expression  f  (a  +  1)  =  a  F  (a)  is  applied  repeatedly  to  increase  the 
value  of  the  argument  until  is  it  greater  than  9,  when  Stirling's 
formula  can  be  applied.  If  the  argument  is  initially  greater  than  9, 
Stirling's  formula  is  used  at  once. 

Stirling's  formula  is  given  in  Abramowitz  w 

For  a  >  9  : 

DS  V  =  W  ■+  I 

X  ( -+•  i  C\. 

5  --  0.  ?  A9rz5352  -  Y  ^  (v  +-X5J  Jlr.Y 

_ ; _ r  __  ,  f 

C:2  v'  ;  !  *  50  Y 2  +  05  Y4 
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EXAMPLE 


FIND:  T  (10.4) 

SOLUTION:  Y  =  11.4  ;  X  =  (11.4)  (10.4) 

S  -  0.9189385332  -  11.4  +  (11.9)  Sin  11.4 

i  r ,  _  i  +  ' 

OzXh.4)  30(j|.4)2  I05(ll.4)4 

— 

S  =  18.4  8624553 
e5=  1. 06776  191 1  *  10 
X  =  I  I  8.56 

5 

r0o.4)=  -y-  =  900608.9 

For.  o<  <  9  * 

de  fimf  Y=  9  fractional  I  part  of  cx 

X  —  C<  C  cx  4-  |  )(cX  -f-  £  )  *  ...  Y 

5  =  O.  9189385332  -Y  ^(y+0.5)A  Y 
I  T,  _  _J _  I 

+  (IZ  Y)  3 0  YZ  150  Y4 

s 

e 

X 


V  (erf) 


I 


<3 


♦4 


EXAMPLE 


FIND:  T  (2.4) 


SOLUTION: 


Y  =  9.4 


X  =  9.4  x  8.4  x  7.4  x  6.4  x  5.4  x  4.4  x  3.4  x  2.4 


X  =  725029.0842 


S  =  0.9189385332  -  9.4  +  9.9  *n  9.4 


(12.  )(  9.4) 


5  =  13.  71052  657 

rM=  C  = 


150  (9.4)' 


9.  00608902/  *10 
7.  25  02  90342.  *  lo' 


(2.  4)  —  1.24216954-6 


As  A  CHEO-T 


P  (P„4)  =1.4-  P  0.4^  =■ 


P  (2.4)=  I.  242169345 


1.4  *  0.  8872  6  38175 

See  Table  J 
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BASIC  LANGUAGE  PROGRAM 


The  previous  calculator  method  is  shown  as  a  program  listing  in 
Basic  Language,  PL-1. 


PL-1 

PROGRAM  LISTING  FOR  T  V'  IN  BASIC 


3  0  RIM  A-  ALPHA 

!'7'  M  G  '  GAMMA  FUNCTION  Will  I  ARGUN!  NT  Al  F'HA 
30  A~5 .  6*1 
•1 0  Y  -  A 
50  X  =A 
60  Y  r  Y+ 1 
70  X  =  X  *  Y 
80  D  =  Y  - 

90  IF  D > ~ 0  THIN  110 

1 00  GO  TO  60 

1 10  8- .918938533 2 

1 2  0  S  ~  8  4  <  Y  +  .  5  >  *  L  0  G  ( Y  )  -  Y 

'.!  30  V-  1-  ( 1  /  <30  *  Y-i-2 )  )  4  ( 1  /  ( 105*  Y  M  )  ) 

140  V= < 1 / ( 12* Y) ) #V 
150  8-8 i V 
3  60  G  -  E  X  P  ( 8 )  /  X 
370  PRINT  ArG 
180  END 
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TI-59  Methods 


A.  USER  ENTERED  PROGRAM 

The  following  program  PL-2  is  for  use  with  the  TI-59  Programmable  Calculator 
and  is  part  of  Texas  Instrument’s  Math  39  program  exchange.  It  calculates 
the  Gamma  Function  T(x)  for  integer  and  non-integer  values  of  the 
ag rumen t  x. 


An  internal  accuracy  check  using  a  fractional  argument  can  be  made  by 
calculating  1(0.5),  then  squaring  the  answer.  The  result  should  be  tt. 
Subtracting  the  stored  value  of  it  from  the  previously  calculated  value 
gives  4x10  Thus,  the  error  is  4  in  the  11th  digit.  Other  non-integer 

values  have  been  compared  with  the  National  Bureau  of  Standards  Tables  rln 
and  almost  all  agree  within  +/-  one  unit  in  the  10th  digit.  PL-2  accuracy 
is  considered  accurate  for  all  practical  purposes. 


To  use  enter  the  value  of  x;  then  press 
will  ultimately  be  displayed, 
x  <  69.5 


.  The  computed  value  of  T(x) 


Calculation  time  is  5  to  15  seconds,  depending  on  x. 


EXAMPLE:  Compute  T(6.5) 

Enter  6.5.  PRESS  jlTj  . 

The  displayed  output  is  287.8852778 
Thus,  T(6.5)  =  287.8852778 


I 
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TI-59  Methods  (Cont’d) 


B.  MATH  UTILITY  MODULE  PROGRAM  MU-11 


The  new  math  utility  module  can  be  used  for  computing  T(x)  directly. 


EXAMPLE:  Compute  T (6 - 3) 


ENTER 


PRESS 


[2nd]  [Pgm]  11 


DISPLAY 


201.8132752 


Thus,  T (6. 3)  =  201.8132752 


GAMMA  FUNCTION 


LISTING  FOR  TI-59  PROGRAMMABLE  CALCULATOR 


r_  -  / 


LOC 

CODE  KEY 

000 

7b  LBL 

n  n  i 

14  U 

fj  f“|  J 

42  STG 

ij  |J . 

0  0  *■+ 

42  810 

nns 

f ;  1  \ 

111!  f, 

b8  HOP 

007 

b9  OP 

009 

43  RCL 

o  i  o 

00  00 

on 

b5 

012 

43  RCL 

H 1 :-! 

0 1  0  1 

014 

0 1  5 

42  STD 

0 1  b 

00  00 

0  1  7 

43  ROL 

0 1  8 

0 1  0 1 

019 

“  CT  _ 

i  J 

020 

MS  9 

U  9  I 

95  = 

0  2  2 

77  GE 

023 

00  00 

024 

99  99 

025 

b  1  GTO 

0  2  6 

0  0  0  U 

027 

07  07 

028 

bS  NOP 

029 

53  < 

030 

53  < 

031 

02  2 

032 

b5 

033 

o  9  -fi 

034 

54  ) 

035 

23  LHM 

03b 

05 

037 

9 9 

038 

05  5 

039 

85  + 

040 

53  < 

041 

43  RCL 

042 

0 1  0 1 

043 

85  + 

044 

9  9 

045 

05  5 

04b 

54  > 

LOC 

CODE  KEY 

047 

85  x 

048 

4  3  H  0  L 

0  4  9 

0  1  0  1 

050 

23  LNM 

051 

7c, 

052 

43  RCL 

j~,  cr  j 

0 1  0 1 

054 

35  + 

055 

43  RCL 

05b 

0 1  0 1 

n  cr  "7 

It.  i  : 

•j  cr.  i  ..  «v« 

058 

cr  cj  _ 

059 

U  i  i 

ObO 

02  2 

Ob  1 

85  x 

0  6  iii 

53  < 

Ub3 

U I  1 

0  b  4 

75  - 

n  b  5 

5  8  f 

0  b  b 

43  RCL 

0  b  7 

01  01 

0  b  8 

: K  ± 

0  b  9 

85 

070 

0  3  3 

071 

00  0 

072 

54  } 

0  7  3 

35  1  ■■■"  X 

074 

85  + 

075 

58  (' 

0  7  b 

43  RCL 

077 

0 1  0 1 

U  i  o 

079 

y  2; 

0  8  0 

85  X 

081 

01  1 

082 

00  0 

083 

05  5 

084 

54  > 

085 

35  1  'N 

0  8  b 

54  > 

087 

54  ) 

0  8  8 

22  INV 

089 

23  LHN 

090 

55 

091 

43  RCL 

092 

00  00 

093 

95  = 

094 

91  R--S 

15 


HP-67  Method 


A  program  listing  is  shown  (PL- 3) for  computing  T(a)  on  an  HP-67  Programmable 
Calculator  using  the  previously  described  Stirling's  approximation 
formula. 
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i. 


iffipmrw*  «■*'*  *  TW*r  ■  :■****■ 


1 


I 


PL  -  3  C Coh,T  >p^ 

Ef.IR*  KEY  COLfE  LM5  S»£P  KEY  ENTRY  KEY  CODE  COMMENTS 


HP-34C  METHOD 


The  key  x!  can  be  used  directly  to  calculate  T(x)  using  the 
following  relationship: 

r  (x)  =  (x  -  1)  : 

EXAMPLE:  Calculate  T(6.3) 

KEYSTROKES  DISPLAY  COMMENT 

6.3  6.3  x 

ENTER  +  1[-]  5.3  x-1 

[h]  [x!]  201.8132  T(6.3) 


I 


* 

ii 

g 


B.  INCOMPLETE  GAMMA  FUNCTIONS  COMPUTATION 


Def init ion: 


The  two  Incomplete  Gamma  Functions  y(a,r)  and  Q  (u,t)  are  defined 


as  follows: 


Y(a,T)  = 


/r 

C  o-l 

y x 


e  dx 


Q(a»  t) 


oO 

r  a- 1  -x 

x  e  dx 


The  Incomplete  Gamma  Functions  are  related  to  the  Complete  Gamma 
Function  f  (a)  as  follows: 


f(a)  =  Y  (a.t)  +  Q  (cx,t) 


That  is 


v  dx  = 


e 

*  J 


a-1  -x 

where  v  =  x  e 


1 IT*  '"*2* 


TABULAR  METHOD 


On  page  941  of  Abramowitz 

Y(a,r)  =  r(a)  P  (x2|v) 

2 

where  x  =  2i  ;  v  =  2a 

Q  (a,x)  =  T (a)  -  y(a,T) 

2, 

Also  Q  (a, t)  =  r(a)  Q  (x  |  v) 

NOTE:  Q  (a,i)  4  Q  (  x^v) 

P  ( X 2  I v )  =  1  -  Q  (x2|v) 

2 

Values  of  Q  (x  | v)  are  tabulated  in  Table  III,  on  pages  978-983 
a-1  -x 

where  y  =  x  e 

A  pictorial  representation  is  shown  in  figure  2.  Note  that  v 

m3X 

at  x  =  a-1. 


[i] 

occurs 
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EXAMPLE 


FIND:  y(1.5,0.9),  Q  (1. 5,0.9) 

SOLUTION:  ct=  1.5;  t  =  0.9 

T(1.5)  =  0.8862269255  from  Table  I  (Table  6.1  page 
268  flj  ) 

=  2t  =  1.8;  v  =  2a  =  3 
In  Table  26.7  page  978  [1] 

Q(x2|v)  «  0.61493 
P  ( X 2  1  v)  =  1  -  Q  (x2  !  v)  =  0.38507 
Y(1 . 5,0.9)  =  F ( 1 . 5)  P(x2|v  )  =  0.34126 
Q(l.  5,0.9)  =  r(l .  5)  0(x2^  )  =  0.54497 

As  a  check 

y(1.5,0.9)  +  Q(l. 5,0.9)  =  0.88623  to  five  places  =  T(1.5) 
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Interpolation  formulas  are  shown  on  the  bottom  of  p.  981  j_l__  (Table  III). 

2 

One  formula  is  for  interpolation  on  x  alone.  The  Double  Entrv  Enter- 

2 

polation  formula  is  for  interpolation  on  both  x  and  v. 

2  2 

where  4  =  h  (x  -  X  ) 

o 

w  =  v  -  v0  '  o 

For  interpolation  on  v  alone,  the  Double  Entrv  formula  can  be  modified 

2  2 

bv  letting  4=0  (i.e.  x  =  X  )  as  follows: 

o 

2,  2,  r  2  w' 

Q  (X  |v)  =  Q  (x  V  -1)  !'iw  -  7 
oo  L  z  J 

+  Q(x  |  v  )  !  l  -  w ^  i 

0  0-  — 

,  2 .  ;  2  w"' 

+  Q  (x0  I  v0  +1>  :  3  w  +  2 

EXAMPLE: 

FIND:  Y(5 . 64 ,  8)  ,  Q  (5.64,  8) 

SOLUTION: 

a  =  5.64  ;  t  =  8 

T(5.64)  =  4.64  x  3.64  x  2.64  x  1.64  x  F(1.64) 

r(1.64)  =  0.8986420302  from  Table  I  (Table  6.1  page  269  [_l]  ) 

K5.64)  =  65.71338911 
2  2 

XQ  =  X  =  2x  =  16  ;  v  =  2a  =  11.28 

In  this  example 

v  =  11  ;  v  -1  =  10  ;  v„  +1  =  12 

o  o  0 

w  =  v  -  v  =  11.28  -  11  =  0.28 
0 


Thus 


2  - 
Q(x  |v)  =  (0.09963)  _1_28  -  ^ 

2  2 

i 

+  (0.14113)  ! _ 1  -  .28 

+  (0.19124)  I  ,  282  +  _^28  j 

2  2 

=  -  0.01004  +  0.13006  +  0.03427 

2 

Q(x  h)  =  0.15429 

P(x2fv)  =  1  -  Q  (x2|v)  =  0.84571 

y(5.64,  8)  =  T (5.64)  P(x2|v)  =  55.574 

2 

Q(5.64,  8)  =  (5.64)  Q(x  |v)  =  10.139 

As  a  check 

'y(5.64,  8)  +  Q(  5 . 64 ,  8)  =  65.713  to  five  places  =  r(5.64) 
EXAMPLE 

FIND:  Y (4. 3,  3.77)  ,  Q  (4.3,  3.77) 

SOLUTION: 

a  =  4.3  ;  t  =  3,77 

T  (4.3)  =  3.3  x  2.3  x  1.3  x  T( 1.3) 

r  (4.3)  =  8.855343359 
2 

X  -  2t  =  7.54  ;  v  =  2a  =  8.6 

For  this  example 

2  2 
X  =  7.54  ;  x0  =  7.4 

•<!>  =  ‘5  (7.54  -  7.4)  =  0.07 

v  =  8.6  ;  v  =8  ;  v  +1=9 

o  o 

v  -4  =  4 

o 

VQ  -  2  =  6 

v  -1  =  7 

0 
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EXAMPLE  (Cont'd) 


w=8.6-8=0.6 


The  Double  Entry  Interpolation  expression  on  page  981  (Table  III) 


will  be  used. 


Q(XoN)  =  Q<X„  I'-'  -4) 


1  2 


+Q(x2No  ~2)  I  <(>  -  4>2  ~  w  4>  J 

o  •- 

fl  2  I 

+Q(x0|v0  -1)  2  w  -  2  w  +  w<!> 


v  r,  2  i 

)  l-w  -  4>+‘j 


+Q(x0h0 


l  2  i 

-hKXqNo  +  ^  2™  +  2w-w‘*> 


Q(x2|v)  =  (0.11620) j  -  (.07) 2  I 


+  (0.28543)  i  .07  -  .07  -  (.6)(.07) 

+  (0.38845)  \  (.6)2  -  \  (.6)  +  (.6)(.07) 

+  (0.49415)  i"l  -  •  62  -  .07  +y  (.07)2  +  (.6)(.07)_ 

r  i  2i  •• 

+  (0.59555)|  2  (.6)  +  J  (•&)  ~  ( - 6) ( .07) 


Q(xz|v)  =  0.00028469  +  0.006593433 

-0.0302991  +  0.30363046  +  0.2608509 


0.5410603905 


Q(4. 3,  3.77)  =  T(4.3)  Q(x  |v) 

=  4.7913  to  five  places 


4.791276  to  seven  places 


y  (4.3,  3.77)  =  r(4.3) 


[l  -  QCx2 I v)J 


=  4.0641  to  five  places 
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4.06406 8  to  seven  places 
8.85534  to  six  places 


TABLE  III  TABULATION  OF  Q(y2|v) 


M.'i'l: \MI.1 1  V  1  IMTIOXS 


i  I.!..  2«,.t  I'litmutinn  imi  i.iki.  oi  >  nisTitiiu  tiov  im  omi-i  i:it.«;\mm\ h  notion 

(i  nii  vmi  sims  ok  mi:  imhsson  ihmisiiu  tion 


•  0.001 

0.002 

0.003 

0.004 

0.005 

0.006 

0.007 

0.008 

0.009 

0.010 

' 

m  0.0005 

0.0010 

0.0015 

0.0020 

0.0025 

0.0030 

0.0035 

0.0040 

0.0045 

0.0050 

1 

0.97477 

0.96433 

0.95632 

0.  94957 

0.94363 

0.  93826 

0.93332 

0.92873 

0.  92442 

0.  92034 

2 

0.99950 

0.  99900 

0.99850 

0.  99800 

0.99750 

0.99700 

0.  99651 

0.99601 

0.99551 

0.99501 

3 

0.99999 

0.  99998 

0.  99996 

0.  99993 

0.  99991 

0.99988 

0.  99984 

0.  999B1 

0.  99977 

0.99973 

4 

0.99999 

0.99999 

0.  99999 

0.99999 

x-  0.01 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0.09 

0.10 

1 

m  0.005 

0.010 

0.015 

0.020 

0.025 

0.030 

0.035 

0.040 

0.045 

0.050 

1 

0.92034 

0.  88754 

0.86249 

0.  84148 

0.  82306 

0.80650 

0.79134 

0.  77730 

0.  76418 

0.75183 

2 

0.99501 

0.99005 

0.98511 

0.  98020 

0.  97531 

0.  97045 

0.96561 

0. 96079 

0.95600 

0.95123 

3 

0.99973 

0.  99925 

0.99863 

0.99790 

0.  99707 

0.99616 

0.99518 

0.99412 

0.99301 

0.99184 

4 

0.99999 

0.99995 

0.99989 

0.99980 

•0.99969 

0.  99956 

0.  99940 

0.99922 

0.99902 

0. 99879 

5 

0. 99999 

0.99998 

0.99997 

0.  99995 

0.  99993 

0. 99991 

0.  99987 

0. 99984 

6 

0.  99999 

0.  99999 

0.  99999 

0.  99998 

»'->  0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1 

m  0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

1 

0.  75183 

0.65472 

0.  58388 

0.  52709 

0.47950 

0.43858 

0.40278 

0. 37109 

0.  34278 

0.31731 

2 

0.95123 

0.90484 

0.86071 

0.81873 

0.  77880 

0.  74082 

0.  70469 

0.  67032 

0.63763 

0. 60653 

3 

0.99184 

0.  97759 

0.960C3 

0.94024 

0.91889 

0.  89643 

0.87320 

0. 84947 

0.82543 

0.80125 

4 

0.99879 

0.99532 

0.98981 

0.98248 

0.97350 

0.96306 

0.95133 

0.93845 

0.92456 

0.90980 

5 

0.99984 

0.99911 

0.99764 

0.99533 

0.99212 

0.98800 

0.98297 

0.97703 

0.97022 

0.96257 

fa 

0.99998 

0.99985 

0.99950 

0.99885 

0.99784 

0.99640 

0.99449 

0.99207 

0.98912 

0.98561 

7 

0.99997 

0.99990 

0.99974 

0.99945 

0.99899 

0.  99834 

0.  99744 

0. 99628 

0. 99483 

8 

0.99998 

0.  99994 

0.99987 

0.99973 

0.99953 

0.99922 

0.99800 

0.99825 

9 

0.  99999 

0.  99997 

0.  99993 

0.  99987 

0.  99978 

0. 99964 

0.99944 

10 

0.99999 

0.99998 

0.99997 

0.99994 

0. 99989 

0.99983 

n 

0.99999 

0.99998 

0.99997 

0.99995 

12 

0. 99999 

0.99999 

x-  LI 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2.0 

1 

m.  0.55 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0.95 

1.00 

1 

0.29427 

0.27332 

0.  25421 

0.  23672 

0.22067 

0.20590 

0.19229 

0.17971 

0.16808 

0.15730 

2 

0.  57fa95 

0.  54881 

0. 52205 

0.49659 

0.47237 

0.44933 

0.42741 

0.40657 

0.  38674 

0. 36788 

3 

0.  77707 

0.75300 

0.  72913 

0.  70553 

0.  68227 

0.65939 

0.  63693 

0.61493 

0.  59342 

0.  57241 

4 

0.89427 

0.87810 

0.  86138 

0.84420 

0.  82664 

0. 80879 

0.  79072 

0.  77248 

0.  75414 

0.  73576 

5 

0.95410 

0.94480 

0.93493 

0.92431 

0.91307 

0.90125 

0.  88890 

0.87607 

0.86280 

0.84915 

fa 

0.98154 

0.97689 

0.97166 

0.  96586 

0.95949 

0.95258 

0.94512 

0.93714 

0.92866 

0.91970 

7 

0.99305 

0.  99093 

0.  98844 

0.98557 

0.98231 

0.97864 

0.97457 

0.97008 

0.96517 

0.95984 

8 

0.99753 

0.99664 

0.99555 

0.99425 

0.99271 

0.99092 

0.  98887 

0.  98654 

0.98393 

0.98101 

9 

0.99917 

0.99882 

0.99838 

0.99782 

0.99715 

0.99633 

0.99537 

0.99425 

0.99295 

0.99147 

10 

0.99973 

0.99961 

0.99944 

0.99921 

0.  99894 

0.  99859 

0.99817 

0.99766 

0.99705 

0.99634 

11 

0.99992 

0.  99987 

0.99981 

0.99973 

0.  99962 

0.99948 

0.99930 

0.  99908 

0. 99882 

0. 99850 

12 

0.99998 

0.99996 

0.  99994 

0.  99991 

0.  99987 

0.  99982 

0.99975 

0.  99966 

0.  99954 

0. 99941 

13 

0.  99999 

0.  99999 

0.  99998 

0.  99997 

0.  99996 

0.  99994 

0.  99991 

0.  99988 

0.  99983 

0.99977 

14 

0.  99999 

0.99999 

0.  99999 

0.  99998 

0.99997 

0.  99996 

0. 99994 

0. 99992 

li 

0.  99999 

0.  99999 

0.  99999 

0.  99998 

0. 99997 

lb 

0. 99999 

0. 99999 

v< 

•-M  l  /’(» 

*•>■[*  r 

©]  '/■> 

i 

V  ill  - 

m  j 

f  ”  r  '/• 
'*»2 

*>«» 

II 

S.  *» 

1 

’"mJ  '_/’!(►  i 

even,  e-i*. 

m-  lx2) 

Compiled  from  E.  S.  Pearson  and  H.  0.  Hartley  (editors),  Biometrika  tables  for  statisticians,  vol.  1. 
Cambridge  Univ.  Press,  Cambridge,  England,  1954  (with  permission). 
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TABLE  III  (Cont'd) 


PROBABILITY  FUNCTIONS 


PKOltMtll.lTY  IMKUtll.  OK  x--l>ISTKIItlTIMN.  INCOMIM.K.TK  GAMMA  FUNCTION  Table  26.7 
Cl  Ml  I.VTIYK  SIMS  OK  TIIK  POISSON  UlSTKIlilTION 


2.2 

2.4 

2.6 

2.8 

3.0 

3.2 

3.4 

3.6 

3.8 

4.0 

9 

m  1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2.0 

1 

0.  13801 

0. 12134 

0.  10686 

0. 09426 

0.  08327 

0.  07364 

0.  06520 

0.05778 

0.  05125 

0.  04550 

2 

0.  33287 

0.  30119 

0.  27253 

0.  24660 

0.  22313 

0.  20190 

0.18268 

0. 16530 

0.  14957 

0.13534 

3 

0.  53195 

0.  49363 

0.45749 

0. 42350 

0.  39163 

0.  36181 

0.33397 

0.30802 

0.28389 

0.26146 

4 

0. 69903 

0.  66263 

0.62682 

0. 59183 

0.  55783 

0.  52493 

0.  49325 

0.  46284 

0.  43375 

0.  40601 

5 

0.82084 

0.  79147 

0. 76137 

0. 73079 

0.  69999 

0.  66918 

0.63857 

0.60831 

0.  57856 

0.54942 

6 

0. 90042 

0.  87949 

0.85711 

0.83350 

0.  80885 

0.  78336 

0.  75722 

0.  73062 

0.  70372 

0.  67668 

7 

0.  94795 

0. 93444 

0.91938 

0.90287 

0.  88500 

0. 86590 

0.  84570 

0.  82452 

0.  80250 

0.  77978 

8 

0.  97426 

0.96623 

0.95691 

0.  94628 

0.  93436 

0.  92119 

0.  90681 

0.89129 

0,87470 

0.85712 

9 

0.  98790 

0.  98345 

0.97807 

0.97170 

0. 96430 

0.  95583 

0.94631 

0.93572 

0.  92408 

0.91141 

10 

0.  99457 

0.  99225 

0.98934 

0.  98575 

0.  98142 

0.  97632 

0.97039 

0.96359 

0.95592 

0.94735 

n 

0.99766 

0.  99652 

0.  99503 

0.99311 

0.99073 

0.98781 

0.98431 

0.98019 

0.  97541 

0.  96992 

12 

0.  99903 

0.  99850 

0.99777 

0.99680 

0.  99554 

0.  99396 

0.99200 

0.98962 

0.  98678 

0. 98344 

13 

0.  99961 

0.  999  38 

0.  99903 

0.  99856 

0.  99793 

0.99711 

0.  99606 

0.99475 

0. 99314 

0.99119 

14 

0.  99985 

0.  999  75 

0.99960 

0.99938 

0.99907 

0.99866 

0.99813 

0.99743 

0.  99655 

0.99547 

15 

0.99994 

0.  99990 

0.  99984 

0.  99974 

0.  99960 

0.  99940 

0. 99913 

0.  99878 

0.  99832 

0.  99774 

16 

0.99998 

0. 99996 

0.99994 

0.  99989 

0.  99983 

0.  99974 

0.  99961 

0.99944 

0. 99921 

0.99890 

17 

0.  99999 

0.  99999 

0. 99998 

0.  99996 

0.  99993 

0.  99989 

0.  99983 

0.  99975 

0.  99964 

0.  99948 

18 

0. 99999 

0.99998 

0.99997 

0.99995 

0.  99993 

0.  99989 

0.99984 

0.99976 

19 

0.99999 

0.99999 

0.  99998 

0.  99997 

0.  99995 

0.  99993 

0.99989 

20 

0. 99999 

0.  99999 

0.99998 

0. 99997 

0.99995 

21 

0.  99999 

0.  99999 

0.  99998 

22 

0.99999 

4.2 

4.1 

4.6 

4.8 

5.0 

5.2 

5.4 

5.6 

5.8 

6.0 

r 

m  2.1 

2  2 

2.3 

2.4 

2.5 

2.6 

2.7 

2.8 

2.9 

3.0 

l 

0.  04042 

0.  03594 

0.03197 

0.  02846 

0.02535 

0.02259 

0.  02014 

0.01796 

0.  01603 

0.01431 

2 

0.  12246 

0.11080 

0.  10026 

0.  09072 

0.  08209 

0.07427 

0.  06721 

0.  06081 

0.  05502 

0.  04979 

3 

0.  24066 

0.  22139 

0.20354 

0.  18704 

0.17180 

0.15772 

0.14474 

0.  13278 

0.12176 

0.11161 

4 

0.  37962 

0.  35457 

0.  33085 

0.30844 

0.  28730 

0.26739 

0.  24866 

0.23108 

0.21459 

0.19915 

5 

0.  52099 

0.49337 

0.46662 

0.44077 

0.41588 

0.  39196 

0.  36904 

0.  34711 

0.  32617 

0.  30622 

6 

0,  64963 

0.  62271 

0.  59604 

0.  56971 

0. 54381 

0.51843 

0.49363 

0.46945 

0.44596 

0.42319 

7 

0.  75647 

0.  73272 

0.  70864 

0.68435 

0.  65996 

0.  63557 

0.61127 

0.58715 

0.  56329 

0.53975 

8 

0.83864 

0.81935 

0.  79935 

0.77872 

0. 75758 

0.  73600 

0.  71409 

0.69194 

0.66962 

0.64723 

9 

0.  89776 

0.88317 

0.  86769 

0.85138 

0.83431 

0.81654 

0.  79814 

0.  77919 

0.  75976 

0.  73992 

10 

0.  93787 

0.92750 

0.91625 

0.90413 

0.  89118 

0.  87742 

0.  86291 

0.84768 

0.83178 

0.81526 

11 

0.96370 

0.95672 

0.94898 

0.  94046 

0. 93117 

0.  92109 

0.91026 

0.  89868 

0.  88637 

0.87337 

12 

0.  97955 

0.97509 

0.  97002 

0.96433 

0.95798 

0.95096 

0.  94327 

0.93489 

0.92583 

0.91608 

13 

0.98887 

0.98614 

0.98298 

0.97934 

0.97519 

0.97052 

0.96530 

0.95951 

0. 95313 

0. 94615 

14 

0.  99414 

0.  99254 

0.99064 

0.  98841 

0.  98581 

0.98283 

0.97943 

0.97559 

0.97128 

0.96649 

15 

0. 99701 

0.99610 

0.99501 

0.99369 

0.99213 

0.  99029 

0.98816 

0.98571 

0. 98291 

0.  97975 

16 

0.  99851 

0.  99802 

0.  99741 

0.  99666 

0.  99575 

0.  99467 

0.99338 

0.99187 

0. 99012 

0.98810 

17 

0.  99928 

0.  99902 

0.  99869 

0.99828 

0.99777 

0.  99715 

0.  99639 

0.  99550 

0.99443 

0. 99319 

18 

0.99966 

0.99953 

0.  99936 

0.  99914 

0.  99886 

0.  99851 

0.  99809 

0.99757 

0.  99694 

0. 99620 

19 

0.  99985 

0.  999  78 

0.  99969 

0.  99958 

0.  99943 

0.  99924 

0.  99901 

0.  99872 

0.99836 

0.  99793 

20 

0.99993 

0.99990 

0.  99986 

0.99980 

0.  99972 

0.99962 

0.99950 

0.99934 

0.  99914 

0.  99890 

21 

0.  9999  7 

0.99995 

0.  99993 

0.  99991 

0.  99987 

0.99982 

0.99975 

0.  99967 

0.99956 

0.99943 

22 

0. 99999 

0.99998 

0.99997 

0.99996 

0.  99994 

0.  99991 

0.  99988 

0.  99984 

0.  99978 

0.  99971 

23 

0.  99999 

0.  99999 

0.99999 

0.  99998 

0.  99997 

0.99996 

0. 99994 

0.  99992 

0.  99989 

0.99986 

24 

0.99999 

0.99999 

0.99999 

0.99998 

0.99997 

0.  99996 

0.  99995 

0. 99993 

25 

0.  99999 

0.  99999 

0.  99999 

0.  9999g 

0. 99998 

0.  99997 

26 

4 

!('■’  '«) 

If  . 

'n  0 

0.  99999 

0.  99999 

0.  99998 

27 

Interpolation  on 

0.  99999 

0.99999 

VI  »-'*■)  v( 

4)[J 

,'o-2)[* 

,>>[> 

X] 

Double  Entry  Interpolation 

V(xV)  V 

<V«  4  ) 

>M  *1.  2)[< 

+  'a- 

«•  +  H VJ 

i[i- « • 

!  *  ’  *■  "  «]  * 
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TABLE  III  (Cant'd) 


9S0 


PltllHAllIl.li  V  1  l  :•  rm.NS 


I..I.I. 

no.:  riMiuutiim 

IN  i  14. I<  31.  ol  1  i;il:i  1  ION.  1\»  o'IH.I  IT.  LAMMA  HVIImN 

u  mi  i.vim-:  vi >  in  1111  hiismin  im-ikiiu  iton 

X2  6.2 

6.4 

6.6 

6.8 

7.0 

7.2 

7.4 

7.6 

7.8 

8.0 

V 

m  3.1 

3.2 

3.3 

34 

3  5 

3.6 

3.7 

3.8 

3.9 

4.0 

1 

0.  01278 

0.01141 

0.01020 

0.00912 

0.  00815 

0.  00729 

0.  00652 

0.00584 

0.00522 

0.  00468 

2 

0. 04505 

0.  04076 

0.  03688 

0.  03337 

0.  03020 

0.  02732 

0.  024  72 

0.  02237 

0.  02024 

0. 01832 

3 

0.10228 

0.  09369 

0.  08580 

0. 07855 

0.  07190 

0.  06579 

0.06018 

0.05504 

0.  05033 

0.  04601 

4 

0.  18470 

0.17120 

0. 15860 

0. 14684 

0. 13589 

0.  12569 

0.11620 

0. 10738 

0.  09919 

0,09158 

5 

0.  28724 

0.26922 

0.  25213 

0.23595 

0.22064 

0.  20619 

0.19255 

0.  17970 

0.16761 

0. 15624 

6 

0.40116 

0.  37990 

0. 35943 

0.  33974 

0.  32085 

0.  30275 

0.28543 

0.  26890 

0.25313 

0.  23B10 

7 

0.  51660 

0.  49390 

0.47168 

0.  45000 

0.42888 

0.  40836 

0.  38845 

0.  36918 

0. 35056 

0. 33259 

8 

0.62484 

0.60252 

0.  58034 

0.  55836 

0.53663 

0.51522 

0.49415 

0.47349 

0.45325 

0.43347 

8 

0.71975 

0.69931 

0.67869 

0.65793 

0.63712 

0.61631 

0.59555 

0.  57490 

0.  55442 

0.  53415 

10 

0.  79819 

0.  78061 

0.  76259 

0.  74418 

0.  72544 

0.  70644 

0. 68722 

0.66784 

0.64837 

0.62884 

11 

0.  85969 

0.  84539 

0.83049 

0.81504 

0.  79908 

0.  78266 

0.  76583 

0.  74862 

0.  73110 

0.  71330 

12 

0.  90567 

0.  89459 

0.  68288 

0.  87054 

0.85761 

0.  84412 

0.  83009 

0.81656 

0.  80056 

0.  78513 

13 

0.  93857 

0.93038 

0.  92157 

0.91216 

0.90215 

0.89155 

0.  88038 

0.  86865 

0. 85638 

0.84360 

14 

0.96120 

0.95538 

0.  94903 

0.  94215 

0.  93471 

0.  92673 

0.91819 

0.  90911 

0.  89948 

0.  88933 

IS 

0.  97619 

0.97222 

0.96782 

0.  96296 

0.  95765 

0.  95186 

0.94559 

0.93882 

0.  93155 

0.  92378 

16 

0.  98579 

0.  98317 

0.98022 

0.  97693 

0.97326 

0.96921 

0.  96476 

0.  959B9 

0.95460 

0.94887 

17 

0.  99174 

0.  99007 

0. 98816 

0.  98599 

0.  98355 

0.98081 

0.97775 

0.97437 

0. 97064 

0.  96655 

18 

0.99532 

0.  99429 

0.99309 

0.99171 

0.99013 

0.  98833 

0.  98630 

0.98402 

0.98147 

0.97864 

19 

0.  99741 

0.99679 

0.  99606 

0.  99521 

0. 99421 

0.  99307 

0.  99176 

0.  99026 

0.  98857 

0.  98667 

20 

0.99860 

0.99824 

0.  99781 

0.  99729 

0.  99669 

0.99598 

0.  99515 

0.  99420 

0.  99311 

0.  99187 

21 

0.  99926 

0.  99905 

0.99880 

0.  99850 

0.99814 

0.  99771 

0.  99721 

0.  99662 

0.  99594 

0.99514 

22 

0.99962 

0.  99950 

0.  99936 

0.  99919 

0.99898 

0.  99873 

0.  99843 

0.  99807 

0.99765 

0.  99716 

23 

0.  99981 

0.  99974 

0.  99967 

0.99957 

0.  99945 

0.  99931 

0.99913 

0. 99892 

0.  99867 

0.  99837 

24 

0.  99990 

0.  99987 

0.99983 

0.99978 

0.99971 

0.  99963 

D.  99953 

0.99941 

0.  99926 

0.  99908 

25 

0.  99995 

0.  99994 

0.  99991 

0.  99989 

0.99985 

0.  99981 

0.99975 

0.  99968 

0.  99960 

0. 99949 

26 

0.  99998 

0.99997 

0.  99996 

0.  99994 

0.  99992 

0.  99990 

0.99987 

0.99983 

0.  99978 

0.  99973 

27 

0. 99999 

0.  99999 

0.  99998 

0.  99997 

0.  99996 

0.  99995 

0.  99993 

0.  99991 

0.  99989 

0. 99985 

28 

29 

30 

0.  99999 

0.  99999 

0.  99999 
0. 99999 

0.  99998 
0.  99999 

0.99998 
0. 99999 
0.  99999 

0.  99997 
0.  99998 
0.  99999 

0.  99996 
0. 99998 
0.99999 

0.  99994 
0.  99997 
0.  99999 

0. 99992 
0.  99996 
0.  99998 

x'2  8.2 

8.4 

8.6 

8.8 

9.0 

9.2 

9.4 

9.6 

9.8 

10.0 

V 

m  4.1 

4.2 

4.3 

4.4 

4.5 

4.6 

4.7 

4.8 

4.9 

5.0 

1 

0.  00419 

0.00375 

0.  00336 

0.  00301 

0.  00270 

0.  00242 

0.  00217 

0.00195 

0.  00175 

0.00157 

2 

0.  01657 

0.  01500 

0.01357 

0.  01228 

0.01111 

0.01005 

0.  00910 

0.  0082) 

0.  00745 

0.  00674 

3 

0.  04205 

0.  03843 

0.  0)511 

0.  03207 

0.  02929 

0.  02675 

0.02442 

0.  02229 

0.  02034 

0.  01857 

4 

0.  08452 

0.  07798 

0.07191 

0.  066)0 

0.  06110 

0.  05629 

0.05184 

0.  04773 

0.  04394 

0.  0404) 

5 

0.  14555 

0. 13553 

0.12612 

0.  11731 

0.  10906 

0.  101)5 

0.  0941) 

0.  08740 

0.  08110 

0.  07524 

6 

0.  22381 

0.  21024 

0.  197)6 

0.  18514 

0. 17)58 

0. 16264 

0.  15230 

0.  14254 

0.  13333 

0.  V465 

7 

0.31529 

0.29865 

0.  28266 

0.  26734 

0.  25266 

0.  23861 

0.  22520 

0.21240 

0.20019 

0.  18857 

8 

0.41418 

0.  39540 

0.  37715 

0.  35945 

0.  34230 

0.  32571 

0.  30968 

0.  2942) 

0.  279)5 

0.  2650) 

9 

0.51412 

0.  49439 

0.  47499 

0.  45594 

0.  4)727 

0.  41902 

0.40120 

0.  38383 

0.  36692 

0.  35049 

10 

0.  60931 

0. 58983 

0.  57044 

0.  55118 

0.  53210 

0.  51)2) 

0.  49461 

0.  47626 

0.  45821 

0.  44049 

11 

0.  69528 

0.  67709 

0. 65876 

0.  640)5 

0. 62189 

0.  60344 

0.  58502 

0. 56669 

0.  54846 

0.  5)0)9 

12 

0.  76931 

0.  75314 

0.  73666 

0.  71991 

0.  70293 

0.  68576 

0.  66844 

0.65101 

0.63350 

0.  61596 

13 

0.  83033 

0.  81660 

0.80244 

0. 78788 

0.  77294 

0.  75768 

0. 74211 

0.  72627 

0.  71020 

0.  6939) 

14 

0.87865 

0.  86746 

0.  85579 

0.84)65 

0.  83105 

0.  81803 

0.  80461 

0.  79081 

0.  77666 

0.  76218 

15 

0.91551 

0.90675 

0. 89749 

0.  88774 

0.  87752 

0.  8668) 

0.  85569 

0.  84412 

0.  83213 

0.  81974 

16 

0.  94269 

0.93606 

0.92897 

0.92142 

0.91341 

0.  90495 

0.  8960) 

0.  88667 

0.  87686 

0.  86663 

17 

0.  96208 

0.  95  72) 

0.  95198 

0.  946)3 

0.  94026 

0.  93)78 

0. 92687 

0.91954 

0.  91179 

0.  90361 

18 

0.97551 

0.9  7207 

0.96830 

0.  96420 

0.  95974 

0.95493 

0.94974 

0. 94418 

0.93824 

0.  93191 

19 

0.  98454 

0.98217 

0.97955 

0.  97666 

0.  97348 

0.  97001 

0.  9662) 

0. 96213 

0.95771 

0. 95295 

20 

0.  99046 

0.  98887 

0.  98709 

0.98511 

0.98291 

0.98047 

0.97779 

0.97486 

0.97166 

0.96817 

21 

0.  99424 

0.99320 

0.99203 

0.  99070 

0.  98921 

0.  9B7S5 

0.98570 

0. 98365 

0.98139 

0.97891 

22 

0.  99659 

0.  9959) 

0.99518 

0.99431 

0.  99333 

0.  99222 

0.  99098 

0.  98958 

0.98803 

0.986)0 

23 

0.  99802 

0.  99761 

0.  99714 

0.99659 

0.99596 

0.  99524 

0.  99442 

0. 99349 

0.  99245 

0.  99128 

24 

0.99888 

0. 99863 

0.  99833 

0.  99799 

0.  99760 

0.  99714 

0. 99661 

0.99601 

0.99532 

0.  99455 

25 

0.99937 

0.  99922 

0.  99905 

0.  99884 

0.  99860 

0.  998)1 

0.  99798 

0.99760 

0.99716 

0. 99665 

26 

0.99966 

0.  99957 

0. 99947 

0.  999)4 

0.  99919 

0.  99902 

0.  99882 

0. 99858 

0.  99830 

0.  99798 

27 

0.  99981 

0.99977 

0. 99971 

0  99963 

0.  99955 

0.  99944 

0.  99932 

0.99917 

0.  99900 

0.  99880 

28 

0.  99990 

0.  99987 

0.  99984 

0.99980 

0.  99975 

0.  99969 

0.  99962 

0.  9995) 

0.  99942 

0.  99930 

29 

0.  99995 

0.  99993 

0.  99991 

0.  99989 

0.  99986 

0.  99983 

0.  99979 

0. 99973 

0.  99967 

0.  99960 

30 

0.  99997 

0. 99997 

0.  99996 

0.  99994 

0.  9999) 

0.  99991 

0.  99988 

0.  99985 

0.  99982 

0.  99977 

V(*M  1  P( 
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TABLE  III  (Cont'd) 


PROBABILITY  FUNCTIONS 


Table  26.7 

t’KOB  ABILITY  IM  KI.lt  \L  OK  x-'-IHSTUIBl  LION.  INCOMIM.KTE  GAMM  A  FI  NOTION 
Cl  Ml  I.ATIX  F.  SI  MS  OF  THE  POISSON  DISTIIIIM  TION 


xM0.5 

11.0 

11.5 

12.0 

12.5 

13.0 

13.5 

14.0 

14.5 

15.0 

V 

i»i  -  5.25 

5.5 

5.75 

6.0 

6.25 

6.5 

6.75 

7.0 

7.25 

7.5 

l 

0.  00119 

0. 00091 

0.  000  70 

0.  0005) 

0.  00041 

0.  000)1 

0.  00024 

0. 0001a 

0.00014 

0. 00011 

2 

0. 00525 

0. 00409 

0.  00318 

3.  30248 

0.  00193 

0. 00150 

0.  00117 

0. 00091 

0. 00071 

0.00055 

3 

0.  01476 

0.01173 

0.00931 

0.  00738 

0. 00585 

0.  00464 

0.  00367 

0. 00291 

0.00230 

0.00182 

4 

0. 03280 

0.02656 

0.  02148 

0.01735 

0. 01400 

0.  01128 

0.  00907 

0.  00730 

0.  00586 

0.  00470 

0. 06225 

0.05138 

0.  042  32 

0.  03479 

0.  02854 

0.02338 

0.01912 

0.01561 

0. 01273 

0.  01036 

6 

0.10511 

0.  08838 

0.  07410 

0.  06197 

0.05170 

0.  04334 

0.  03575 

0.  02964 

0.02452 

0.  02026 

7 

0. 16196 

0.  13862 

0.  11825 

0.  10056 

0.  08527 

0.  07211 

0. 06082 

0.  05118 

0.  04297 

0.  03600 

B 

0.23167 

0.201 70 

0.  1  749*. 

0.  15120 

0. 13025 

0.  11185 

0.09577 

0.08177 

0.06963 

0.  05915 

3 

0.  >1 1 54 

0.  27571 

0.  24299 

0.  21331 

0.  1865  7 

0.  16261 

0.  14126 

0.  12233 

0.  10562 

0.  09094 

10 

0. 39777 

0.  35752 
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0.98365 

0.98139 

0.  97891 

22 

0.99659 

0.  99593 

0.99518 

0. 99431 

0.99333 

0.  99222 

0. 99098 

0.  98958 

0. 98803 

0.  98630 

23 

0.  99802 

0.  99761 

0.  99714 

0.99659 

0.  99596 

0.  99524 

0. 99442 

0.  99349 

0.99245 

0.  99128 

24 

0.99888 

0.  99863 

0.  99833 

0.  99  799 

0.  99760 

0.  99714 

0.  99661 

0.  99601 

0.  99532 

0.  99455 

25 

0.  99937 

0.  99922 

0.  99905 

0.  99884 

0.  99860 

0.  99831 

0. 99798 

0.99760 

0.99716 

0. 99665 

26 

0.  99966 

0.  99957 

0.99947 

0.99934 

0.99919 

0.  99902 

0. 99882 

0.99858 

0.  99830 

0.  99798 

2  7 

0.  99981 

0.  99977 

0.  99971 

0.99963 

0.99955 

0.  99944 

0.  99932 

0.  99917 

0.  99900 

0.  99880 

28 

0,  99990 

0.  99987 

0.  99984 

0.  99980 

0.  99975 

0.  99969 

0. 99962 

0.99953 

0. 99942 

0.  99930 

29 

0.  99995 

0.99993 

0.  99991 

0.99989 

0.  99986 

0.  99983 

0. 99979 

0.99973 

0.  99967 

Q.  99960 

30 

0.  99997 

0.  99997 

0.99996 

0.  99994 

0.  99993 

0.  99991 

0.  99988 

0.  99985 

0. 99982 

0.  99977 

i{i  1,11  [r  (2)  /i*r  ,3  1//.  ">nu/>!(reven,  '•*{►, 
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Table  26.7 


PI«lK\Bll.m  IMKIi  U.  OK  x-'-DISTItlBl  MON.  INCOMPLETE  GAMMA  FI  NCTl'oN 
Cl  Ml  I.U  INE  St  MS  OK  TIIE  POISSON  IMSTRIBI  TION 


x~ 

- 10.5 

11.0 

11.5 

1-2.0 

12.5 

13.0 

13.5 

14.0 

14.5 

15.0 

V 

If* 

5,25 

5.5 

5.75 

0.0 

6.25 

6.5 

6.75 

7.0 

7.25 

7.5 

l 

0.00119 

0.  00091 

0.  00072 

J.  0005) 

0.  00041 

0.  00031 

0.  00024 

0.  00018 

0.  00014 

0.  00011 

2 

0.00S2S 

0.  00409 

0.  00318 

0.  30248 

0.  0014) 

0.  00150 

0.  00117 

0.  00091 

0.  00071 

0.00055 

3 

0.  01*76 

0.0117) 

0.009)1 

0.  00730 

0. 00585 

0.  00464 

0.  0036  7 

0.  00292 

0.002)0 

0.  00182 

4 

0.  03280 

0.  02656 

0.02140 

0.017)5 

0. 01400 

0.  01128 

0.  00907 

0. 00730 

0. 00586 

0.00470 

t. 

0. 0622S 

0.051)8 

0.  0*2  '2 

0.  0)479 

0. 02B54 

0.  02)38 

0.01912 

0.01561 

0.0127) 

0.  01036 

b 

0. 10S11 

0.  0B8'8 

0.  07410 

0.  06197 

0.05170 

0.  0*304 

0. 03575 

0.  02964 

0.02452 

0. 02026 

7 

0. 16196 

0.  13862 

0.  11825 

0.  10056 

0.  08527 

0. 07211 

0.  06082 

0.  05118 

0. 04297 

0.  03600 

8 

0.23167 

0.20170 

0.  1749* 

0.  15120 

0. 1)025 

0.  11185 

0.  09577 

0.  08177 

0.06963 

0.05915 

4 

0.  311S4 

0.  27571 

0.  24299 

0.21)31 

0. 18657 

0.  16261 

0. 14126 

0. 122)3 

0. 10562 

0.  09094 

10 

0.  39777 

0.  35752 

0.  31941 

0.  28500 

0.  25299 

0.  22367 

0.  19704 

0.  17299 

0. 151)8 

0.1)206 

11 

0.  4860S 

0.  *4>2o 

0.  402)7 

0.  )63b4 

0.  32726 

0.  29)3) 

0.26190 

0.2)299 

0.20655 

0. 18250 

12 

0.57218 

0.  52042 

0.  48662 

0.  44568 

0.  40640 

0.  36904 

0.  33)77 

0.  30071 

0.26992 

0.24144 

13 

0.  65263 

0.61082 

0.  56901 

0.  52  704 

0. *071) 

0.  44781 

0.  40997 

0.  37384 

0. 3)960 

0.  30735 

14 

0.  72479 

0.  68604 

0.64639 

0.  606)0 

0. 56622 

0.  52652 

0.48759 

0.  44971 

0.41316 

0.  37815 

IS 

0.  78717 

0.  75254 

0.  71041 

0.  6790) 

0.  6408e 

0.  60230 

0.  56374 

0.  52553 

0.  48800 

0.45142 

16 

0.  83925 

0.  80944 

0.  777b 2 

3. 74)98 

0.  70890 

0. 67276 

0. 63591 

0.  59871 

0.  56152 

0.  52464 

17 

0.  88135 

0.  8Stco 

0.  82942 

0.  80014 

0.  768*3 

0.  7)619 

0.  70212 

0.  66710 

0.63145 

0.  59*48 

18 

0.91436 

0.  894 ‘*0 

0.  8  719* 

0. 84724 

0.  82C3 

0.  79157 

0. 76106 

0.  72909 

0.  69596 

0.  66197 

19 

0.  93952 

0.  92)04 

0.  7050  7 

0.  88*62 

0.  8o)l o 

0. 8)857 

0.81202 

0. 78)69 

0. 75)80 

0.  72260 

20 

0.95817 

0.94622 

C.  9)221 

0. 91008 

0.  89779 

0.  877)8 

0. 85492 

0.  83050 

0.8042  7 

0.  77641 

21 

0.47166 

0.  96274 

0.  9*214 

0.  9)962 

0.  9251  ) 

0.  90062 

0.  89010 

0.  86960 

a  84718 

0.  82295 

22 

0.  38118 

0.  974  75 

0.  90686 

0.  457)8 

0. 94olfl 

0.93)16 

0.91827 

0. 90148 

0.  88279 

0.  86224 

23 

0.  98773 

O.90n> 

0. 47748 

0.  970*7 

0.  90201 

0. 95199 

0.  94030 

0.92687 

0.91165 

0.8946) 

24 

0.99216 

0.40401 

0.  98446 

0.97991 

0.97367 

0.  96612 

0.  95715 

0.  94665 

a  93454 

0.  92076 

2S 

0.99507 

0. 9924* 

0.  99015 

0.  98057 

0.982Q6 

0.  97650 

0.  96976 

0.9617) 

0.95230 

0.941)8 

2b 

0.  99646 

0. 99*55 

0.  99'66 

0.99117 

0.98798 

0.  98397 

0.  97902 

0.97)00 

0.  96581 

0.  957)3 

27 

0.  99815 

0.  99724 

0.  49548 

0. 99429 

0.99208 

0. 98925 

0.98567 

0.98125 

0.  97588 

0.9694) 

28 

0.  9909Q 

0.  998)1 

0. 99744 

0.996)7 

0.  99487 

0. 99290 

0.990)7 

0.98719 

0.  98)24 

0.  97844 

29 

0.  999  35 

0.  99894 

0.  9984o 

0.  9977) 

0.99672 

0.995)8 

0. 9936) 

0.99138 

0.98854 

0. 98502 

30 

0.  4946  3 

0.99940 

0.  44407 

0. 998b0 

Q. 99794 

0.  99704 

0.99585 

0.99428 

0.9922  7 

0.  98974 

X' 

« 15.5 

16.0 

1().5 

17.0 

27.5 

28.0 

18.5 

19.0 

19.5 

20.0 

V 

m 

7.75 

8.0 

8.25 

8.5 

8.75 

9.0 

9.25 

9.5 

9.75 

10.0 

1 

0.00008 

0.  0000o 

0.  00205 

0.  00004 

0. 00003 

0.00002 

0.00002 

0. 00001 

0.  03001 

0.  00001 

2 

0.00043 

0.  000)4 

0,  0302b 

0.  00020 

0. 00016 

0.  00012 

0. 00010 

0.  00008 

0.  00006 

0.  00005 

3 

0.  00144 

o.  oom 

0.  00090 

0.  00071 

0.  00056 

0.  00044 

0.  00035 

0.  00027 

0.  00022 

0.  00017 

4 

0.  003  7  7 

0.  00)02 

0.  00242 

0.00193 

0.0015* 

0.  0012) 

0. 00099 

0.  00079 

0.  00063 

0. 00050 

S 

0.  00843 

0.  00o84 

0. 005*c 

0. 00450 

0. 0036* 

0.  00295 

0. 002)8 

0.00192 

0.  00155 

0. 00125 

b 

0.01670 

0.01)7* 

0.  01131 

0.  00928 

0. 00761 

0.  00b2) 

0. 00510 

0.  00415 

0.  00)40 

0. 00277 

7 

0.  03010 

0.  02512 

0.  0209 2 

0.017*0 

0.01*44 

0.01197 

0.  00991 

0.00819 

0. 00676 

0.  00557 

0 

0.  05012 

0. 04238 

0.  0)*7o 

0.  03011 

0.025)0 

0.  02123 

0.  01777 

0.  0 1486 

0.  01240 

0.  01034 

4 

0.  07809 

0.  06688 

0.  0571* 

0.  0*872 

Q. 04144 

0.0)517 

0. 02980 

0.02519 

0.  02126 

0.01 791 

10 

0.  11487 

0.  0996) 

0.  08614 

0.  074)6 

0.  06401 

0.0549b 

0. 04709 

0.  04026 

0.  03435 

0.02925 

11 

0.  16073 

0.  14113 

0.  12356 

0.  10780 

0. 0939) 

0.  08158 

0.  07068 

0.06109 

0.0^269 

0. 04S34 

12 

0.  2152? 

0.  19124 

0.  10919 

0. 14960 

0.1)174 

0.  11569 

0.  10133 

0.  08853 

0.  07716 

0.  06709 

13 

0.  27719 

0.  2491  i 

0.  22318 

0. 199)0 

0. 17744 

0.  15752 

0. 1394* 

0.12)10 

0. 10840 

0.  09521 

14 

0.  34485 

0.  313)7 

0.  28)80 

0.25618 

0. 23051 

0.  20670 

0.  18495 

0.  16495 

0.  14671 

0. 13014 

IS 

0.  41604 

0.  38205 

0.  34**2 

0.  31886 

0.  28906 

0.  26267 

0.23729 

0.21373 

0. 19196 

0.17193 

16 

0.  4883  7 

0. 4*246 

0. 41864 

0. 38560 

0. 35)98 

0.  32)90 

0. 2954* 

0.  26866 

0.  24)59 

0.  22022 

1  7 

0.  55951 

0.  52303 

0.  480  71 

0.45*37 

0.  42102 

0.  38884 

0.  35797 

0.  3285) 

o: 30060 

0.27423 

18 

0.  62740 

0.  54255 

0.  5*770 

0.52)11 

0.  48902 

0,  45565 

0.  42320 

0.  39182 

0.  36166 

0. 33282 

t  4 

0.  69033 

0.  65728 

0. 62370 

0.  58987 

0.5560) 

0.  522*4 

0.489)1 

0.  45684 

0.  42521 

0.  39458 

20 

0.  74712 

0.  7166? 

0.  60*16 

0.65297 

0.62031 

0.  58741 

0. 55451 

0.52183 

0.  48957 

0.  45  793 

21* 

0.  79705 

0. 76965 

0.  7409) 

0. 71111 

0.  600)9 

0.  6*900 

0.  61718 

0.  58514 

0.  55310 

0.52126 

23 

0.  0)990 

0. 81589 

0.  790)2 

0.  76)36 

0. 7)519 

0.  70599 

0.67597 

0.645)3 

0.61428 

0.58)04 

2* 

0.  87502 

0.8552  7 

0.83)04 

0.  8092S 

0. 78402 

0.  757*9 

0. 7298) 

0. 70122 

0.67185 

0.64191 

24 

0.  90527 

0.  80808 

0.  80419 

0.  8*866 

0.82657 

0.80301 

0. 77810 

0. 75199 

0.  72483 

0. 69678 

2S 

0. 92891 

0.91483 

0.  84912 

0.  88179 

0.  86287 

0.8*239 

0.  82044 

0.  79712 

0.  77254 

0. 74683 

2h 

0.  9*749 

0.  9)620 

0.  92)41 

0.  90908 

0.  89)20 

0.  87577 

0.  85683 

0.  8)64) 

0.  81464 

0.  79156 

27 

0.96182 

0.  95295 

0.  942  74 

0.9)112 

0.91806 

0.  90352 

0. 86750 

0.87000 

0.85107 

0.83076 

28 

0.  97266 

0.96582 

0.  95782 

0.9*859 

0.9)805 

0.92615 

0.91285 

0.69814 

0.88200 

0.  86446 

29 

0.98071 

0.9  7554 

0.  904)9 

0.  96218 

0.95)83 

0.  94427 

0.  93344 

0. 92129 

0.  90779 

0.  8929) 

30 

0.  98659 

0.9827* 

0.97810 

0.  97258 

0.  96608 

0.9585) 

0. 94986 

0.94001 

0.92891 

0.91654 

* 

«’  *  v0 

o 

Interpolation  on  <’ 


Q 'A-  ►)  y(s>„  4)[>2] •<)(*?,.„  2)[«  >„)[l  *>2] 

Double  Entry  Interpolation 


n,x- 


V(v.  •»>[>-]  vhi 


*„  2)[«  »2  it*]y(v;  i,  i)[X 

*  i.«].y(»2  ii •  i)[.J ,4'2 ■  J 


ie*J 


jt 


i 
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U  Ml  1 

VMM.  M  MS  OK  1  TIE  l*OISNO\  111" 

1  Itlltt  TION 

«-  21 

•>•> 

23 

24 

25 

26 

27 

28 

29 

30 

. 

m  10.5 

11.0 

11.5 

12.0 

12.5 

13.0 

13.5 

14.0 

14.5 

15.0 

1 

0.  00001 

2 

0.  00003 

0.  00002 

0.  00001 

0.00001 

3 

0. 00011 

0.  00007 

0.  00004 

0.  00003 

0.00002 

0.  00001 

0.  00001 

4 

0.  00032 

0.  00020 

0.  00013 

0. 00008 

0.00005 

0.  00003 

0.00002 

0. 00001 

0.00001 

0. 00001 

5 

0.00081 

0.  00052 

0.  00034 

0. 00022 

0.  00014 

0.  00009 

0.00006 

0.00004 

0.  00002 

0. 00002 

6 

0.  00184 

0.00121 

0.00080 

0. 00052 

0.00034 

0.  00022 

0.00015 

0. 00009 

0.  00006 

0.00004 

7 

0. 00377 

0. 00254 

0.  00171 

0. 00114 

0.  00076 

0.  00050 

0.00033 

0. 00022 

0.  00015 

0.  00010 

8 

0.  00715 

0.  00492 

0.  00336 

0.  00229 

0.  00155 

0.  00105 

0.  00071 

0.  00047 

0.  00032 

0.  00021 

9 

0.  01265 

0.  00888 

0.00620 

0. 00430 

0.  00297 

0.  00204 

0.  00140 

0.  00095 

0.  00065 

0.  00044 

10 

0.  02109 

0.  01511 

0.  01075 

0. 00760 

0.  00535 

0.  00374 

0.  00260 

0.  00181 

0.00125 

0.  00086 

11 

0.  03337 

0.  02437 

0.  01768 

0.01273 

0.00912 

0.  00649 

0.  C0460 

0.  00324 

0.00227 

0.00159 

12 

0.  05038 

0.  03752 

0.  02773 

0.  02034 

0.  01482 

0.  01073 

0.  00773 

0.  00553 

0.  00394 

0.00279 

13 

0.  07293 

0.  05536 

0.  04168 

0. 03113 

0.  02308 

0.  01700 

0.  01244 

0. 00905 

0.00655 

0.  00471 

14 

0.  10163 

0.  07861 

0.  06027 

0. 04582 

0.  03457 

0.  02589 

0.  01925 

0. 01423 

0.01045 

0.00763 

15 

0.  13683 

0.  10780 

0.  08414 

0.  06509 

0.  04994 

0.  03802 

0.  02874 

0.  02157 

0.  01609 

0.  01192 

16 

0.  1  7851 

0.  14319 

0.  11374 

0.  08950 

0.  06982 

0.  05403 

0.  04148 

0. 03162 

0.  02394 

0.01800 

17 

0.  22629 

0.  18472 

0.  14925 

0. 11944 

0.  09471 

0.  07446 

0.  05807 

0.04494 

0.  03453 

0.02635 

18 

0.  27941 

0.23199 

0.  19059 

0.15503 

0.  12492 

0.  09976 

0.  07900 

0.  06206 

0.  04838 

0.  03745 

19 

0.  33680 

0.  28426 

0.  23734 

0. 19615 

0.  16054 

0.  13019 

0.  10465 

0.  08343 

0.  06599 

0.05180 

20 

0.  39713 

0.  34051 

0.28880 

0.  24239 

0.20143 

0.  16581 

0.  13526 

0.  10940 

0.  08776 

0.  06985 

21 

0.45894 

0.  39951 

0.  34398 

0.  29306 

0.  24716 

0.  20645 

0.17085 

0. 14015 

0.11400 

0.09199 

22 

0.  52074 

0.  45989 

0.40173 

0.  34723 

0.  29707 

0.25168 

0.21123 

0.  17568 

0.14486 

0.11846 

23 

0.  58109 

0.  52025 

0.  46077 

0.  40381 

0.  35029 

0.  30087 

0.  25597 

0.21578 

0.  18031 

0.14940 

24 

0.  63873 

0.  57927 

0.51980 

0.46160 

0.  40576 

0.  35317 

0.  30445 

0. 26004 

0.22013 

0.  18475 

25 

0.69261 

0.63574 

0.  57756 

0.51937 

0.  46237 

0.  40760 

0.  35588 

0.  30785 

0.  26392 

0.22429 

26 

0.  74196 

0.  68870 

0.63295 

0.  57597 

0.51898 

0.  46311 

0.  40933 

0.  35846 

0.  31108 

0.  26761 

27 

0.  78629 

0.  73738 

0.68501 

0.63032 

0.  57446 

0.51860 

0.  46379 

0.41097 

0.  36090 

0.31415 

28 

0.  82535 

0.78129 

0.  73304 

0.68154 

0.62784 

0.  57305 

0.51825 

0. 46445 

0.  4  1253 

0. 36322 

29 

0.  85915 

0.82019 

0.  77654 

0.  72893 

0.  67825 

0.  62549 

0.57171 

0.51791 

0.46507 

0.41400 

30 

0.  88789 

0.  85404 

0.81526 

0.  77203 

0.  72503 

0.  67513 

0.  62327 

0.57044 

0.  51760 

0.46565 

,7  31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

. 

m  15.5 

16.0 

16.5 

17.0 

17.5 

18.0 

18.5 

19.0 

19.5 

20.0 

5 

0.  00001 

0.00001 

6 

0.  00003 

0.00002 

0.  00001 

0. 00001 

7 

0.  00006 

0.  00004 

0.  00003 

0. 00002 

0.00001 

0.  00001 

8 

0.  00014 

0.  00009 

0.  00006 

0.  00004 

0.  00003 

0.  00002 

0.  00001 

0.00001 

9 

0.  00030 

0.00020 

0.00013 

0.  00009 

0.00006 

0.  00004 

0.  00003 

0.  00002 

0.  00001 

0.00001 

10 

0.  00059 

0.  00040 

0.  00027 

0.  00019 

0.  00012 

0.  00008 

0.  00006 

0.  00004 

0.  00003 

0.  00002 

11 

0.  00110 

0.  00076 

0.00053 

0.00036 

0.  00025 

0.  00017 

0.  00012 

0. 00008 

0.  00005 

0.  00004 

12 

0.  00197 

0.00138 

0.00097 

0. 00068 

0.00047 

0.00032 

0.  00022 

0. 00015 

0.  00011 

0.00007 

13 

0. 00337 

0.  00240 

0.00170 

0.00120 

0.  00085 

0.  00059 

0.  00041 

0.00029 

0.  00020 

0.00014 

14 

0.  00554 

0.00401 

0.  00288 

0. 00206 

0.00147 

0.  00104 

0.  00074 

0.00052 

0.  00036 

0.00026 

15 

0.  00878 

0.00644 

0.  00469 

0.  00341 

0.  00246 

0.00177 

0.00127 

0. 00090 

0.  00064 

0.  00045 

16 

0.  01346 

0.01000 

0.  00739 

0. 00543 

0.  00397 

0.  00289 

0.00210 

0.00151 

0.  00109 

0.  00078 

17 

0.01997 

0.01505 

0.  01127 

0.  00840 

0.00622 

0.00459 

0.  00337 

0. 00246 

0.  00179 

0.00129 

18 

0. 02879 

0.02199 

0.  01669 

0.01260 

0.00945 

0.  00706 

0.  00524 

0.00387 

0.  00285 

0.  00209 

19 

0.  04037 

0.03125 

0.  02404 

0.  01838 

0.01397 

0.01056 

0.  00793 

0.00593 

0.00442 

0.  00327 

20 

0.05519 

0.  04330 

0.03374 

0.  02613 

0.  02010 

0.01538 

0.  01170 

0.00886 

0.  00667 

0.  00500 

21 

0.07366 

0.  05855 

0.  04622 

0.  03624 

0.  02824 

0.  02187 

0.  01683 

0.  01289 

0.  00981 

0.  00744 

22 

0.  09612 

0.  07740 

0.06187 

0.  04912 

0.03875 

0.  03037 

0.  02366 

0.01832 

0.01411 

0.01081 

23 

0.  12279 

0.10014 

0.  08107 

0. 06516 

0.  05202 

0.04125 

0.  03251 

0.  02547 

0.  01984 

0.  01537 

24 

0.  15378 

0.  12699 

0.  10407 

0.  08467 

0.  06840 

0.05489 

0.  04376 

0.  03467 

0.  02731 

0.  02139 

25 

0.18902 

0.15801 

0.13107 

0.  10791 

0.  08820 

0.07160 

0.  05774 

0. 04626 

0.  03684 

0.  02916 

26 

0. 22827 

0.19312 

0. 16210 

0.  13502 

0.11165 

0.09167 

0.07475 

0. 06056 

0.  04875 

0.  03901 

27 

0.27114 

0.23208 

0.  19707 

0.  16605 

0.13887 

0. 11530 

0.  09507 

0.  07786 

0.  06336 

0.05124 

28 

0.  31708 

0.  27451 

0.  23574 

0.  20087 

0. 16987 

0. 14260 

0. 11886 

0.  09840 

0.  08092 

0.  06613 

29 

0.  36542 

0.31987 

0.  27774 

0. 23926 

0.  20454 

0.  17356 

0. 14622 

0.  12234 

0.10166 

0.  08394 

30 

0.41541 

0.  36753 

0.  32254 

0.  28083 

0.  24264 

0.  20808 

0.17714 

0.14975 

0.  12573 

0. 10486 

31 
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TABLE  III  (Cont'd) 


PROBABILITY  FUNCTION'S 
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PROBABILITY  INTEGRAL  OF  x-’-DISTRIIUHON,  INCOMPLETE  GAMMA  FI  ACTION  Table  26.7 
CUMULATIVE  SUMS  OF  TIIE  POISSON  DISTRIBUTION 


i 

i 


*2=42 

44 

46 

48 

50 

52 

54 

56 

58 

60 

V 

m -21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

10 

0.  00001 

11 

0.  00002 

0.00001 

12 

0.  00003 

0.  00002 

0.00001 

13 

0.00006 

0.  00003 

0.  00001 

0.00001 

14 

0.00012 

0.00006 

0.00003 

0.  00001 

0.  00001 

15 

0.  00023 

0.00011 

0.  00005 

0.  00003 

0.  00001 

0.  00001 

16 

0.00040 

0.  00020 

0.  00010 

0. 00005 

0.  00002 

0.  00001 

0.  00001 

17 

0.00067 

0.  00034 

0.00017 

0.  00009 

0.  00004 

0.  00002 

0.  00001 

0.  00001 

18 

0.  00111 

0.  00058 

0.  00030 

0. 00015 

0.  00008 

0.  00004 

0.  00002 

0. 00001 

19 

0.  00177 

0.00094 

0.00050 

0.  00026 

0.  00013 

0.  00007 

0.00003 

0.  00002 

0.  00001 

20 

0.00277 

0.00151 

0.00081 

0.  00043 

0.  00022 

0.  00011 

0.  00006 

0.  00003 

0.  00001 

0.  00001 

21 

0.  00421 

0.  00234 

0.  00128 

0.  00069 

0.  00036 

0.  00019 

0.  00010 

0.  00005 

0. 00003 

0.  00001 

22 

0.  00625 

0.00355 

0.00198 

0.  00109 

0.  00059 

0.  00031 

0.  00016 

0.  00009 

0.  00004 

0.  00002 

23 

0.  00908 

0.  00526 

0.  00299 

0.  00167 

0.  00092 

0.  00050 

0.  00027 

0.  00014 

0.  00007 

0.  00004 

24 

0.  01291 

0.  00763 

0.  00443 

0.  00252 

0.00142 

0.  00078 

0. 00043 

0. 00023 

0.  00012 

0.  00006 

25 

0.  01797 

0.01085 

0.  00642 

0.  00373 

0.  00213 

0.  00120 

0. 00066 

0.  00036 

0.  00020 

0. 00011 

26 

0.02455 

0.  01512 

0.  00912 

0.  00540 

0. 00314 

0.  00180 

0.  00102 

0.  00056 

0.  00031 

0.  00017 

27 

0.  03292 

0.  02068 

0. 01272 

0.  00768 

0.  00455 

0.  00265 

0.  00152 

0.  00086 

0.  00048 

0.  00026 

28 

0.  04336 

0.  02779 

0.  01743 

0.  01072 

0.  00647 

0.  00384 

0.  00224 

0. 00129 

0.  00073 

0.  00041 

29 

0,  05616 

0.  03670 

0.  02346 

0.  01470 

0.  00903 

0.  00545 

0. 00324 

0.  00189 

0. 00109 

0.  00062 

30 

0.07157 

0.  04769 

0.  (>3107 

0.  01983 

0.  01240 

0.  00762 

0.  00460 

0.  00273 

0.  00160 

0.  00092 

x2  -  62 

64 

66 

68 

70 

72 

74 

76 

V 

m  31 

32 

33 

34 

35 

36 

37 

38 

21 

0.  00001 

22 

0.  00001 

0.00001 

23 

0.  00002 

0.  00001 

0.  00001 

24 

0.  00003 

0.00002 

0.  00001 

25 

0.  00006 

0.  00003 

0.  00002 

0.  00001 

26 

0.  00009 

0.00005 

0.  00003 

0.  00001 

0.  00001 

27 

0.00014 

0.00008 

0. 00004 

0. 00002 

0.  00001 

0.  00001 

28 

0.  00023 

0.00012 

0.00007 

0.  00004 

0.  00002 

0.  00001 

0.  00001 

29 

0.00035 

0.00019 

0.  00011 

0. 00006 

0.  00003 

g.  00002 

0. 00001 

30 

0.00052 

0.00029 

0.  00016 

0.  00009 

0.  00005 

0.  00003 

0.  00001 

0.00001 

V(*l  O'l- /’(**!») 


Interpolation  on  *-’ 


i  r  i 

<//  i:  <■ 
y-<> 


my  j\(>  even, 


--P,  »'  D2) 


V(x2i")  V(\2  '<>  4)[‘d2]+V(x^o-2)[«»--*2]+V(.*o;‘o)[l- 


Double  Entry  Interpolation 


V  (x*!*)"*?  4)[y  *']f '■>  (  xn  -2)[* (  xjjj  -o-l)[i  "S  “■  +  "*] 

t-V  ( x2'  “o)[l  -  ‘  j  f  "'*  ]  +  V  ( 1  “d + 1 )  ['.  "'2e.', 


^  ► 
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TI-59  METHOD  (ML-09) 


The  TI-59  Programmable  Calculator  has  a  numerical  integration  program 
it  its  Master  Library,  ML-09.  Included  in  PL-4  are  pages  29-31  of  the 
TI-59  User's  Manual.  The  method  uses  Simpson's  Rule.  and 

Q(a,i)  can  be  computed  using  ML-09.  The  integration  limits  will  be 
different  for  the  two  Incomplete  Gamma  Functions. 

Y(a,T) 

Lower  Limit  XQ  0 

Upper  Limit  3^  t 

*A  suitably  large  real  number  must  be  found  to  represent  °°. 

This  would  be  a  value  such  that  a  large  value  would  result  in 

a  negligible  increase  in  the  value  of  the  integral.  A  value 

of  X  >  227  will  cause  underflow, 
n 

A  sufficiently  large  value  for  n  must  be  established.  The 
larger  the  value  of  n  used,  the  better  the  accuracy  of  the 
integral.  The  computation  time  also  increases. 

A  listing  for  f(x)  is  also  included. 


Q(a,T) 

T 

co"k 
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PL  -4 


ML-09 


SIMPSON'S  APPROXIMATION 
(CONTINUOUS) 


This  program  may  be  used  to  approximate  the  integral,  I,  of  a  function  defined  by  the  user,  over  an 
interval  x0  to  xn,  using  Simpson's  Rule. 

*n 

I  = J  f(x)dx 

*o 

The  function  f(x)  must  be  expressed  as  a  sequence  of  keystrokes  in  the  user  program  memory. 


SoOd  State  Softwares 
SIMPSOtfS^APPROXlMATIOfirtCOMINUOUSHp.  ML-09  >• 


procedure 


Imti 


S-t«  *  It-vn  mode 
Us-  A*  as  k*u.»i 

Ent  •  f ( x )  a-  3  senes  of  keystrokes.  Do  not  use  [=J 
or  CIO  .  Do  no'  us-  reijist-i  >  0  S 

E mi  f ( x )  with  [ ir» v 

L P3vi>  I,: jni  mode 
Select  program 
Ent-r  low-r  limit 
Enl-r  up))-!  limit 

Enter  n(n  2.4.6,  •  •  •.displ'.y  flashes  if  not 
legal  entry) 

Con'(i -re  integral 

Eoi  an'  ,v  inter ,  ’*  or  a  nr ..  n,  rep-.'*  Steps  7  11. 


ENTER 


PRESS 

gSTj' 

tlRNj 

[M  pi 

[?nd  1  LL.j 


[mv;  fs^ej 


[?•«*]  ti.’loA 

[*  I 

[B] 

Cc  ! 


DISPLAY 

0 

000  00 

001  00 
002  00 


xxx  00 
0 


N  Ht  1.  Evaluate  expressions  using  parentheses  only 
2  Running  time  is  dependent  on  input  data. 
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PL-  + 


r  •'...nple: 


Evaluate 


f\ 

**  r\ 


--  —  dx  using  two  subintervals, 
cosx  +  2 


REF. 

ENTER 

PRESS 

DISPLAY  COMMENTS 

[rst] 

0. 

s' 

irn] 

000  00 

[2ndj  m 

001  00 

[  2nd)  na 

002  00 

[2nd] 

003  00 

l  <  ]  [M  Su3 

005  00  Key  in  <(x). 

r+  ]  L  2  j  LO 

008  00 

( VxJ  :  invJ  issg 

010  00 

lrn] 

0. 

!2"dJ  pj^3o9 

0.  Select  piogram 

1 

0 

(  *J 

0.  Xq 

>d]  IZi  [+J 

3.14  1592654 

2 

2 

[  =3  1-  bj 

1  5/0/96327  Xj  |it/2) 

3 

2 

[C] 

.7853981634  h 

4 

1  0 J 

0  004998903  1 

'  V-i.  p.  • 

■  tout  shown  is  obljined  using  the  print 

routine  of  Program  01 . 

Register  Contents 

R  00 

Ros 

n 

R  io 

R  o  i 

*0 

R  06 

R . . 

R  0  2 

R  o  i 

r,2 

R  03 

h 

R  o» 

R  13 

R  04 

1 

R  09 

R  i« 

OPTIONAL 

PRINTOUT* 

REF.  PRINT 

1  0 

0. 

2  1670/96327 
1670796  *27 

3  2 

7S6j:)31i'34 

4  7  as  a.- Hi  m 
or.cjre 
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EXAMPLE 


FIND:  y(1.5,  0.9)  and  Q  (1.5,  0.9) 
a  =  1 . 5 

x  =0  ;  x  =0.9 

o  n 


n 

y(1.5,  0.9) 

16 

0.3401 

32 

0.34087 

6-4 

0.3411193 

128 

0.341207 

512 

0. 341249 

Y(l. 

5,  0.9)  =  0.341249 

Accuracy  could  be  improved  by  increasing  n  still  further  or  by 
breaking  the  integral  into  two  integrals  with  limits  from  0-0.5 
and  0.5  -  0.9  respectively, 
n  =  512 

Y0-0  5  =  °* 1761333892 

Y  =  0.1651191608 

0. 5-0.9 

Y (1 . 5,  0.9)  =  0.34125255 

n  =  128  for  each  of  the  following  integrals: 

Q0  9_4  -J'*  y  dx  *  0.504195085 
0.9 

Q.  =  0.0397718 
4-o 

Q8_20=  0.001005 
Q20_4o  =  0.0000000094 

Q(1.5,  0.9)  =  0.54497189 

y(1.5,  0.9)  +  Q  (1.5,  0.9)  -  0.88622444 
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EXAMPLE 

FIND:  y(5.64,  8) 

xQ  =  0;  xn  =  8;  n  =  128 
y(5.64,  8)  =  55.57279  (55.5728  to  six  places) 

FIND:  Q  (5.64,  8) 

x  =  8;  x  =  “  (say,  x  =  64) 

o  n  11 

n  =  128 

Q  (5.64,  8)  =  10.14056  (10.1406  to  six  places) 
y(5. 64,  8)  +  Q  (5.64,  8)  =  65.7134 

FIND:  F ( 5 . 64) 

xo  =  0  >  xu  =  64 ’  n  =  256 
r(5 . 64)  '  65.7134 
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iGrORE  4  GfAMMA  FukICTIOMS  For  <X=5.64 


B.  TI-59  METHOD  (MU-17) 


The  new  Math  Utility  module  can  also  be  used  for  computing  the  Gamma 
Functions.  The  numerical  integration  program,  MU-17,  is  much  faster  in 
execution  time  than  ML-09  (i.e.  PL-4)  and  simpler  to  use.  This  method 
uses  the  systematic  procedure  called  Romberg  Integration.  The  subroutine 
for  f(x)  must  be  entered  (see  PL-5)  along  with  the  integration  limits 
and  expected  accuracy  of  the  result.  The  program  then  quickly  computes 
the  integral. 
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f 

I 


I 


1 

i 


■ 


PL-5 


MU-17 

LISTING  FOR  f(x) 

RST 

LRN 

ENTER 

f(x) 

LOC 

CODE 

KEY 

ENTER 

PRESS 

DISPLAY 

42 

-  rn 

"  a  " 

R/S 

"a-1" 

z 

22 

UU 

2nd  Pgts  17 

4  3 

=  4  L 

a,  lower  limit 

A 

i  l  •:  i 

75 

- 

b,  upper  limit 

B 

0  0  7 

5  4 

r  ;r  1 

£ 

C 

U  U  r 
f;  1  f: 

0  6 

0  A 

D 

Integral 

s~i  i  t 

i  W  { 

42 

H 

n  ^  s“; 

.0  A 

5  3 

R  L  L 

~  1  l 

45 

. . ; J 1 

NOTE:  b  <  227 

otherwise 

underflow  will 

020 

0  A 

RUL 

0  6 

occur  in 

computing 

f(x) 

022 

4  2 

STO 

024 

72- 

n  >■  5 

.  i .  *i 

0  2 1 
f\  v  7 

m- 

R  L-  L 

0  2  8 

94 

•+•  — 

0i"4 

21  cl 

T  k? •' 

030 

I  h  y 

03 1 

S  cl 

032 

tifi 

03  2 

4  -I* 

RCL 

0  2  4 
035 

02 

54 

0  v 

036 

92 

RTN 

LRN 

RST 


I 


«  »: 


r* 


EXAMPLES : 

FIND:  y 

(5.64,  8) 

c 

y  (5.64,  8) 

APPROXIMATE 

EXECUTION 

TIME 

.01 

55.57467363 

45  sec 

IE-3* 

55.574674  , 

1  min 

IE-4 

55.572748 

1  min  20  sec 

IE-5 

55.5727 9 

2  min  30  sec 

IE-6 

55.57279 

4  min 

FIND:  Q  (5.64,  8) 

b  =  upper 

limit 

G  b 

0  (5.64,  8) 

10.140547 

appropriate 

EXECUTION 

TIME 

IE- 5  25 

2  min  30  sec 

IE-6  50 

10.140599 

5  min 

IE-6  200 

10.140599 

30  min 

*  IE-3  =  10"3 

;  IE-4  =  10'4 

;  etc. 

NOTE:  T(5.64)  =  65.71338  =  55.57279  + 

10.140599 
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PL-& 


BASIC  LANGUAGE  PROGRAM  LISTING  FOR  COMPUTING  THE 
INCOMPLETE  GAMMA  FUNCTIONS  USING 
SIMPSON’S  NUMERICAL  INTEGRATION  TECHNIQUE 


10  REM  A  =  ALPHA 
20  REM  T  =  TAU 

30  REM  I  =  INCOMPLETE  GAMMA  FUNCTION  WITH  ARGUMENT  ALPHA,  TAU 
40  REM  X0  =  LOWER  INTEGRAL  LIMIT 
50  REM  XI  =  UPPER  INTEGRAL  LIMIT 

60  REM  N  =  NUMBER  OF  INTEGRATION  SUB  INTERVALS (MUST  BE  EVEN  NUMBER ) 
70  REM  INTEGRATION  CALCULATIONS  ARE  BASED  ON  SIMPSON’S  RULE 


80  A 

i  =  5.64 

90  T 

'  =  8.0 

100 

X0  =  0 

110 

X  1  =8 

120 

N  =  32 

130 

H= ( X 1  - X0) /N 

140 

S1=0 

150 

S2  =  0 

160 

DIM  Y (500) 

170 

X  =  X0 

180 

Y  (0)  =  <  X  t  ( A  -  1 )  )  •*  ( EXP  (  - 

X)  ) 

190 

N2=N- 1 

200 

Y (N)  =  (XI t  (A-l) ) * ( EXP ( 

-XI ) 

210 

FOR  J=1  TO  N2 

220 

X  =  X  +  H 

230 

Y ( J ) =  <  X  *( A - 1 ) ) * (EXP  < - 

X)  ) 

240 

NEXT  J 

250 

FOR  J  =  2  TO  (N-2)  STEP 

2 

260 

S 1 =  S 1 +  Y (J) 

270 

NEXT  J 

280 

S 1 =  Y ( 0 ) +  Y ( N )  +  ( 2*S l ) 

290 

FOR  J=1  TO  N2  STEP  2 

300 

S2=S2+ Y (J) 

310 

NEXT  J 

j20 

S2=4*S2 

330 

I  =  (H/3) # ( S 1  +  S  2 ) 

340 

PRINT  "ALPHA  =  " i A 

350 

PRINT  "TAU  *  "»T 

360 

PRINT  "X (LOWER)  =  " l 

X0 

370 

PRINT  "X (UPPER)  =  "» 

XI 

380 

PRINT  "N  =  " 5  N 

390 

PRINT 

400 

PRINT 

410 

PRINT  "INTEGRAL  =  "i 

I 

420 

END 

45 


V 


HP-34C  METHOD 


The  program  for  computing  the  Incomplete  Gamma  Functions  on  the  HP-34C 
is  listed  in  PL-15  which  also  computes  other  functions.  The  following 
listing  for  y(cx,t)  and  Q(a,r)  are  included  for  information  only.  This 
computer  has  a  very  fast  and  accurate  numerical  integration  scheme 
built  into  it. 


46 


HP34C 

INCOMPLETE  GAMMA  FUNCTIONS 
Y(a.t)  ,  Q(a ,t) 


Y(a,t)  : 

Put  value  of  a  in  X-register 
Press  ENTER  i 

Put  value  of  t  in  X-regiscer  * 

Press  GTO  4 
Press  R/S 

Display  will  be  y(ci»'r) 


Q(a,r)  : 

Put  value  of  a  in  X-register 
Press  ENTER  t 

Put  value  of  t  in  X-register  * 
Press  GTO  3 
Press  R/S 

Display  will  be  Q(a,r) 


NOTE:  The  upper  integration  limit  that  represents  «  is  stored  in  R  ^  . 

For  example:  upper  limit  =  50.  Put  50  in  the  X-register.  Press 
STO  .1  .  The  upper  limit  should  be  less  than  227;  otherwise  under¬ 
flow  will  occur  in  computing  f(x).  An  upper  limit  of  50  has  proven 
to  give  accurate  results.  See  page  44. 

*  Press  f  SCI  7  (for  accuracy  of  7  significant  figures) 
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EXAMPLES 


FIND:  y(5.64,  8) 


APPROXIMATE 

EXECUTION 

ACCURACY 

y(5.64,  8) 

TIME 

f  SCI  5 

55.5728 

4  min 

f  SCI  7 

55.57279 

7  min 

FIND:  Q(5 . 64 ,  8) 

Upper  limit  (i.e.  °°)  =  50 

APPROXIMATE 

EXECUTION 

ACCURACY  Q(5.64,  8)  TIME 


SCI 

5 

10.1406 

5 

min 

SCI 

6 

10.14060 

6 

min 

SCI 

7 

10.14059 

7 

min 

i 


4 


a 

•* 


HP-67  Method 


Simpson's  Numerical  Integration  Technique  is  most  conveniently  found  as 
program  19-01  of  Hewlett-Packard's  Math  Pac  1  package. 

Included  also  in  PL-7  is  a  program  listing  for  computing  the  Incomplete 
Gamma  Functions  on  the  HP-67. 
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Iser  I nsi ructions 


PL-7 


Program  Description 


••v 


DOS  NOT  USE  THIS-  SPACE* 


•  W-V,'*'  ■  "C ‘  -T  CM*"" 


C.  GAMMA  FUNCTIONS  COMPUTATION  PROGRAM  (BASIC  LANGUAGE) 


Included  in  PL- 8  is  a  computer  program  written  in  Basic  Language 
which  quickly,  conveniently  and  accurately  computes  the  three  Gamma 
Functions:  F(a),  y(a,t),  Q(a,x). 

The  method  for  computing  F(a)  is  as  previously  described.  The 
method  for  computing  the  two  Incomplete  Gamma  Functions  has  not 
been  previously  described.  This  method  uses  an  approximation  formula 
as  described  in  Abramowitz  1  to  calculate  a  "Normalized  Incomplete 
Gamma  Function:  Qo  (a,x).  Qq  (oi,t)  is  extremely  accurate  for 
a  =  integer.  i  may  be  either  an  integer  or  non-integer. 


Q  (a, x)  =  Q(a. i) 

°  (a-1) 1 

(a-1)  !  =F(a)  for  a  =  integer 

The  approximation  formula  is  as  follows: 


Q  (« 


/ 

r<r  -  Z 


GT-  I 


(nr- ') 

1  1  (  IA  -  cx) 


1  £  a  _<  350 

NOTE:  For  most  fatigue  analysis  applications  5  <  a  <  15. 
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C.  GAMMA  FUNCTIONS  COMPUTATION  PROGRAM  (Cont'd) 


i 

i 


k 


For  a  =  non-integer  the  interpolation  method  in  Abramowitz  1 
used . 

Def ine: 


a0  =  integer  part  of  a 
w  =  a  -  a 


I  I  2  1 

CL  (a,x)  =  Q_  (a  -  1,  t)  j  w  -  2  w 


+Q 


,  r  2~j 

(a„  ,t)  '  1  -  w  i 


o  o 


+Qo  (a0  +  1,  t) 


1  2  1  : 

-w  +2W; 


Then 

Q(a , t)  =  I (a)  Qq  (a,x) 

Y(a,t)  =  T (a)  |^1  -  Q  (a>T)j 


Several  computation  runs  are  included.  It  can  be  seen  that  the 
accuracy  is  sufficient  for  practical  applications.  Refer  to  Tables  IV  and  V 
In  PL- 8 
N  =  a 
Z  =  T 

G0  =  Q  (o,t) 
o 

G1  =  y(a, t) 

G2  =  Q  (a, t) 

G5  =  T(a) 
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GAMMA  FUNCTIONS  COMPUTATION  PROGRAM 
(BASIC  LANGUAGE) 

10  REM  G0  IS  THE  NORMALIZED  INCOMPLETE  GAMMA  FUNCTION 
20  REM  iWITH  INTEGRATION  LIMITS  TAU  TO  INFINITY  AND 
30  REM  WITH  ARGUMENT  CH I -SQUARED > NU 
40  REM 

50  REM  ALPHA  MUST  BE  >=  UNITY. 

60  REM  ALPHA=  N 
70  REM  TAU=  Z 
80  REM 

90  DEF  FNG  <Nt  Z) 

100  IF  N>1  THEN  130 
110  Q0=EXP ( - Z ) 

120  GO  TO  480 
130  H  =  Z/6 
140  A=0 
150  M=N 

160  L=  INT ( (N-4) /6) 

170  FOR  1=2  TO  M 
180  Y  =G= 1 

190  FOR  K= 1  TO  1-1 
200  Y - (K-Nt /Z*Y 
210  V1=ABS(Y) 

220  IF  V 1 < 1 E - 30  THEN  240 

230  NEXT  K 

240  A  =  A  +  Y*  (  -  1 )  •*  ( I  +  l ) 

250  NEXT  I 
260  Q=A+1 

270  IF  N<10  THEN  430 
280  J=0 

290  FOR  B= 1  TO  N-l 
300  IF  GC1E34  THEN  330 
310  J=J+l 
320  G=G/ 10 
.330  G  =  Z/B*G 
340  IF  B=L  THEN  400 
350  IF  B  =  2*L  THEN  400 
360  IF  B=3*L  THEN  400 
370  IF  B  =  4*L  THEN  400 
380  IF  B=5*L  THEN  400 
390  GO  TO  410 
400  G=G*EXP(-H) 

410  NEXT  B 
420  GO  TO  470 
430  FOR  B= 1  TO  N-l 
440  G=Z/B*G 
45#  NEXT  B 


t 


pl-  a 


GAMMA  FUNCTIONS  COMPUTATION  PROGRAM 
(BASIC  LANGUAGE) 

CONTINUED 


G  =  G*EXF‘  ( -H*5 ) 

Q  0  =  G  *  Q  *  1  0  M  J  /  2  >  *  £  X  F‘  <  -  H  >  *  1 0  M  J  /  2  > 

FNG=Q0 
FNEND 
N  =  1 2 . 5 
Z= 12 . 5 
W  =  N- I  NT (N) 

IF  W>0  THEN  560 
G0=FNG <N  >  Z ) 

CO  TO  640 
L 6  =  .  5 *  (U-t  '2-W) 

L7  =  1  -  W  t-2 
LO  =  .  5  *  ( W  *  2  -*•  W ) 

N4  =  I NT <N) -  1 
N5= INT  <N) 

N6= I NT (N )  + 1 

G0  =  FNG ( N  4  >  Z ) *L6+FNG <N5 >  Z  > *L ,7  +  FNG <  N6 , Z  >  *L8 

GO  TO  640 

PRINT  "ALPHA*"# N 

PRINT  "TAU= " ! Z 


PRINT  "ALPHA*"# N 
PRINT  "TAU* " #  Z 
PRINT 

PRINT  "NORM’D  GAMMA  FNC*"»G0 
PRINT 
G4= 1 -G0 
R  =  N 
X  =N 
R  =  R+ 1 
X  =  X*R 
D  =  R  -  9 

IF  D  >  =  &  THEN  770 

GO  TO  720 

b  * . 9 1 89385362 

S=S+ (R+ .  5 ) *LOG <R )  -  R 

V  =  1- <  1/  <30*Rt‘2) )  + ( 1 / < 105*R>4: 

V*  <  1  /  ( 12*R)  ) *0 

S  =  S-*- V 

G5  =  EXP  <8)7 X 

PRINT  "GAMMA (ALPHA) =" #G5 

PRINT 

G1 =G5*G4 

PRINT  "FIRST  GAMMA  FNC="#G1 
G2=G5*G0 

PRINT  "SECOND  GAMMA  FNC*"#G2 
END 


<  105*RM>  ) 
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TABLE  IV 

COMPUTATION  EXAMPLES 
(a  *  non-integer) 


ALPHA=  14.5 
T  AU  =  18 

NORM'D  GAMMA  FNC  =  .1735187  (ACTUAL  .17356) 

GAMMA (ALPHA) =  2.30923e  10 

FIRST  GAMMA  FNC  =  1.90854e  10 
SECOND  GAMMA  FNC=  4.00695e  09 


read  y 

*500  N=5.64 
*510  Z=8 
*RUN 


ALPHA=  5.84 
TAU=  8 

NORM'D  GAMMA  FNC=  .1543015  (ACTUAL  .15429 

GAMMA(ALPHA) =  65.71338 

FIRST  GAMMA  FNC  =  55.57371 

SECOND  GAMMA  FNC=  10.13967 


ALPHA  =  2.5 

TAU=  2 

NORM'D  GAMMA  FNC=  .5413411  (ACTUAL  .54924) 

GAMMA (ALPHA) =  1.32934 

FIRST  GAMMA  FNC=  .6097138 
SECOND  GAMMA  PNC*  .7196266 


TABLE  ? 

COMPUTATION  EXAMPLES 
(a  “  integer) 


-flLPHflg - r 

TAU  =  1 

NORM  ’  D  GAMMA  FNC  =  .3678794  (ACTUAL  -36788) 

GAMMA (ALPHA) =  .9999999 

FIRST  GAMMA  FNC=  .6321205 
-"SECOND  GAMMA  FNC  = - .3~S73794 - - - 


ALPHA =  Z 
T  AU=  1.5 

NORM  ’  D  GAMMA  FNC=  .5578254  (ACTUAL  .55783) 

GAMMA (ALPHA) =  .9999999 

FIRST  GAMMA  FNC=  .4421746 
SECOND  GAMMA  FNC=  .5578254 


NORM'D  GAMMA  FNC=  .5987138  (ACTUAL  .59871) 


GAMMA (ALPHA) =  5040 


FIRST  GAMMA  FNC=  2022.482 
SECQNO  GAMMA  FNC  =  - ~  ZmtttT - 


ALPHA=  13 
T  AU*  17.5 

NORM'D  GAMMA  FNC=  .  1 1 1 6488  (ACTUAL  .11165) 

GAMMA (ALPHA) =  4. 79001e  08 

FIRST  GAMMA  FNC  =  4. 25521a  08 
SECOND  GAMMA  FNC=  5.34800a  07 
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D.  ERROR  FUNCTION  COMPUTATION 


Def initlon 

The  Error  Function  used  is  that  of  Papoulis  [_2^  : 


Cu  ■  2 

v  / 


(  =* .■  - 


/  p  rr 


erf  (0)  =  0  ;  erf  (“>)  =0.5 

P  P 


erfp  (-a  )  =  -  erfp  (a) 


Hereafter  the  above  subscript  p  will  be  dropped  for  convenience. 
It  will  be  understood  that  erf  (a)  =  erfp  (ot) . 


Figure  5  shows  a  graph  of  erf  (a)  versus  a. 


TI-59  METHOD 


The  TI-59  Programmable  Calculator  can  be  used  as  the  program  listing  PL-9 
to  calculate  erf (a). 

EXAMPLE: 

FIND:  erf  (1.5) 

SOLUTION:  Use  PL-9  Display 

Enter  1.5  1.5 

Press  jT]  0.4331927713 

erf (1.5)  =  0.4331927713 

EXAMPLE: 

Find:  erf(-0.125) 

Display 

Enter  .125  .125 


Press 


Press 


erf (-0. 125)  =  -0.0497382658 


-.125 


-.0497382658 


PL-9  (Cont’ 


ERROR  FUNCTION  COMPUTATION  PROGRAM 
(TI-59) 


LOC 

CODE  KEY 

LOC 

CODE 

KEY 

084 

65 

1 

26 

0 1 

01 

n  sf. 

fl  1 

i 

1 

54 

) 

fl  R  f. 

Cf  *"* 

1 

il  O 

42 

STD 

n  r  ? 

07 

7 

1 

i 

cL  9 

04 

04 

llRR 

fl  o 

o 

i 

30 

o  7 

IFF 

n  r  r 

fi  1 

i 

i 

3  I 

01 

01 

090 

04 

4 

i 

33 

M  1 

01 

09  i 

07 

7 

i 

•  "«  •  "s 

4 1 

4  i 

fl  9  2 

i  1 

- 

1 

34 

53 

( 

n  9  r 

fjQ 

Q 

1 

•“  cr 

94 

+  /  - 

fl  9  4 

U 

1 

36 

Q  cr 

+ 

fl95 

07 

7 

1 

”7 

9  3 

0  9  f. 

*7  c; 

_ 

i 

i— i 
i_> 

05 

097 

43 

RCL 

1 

39 

54 

) 

0  9  R 

no 

fi  £ 

1 

40 

91 

R/S 

fi  9  9 

65 

i 

41 

53 

i  n  o 

9 

1 

42 

43 

RCL 

101 

0  3 

o 

1 

43 

04 

04 

102 

05 

i 

44 

75 

- 

103 

06 

b 

l 

45 

9  3 

■ 

104 

05 

ET 

i 

46 

05 

r 

105 

0  6 

f: 

l 

47 

54 

106 

03 

•  j 

1 

48 

91 

R/S 

1 07 

07 

7 

108 

08 

o 

109 

02 

1 10 

per 

-r 

1 1 1 

9  3 

e 

112 

03 

•“» 

113 

01 

114 

09 

9 

1  15 

03 

•~! 

1  16 

08 

R 

1  17 

01 

1 

118 

05 

.J 

119 

03 

-“i 

120 

54 

> 

121 

65 

122 

43 

RCL 

123 

02 

02 

124 

65 

125 

43 

RCL 

64 


HP-34C  METHOD 


The  program  listing  for  erf  (a)  is  given  in  PL-15  and  will  not  be 
shown  here.  The  program  is  used  as  follows  to  calculate  erf  (a): 


HP-34C 


ERF (a) 


GIVEN:  Positive  or  negative  value  of  a 
FIND:  ERF (a) 

Put  the  value  of  a  into  X-register  (i.e.  Display  register) 

Key  in  desired  accuracy  * 

Press  B 

Final  display  will  be  the  value  of  erf(a) 

*Accuracy:  In  scientific  notation  set  the  display  for  the  value  of 

a  to  the  desired  number  of  significant  digits  desired  for 
the  final  value  of  erf (a). 


EXAMPLE: 

GIVEN:  o=l. 

FIND:  ERF (a) 

Put  "1"  in  X-register 

Press  f  SCI  7  (for  accuracy  of  7  significant  figures) 
Press  B 


ERF(a)  =  0.3413447  will  be  displayed 


HP -6 7  METHOD 


A.  MATH  PAC  1 

An  Error  Function  Program  is  most  conveniently  obtained  as  program  18-02 
of  Hewlett-Packard's  Math  Pac  1  package.  The  Error  Function  of  program 
18-02  is  defined  slightly  different  than  Papoulis'  Error  Function. 
Program  18-02  can  still  be  used  to  calculate  erf^  (a)  as  follows: 

a)  Use  a  divided  by  /2~  instead  of  a  as  the  argument  value. 

b)  Divide  the  final  answer  by  2. 

B.  PL-10 

Program  listing  PL-10  computes  both  erf  (a)  and  its  inverse. 

TABULAR  METHODS 

Table  VI  tabulates  values  of  erf  (a)  versus  a  for  values  of  a  from  zero 
to  4.0.  These  values  were  computed  using  the  BASIC  LANGUAGE  Program  as 
listed  in  PL-11. 

EXAMPLE: 

FIND  erf  (2.05) 

erf (2 . 05)  =  0.4798178 

FIND  erf  (-1.72) 

erf (1. 72)  =  0.4572838 
erf (-1.72)  =  -0.4572838 

FIND  erf  (2.044) 

erf (2.045)  =  0.4795726 

erf (2.040)  =  0.4793249 

Using  linear  interpolation: 

erf (2.044)  =  0.4795230 

erf (2.044)  ,  =  0.4795232 

actual 
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TABULAR  METHOD  (Cont'd) 

FIND:  erf  (0.068) 

erf (0.065)  =  0.0259131 
erf (0.070)  =  0.0279032 
Using  linear  interpolation 
erf (0.068)  =  0.0271071 

actual 


erf (0.068) 


=  0.0271072 


User  InsiriK'lions 


PL -10 


Program  Description 


Program  Tliie  Papoulis  Error  Function  and  Its  Inverse  -  HP-67 


Name 

Address 

City 


Program  Description,  Equations,  Variables,  etc. 
Evaluates  erf (a)  = 


State 


dx 


Date 

Zip  Code 


As  well  as  the  inverse  function. 


Operating  Limits  and  Warnings  The  computing  time  for  the  inverse  error  function  is  long. 
It  is  a  function  of  the  maximum  error  e  desired. 

For  an  error  <.0005  the  running  time  is  about  10  minutes. 

Range  of  a  for  inverse  error  function  -5  -*■  +5  . 


ERROR  FUNCTION  PROGRAM  LISTING 
(BASIC  LANGUAGE) 


10  DEF  FNE(X> 

20  F'=  .  23 164  I  9 
30  T  =  1  /  ( 1  +1  '*  X  > 

•10  El-  .31938153 

■>0  EE  =  -  .  356563782 

60  E 2=1 .781477937 

70  F; 4  =  -  1  . 82 1255978 

:  0  B5= 1 .  3 30 27  4  4 2 9 

70  H- . 79894228 

t  00  2  -  H  '  E  XP  <  -  ■  5*X  +  2 ) 

i  *  0  -  E 1  ■»  T  AE2*T  t'Z  +  B3  »T  r-3 

j  ;.0  s  =  S  +  E  4  *  T  t  4  +  B  5  *  T  1 5 

120  Ci r  Z  4  3 
i40  FNE- .5  Q 
150  f NEND 

U0  PRINT  "  A"  »  "ERF  (A)" 
170  FOR  X  = . 005  TO  1  STEP 
;  :J0  F  =FNF  <  X ) 

190  PRINT  XfF 
200  NEXT  X 
210  END 


.005 


TABLE  VI  ERROR  FUNCTION  VERSUS  a 


t 


a 

ERF (a) 

a 

ERF  (a) 

0 

.  005 

b 

.0019947 

.0928956 

.01 

.  0039394 

.24 

.0948348 

.  0 1 5 

. 005984 

.  245 

.0967717 

. 

.0079784 

■  4-  -«5 

.0987063 

.  @2:5 

.  0099726 

.255 

.  100633  4 

.03 

. 0 1 1 9665 

.26 

.  1 0  2  5  6  8 

.035 

.04 

.  0 1 3  9  6  0  2 

.  265 

.  i  044952 

.0159535 

.27 

. 1064193 

.  045 

.0179464 

.  275 

.  1083418 

.05 

.0199389 

« 'o 

.  1102612 

.  055 

.0219308 

i -i  cr 

m  <2  O  Jr 

.  1 121779 

.  03 

. 0239222 

.29 

.  1 1 409 1 8 

.  065 

. 0259 1 3 1 

.  295 

.  1  1 6003 

.07 

.0279032 

.  1179113 

.  075 

. 0298927 

.  305 

. 1 198168 

.  06 

.0318814 

.31 

.  1217194 

055 

. 0 338694 

.315 

. 1236191 

.09 

.0358565 

.  Sl. 

.  1255158 

.  095 

.0378427 

cr 

•  •!' J 

. 1274094 

1 

.0398279 

. 

.1292999 

.  10  5 

.0418122 

.  1311874 

.  1  1 

.0437954 

.34 

. 1330717 

.  1 1 5 

.0457775 

.345 

•  1  4  V  5  L. 

.  I . 

. 0477585 

.65 

.  1 368306 

.0497383 

cr  cr 

.  1  38705 1 

.  1  3 

.05171 68 

.36 

.1405764 

.  125 

.0536941 

.  365 

.  1424442 

.  14 

. 05567 

-i  -7 

.  •  / 

. 1443087 

.  145 

. 057 6  446 

.  37  5 

.  1461697 

.  15 

. 0596 1 77 

.  1 4  y  0  l  7  ti 

.  155 

.0615893 

”!!  !,=i 

.1498812 

.  16 

.0635595 

.39 

.  1517317 

.  165 

. 065528 

.  395 

. 1 535785 

.  1  7 

.067495 

.4 

. 1554217 

.  175 

.0694602 

.  405 

. 1572612 

.  IS 

.  0  7  i  4  L  i  7 

.41 

. 1 59097 

.  185 

.0733855 

.415 

.  160929 

.  1 9 

.0753454 

.42 

. 1627572 

.  195 

.0773035 

.425 

.1645816 

.0792597 

.43 

.1664021 

.  205 

.0812139 

.  435 

. 1682188 

.21 

. 083 1661 

.44 

.1700314 

.  2  1 5 

.0851163 

.  445 

.1718401 

.0870644 

.45 

. 1736447 

.  225 

.0890103 

.  455 

. 1754454 

•  o 

.090954 1 

.46 

.1772419 
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TABLE  VI  ERROR  FUNCTION  VERSUS 


a 


a 

ERF  (a) 

46‘:. 

. 1790343 

.47 

.1808225 

475 

. 1826065 

A 

•  "i  O 

. 1843863 

4:75 

. 1 86 1 6 1 8 

.49 

. 187933 

495 

.1396999 

c 
.  J 

. 1914624 

505 

. 1932206 

•  b  1 

. 1949743 

5.15 

.  1967235 

C  •".* 

•  -I  £- 

.1984632 

CT  KT 

•  J  /.  J 

.2002084 

i:r 

« 

.2019441 

cr  ■■  j  cr 

.2036751 

.54 

.2054015 

545 

.2071233 

cr  cr 

-  •  J 

. 2088403 

cr  cr  cr 

j  j  j 

.2105527 

.  b  6 

.2122603 

565 

.2139631 

.57 

.2156612 

575 

.2173544 

.  58 

.2190427 

5:  :5 

.2207262 

.59 

.2224047 

595 

.2240783 

.6 

.2257469 

605 

. 227 4 1 06 

.61 

.  2290691 

615 

. 2307227 

.62 

.2323712 

625 

.2340145 

.2356528 

a  J 

635 

.2372859 

.64 

■  L'i  '!■  V  1  J  b 

645 

.2405365 

.65 

.242154 

655 

.66 

.2453731 

665 

.2469748 

.67 

.2485712 

675 

.2501622 

.63 

.2517478 

.  685 

.2533281 

a_ 

ERF  (a) 

.69 

.  254903 

.  695 

.2564724 

.  7 

.2580364 

705 

.259595 

.71 

.261148 

715 

. 2626955 

.72 

.2642376 

725 

.2657741 

.  73 

.267305 

735 

■ 2  6  o  8  303 

.74 

.2703501 

745 

.2718642 

.75 

.2733727 

755 

.2748756 

.76 

.2763728 

765 

.2778643 

.77 

.2793501 

775 

.  2  8  0  3  r  0  l 

.  78 

.2823046 

785 

.2837733 

.79 

.2352362 

795 

.2866933 

•  O 

.2381447 

805 

.2895902 

.81 

.29103 

o  i  c; 

•J  1  J 

.2924639 

.  Oil 

.293892 

825 

.2950.142 

O 

a  'J 

.2967307 

835 

.2981412 

.84 

.2995459 

845 

.3009446 

.85 

.3023375 

855 

.3037245 

.  86 

.3051055 

865 

. 3064806 

.87 

.3078498 

,875 

.3092131 

.88 

.3105704 

,885 

.3119217 

.89 

.3132671 

.895 

.3146065 

.9 

.3159399 

.  905 

.3172673 
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TABLE  VI  ERROR  FUNCTION  VERSUS  a 


a 

ERF  (a) 

..  3 1  y5;383 

.  ?  1 

.  9 1 5 

.3199042 

.  /  i- 

.  82 1 2  5  U 

-  -  ir 

. 3225 1 7 1 

.93 

.3238145 

.935 

.3251059 

.94 

.3230912 

.945 

.3276706 

.95 

. 328 v4o  / 

.  955 

.3302112 

a  7  O 

.'77 

.3339768 

.  975 

.335219? 

■  '■■/*  O 

.3364569 

.  337688 

.  3 9  1  £? 

.3401 319 

.3413447 

1 . 005 

.3425516 

1.01 

.3437523 

1.015 

.3449471 

1  .02 

.3461358 

3  . 025 

.3473184 

1  .05 

.  348495 

1 .035 

.3496655 

1  .04 

.35083 

1  .045 

. 35 1 9385 

1 . 05 

.3531409 

1 . 055 

.3542873 

1  .06 

.3554277 

1 . 065 

.356562 

1  .07 

.3576903 

1 . 075 

.  35*3 13 1  i- & 

1.08 

.3599289 

1 . 035 

.  36 1 039 1 

1 . 09 

.  3  6  2  1  4  3  4 

1.095 

.3632416 

1  .  1 

.3643339 

1  .  105 

.3654201 

1.11 

.  36650(34 

1.115 

.3675747 

1  .  12 

.368643 

1  .  125 

.3697054 

1.13 

.3707613 

_a 

ERF  (a) 

1  .  135 

. 37 1 8 1 23 

1.14 

.3728568 

1  .  1  45 

.3738955 

1  .  1 5 

.  374  928 

1  .  155 

.8759547 

1.16 

.3769755 

1  .  165 

. 3779904 

1.17 

.3789994 

1  . 175 

. 3800026 

1  .  1 8 

.3809993 

1  .  185 

. 38 1 9 9 1 2 

1  .  1 9 

. 3829767 

1  .  195 

.3839564 

1 .2 

.3849302 

1.205 

•  -2'  3  5  3  7  O  2' 

1.21 

. 3868605 

1.215 

.3878169 

1  • 

.3887675 

1 . 225 

.3897123 

1  .23 

.3906513 

1  •  2o5 

.3915846 

1  .24 

.3925122 

1.245 

.393434 

1  . 2  5 

.3943501 

1  . 255 

.  3952605 

1  . 26 

.3961652 

1  .  265 

.3970642 

1  . 27 

.3979576 

1.275 

. 3988453 

5  .  <-  o 

.  3997273 

1 . 285 

.4006037 

1 .29 

.4014746 

1 . 295 

.4023398 

1 .3 

.4031994 

1  . 305 

. 4040535 

1  .31 

.404902 

1 .3.15 

. 4057449 

1  .32 

.4065824 

1 . 325 

.4074143 

1  .33 

.4082408 

1 . 335 

.4090617 

1.34 

.4098772 

1.345 

.4106873 

1  .35 

.4114919 

1 . 355 

.4122911 

1  •  2>  O 

.4130849 

1 . 365 

.4138734 
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TABLE  VI  ERROR  FUNCTION  VERSUS  a 


a_ 

ERF  (a) 

1  .  .37 

. 4146564 

I  .  3  AS 

.  1 1 5434 Z 

1  .  30 

.4 i 62066 

1  .  385 

.4 1697 07 

1  .  ;*'? 

. 4 1 77354 

1 , 395 

.4 1 0497 

1 . 4 

. 41924 32 

1  .  405 

. 4 1 V VO Vo 

1  .  T  1 

. 120730 1 

1  .415 

. 42 1 4656 

1  .  1 2 

. 4221961 

1  .  425 

.42292  1  3 

1  .43 

. 42364  1  4 

1  .  435 

.4243563 

1 . 44 

.4250662 

1  .  4  45 

.425771 

i  .  4  5 

. 4264706 

1  .  435 

.4271653 

1 .46 

.4278549 

1  .  465 

.  4235:374 

1  . 4  7 

.429219 

1  .4  75 

.4290936 

1  .40 

.4305633 

1  . 405 

. 431228 

I  .  49 

.4318878 

t  .4  95 

.4325427 

1 . 5 

.4331927 

1 .505 

. 4338379 

1 . 6 1 

.4344783 

!  .M  5 

. 435 1 1 38 

t  .52 

.4357445 

1 .525 

.4363704 

1  ..  5  •’ 

. 4369916 

1  .  5  35 

.4876031 

1.34 

.4382198 

1 . 545 

.4383268 

1  . 55 

.4394292 

1 . 555 

.4400269 

1 . 56 

.  44062 

1 . 565 

.4412085 

1  .57 

.4417924 

1 . 575 

.4423713 

1 . 50 

.4429466 

1 . 585 

.4435168 

1  . 59 

.4440826 

1 . 595 

.4446439 

ERF  (a) 


1  .. 

..  4452007 

1  .  605 

.4457531 

1  .61 

. 146301 1 

1 . 6  1  5 

.446 3446 

1  .62 

.4473839 

1 . 625 

.4479187 

1  .68 

.4484493 

1 . 685 

.4489755 

1.64 

.4494974 

t.645 

.4500151 

1 . 65 

. 4505285 

1.655 

.4510378 

1  .66 

.4515428 

1 . 665 

. 4520436 

1  .67 

.4525403 

1.675 

.4530329 

1  f4  O 

.4535214 

1 . 685 

.4540057 

1  .69 

, 454486 

1 . 695 

.4549623 

1.7 

.4554345 

1  .  705 

. 4559028 

1.71 

.4563671 

1.715 

.4568274 

1  . 72 

.4572838 

1  .  725 

.4577363 

1  .73 

.4581349 

1 . 735 

.4536296 

1 .74 

.4590705 

1.745 

.  4595076 

1 . 75 

.4599409 

1 . 755 

.4603704 

1 .76 

.4607961 

1 . 765 

.4612181 

1.77 

.4616365 

1.775 

.4620511 

1  .  78 

.462462 

1 . 785 

.4628694 

1  . 79 

.4632731 

1 . 795 

.4636732 

1  .8 

.4640697 

1 . 805 

.4644627 

1  .81 

.4643521 

1.815 

.4652381 

1  . 82 

.4656205 

1 . 825 

.4659995 

1.83 

.4663751 
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a 

Mi 

ffl 


1 . 835 

1  .84 

1.845 

< 

i.  ■ 

1 . 855 
1  . 86 
l  .  865 
1  .::7 
1.875 

<  r  n 
‘  •  *  '•  O 

1 . 885 

1  . 89 

1  . 895 
1  . 7 

1 . 905 
1  . 9 1 

1.91  5 
1 . 92 
1 . 925 
1  .93 
1  . 935 
1  .94 
1 .945 
1  . 95 
1 . 955 
1  .96 
1 . 965 
1  .97 
1.975 
1  .98 
1 . 985 


TABLE  VI  ERROR  FUNCTION  VERSUS  a 

ERF  (a) 

.4667472 
.4671159 
.4674813 
. 4678433 
.4682019 
.4685573 
.  4689093 
.4692581 
.  4  696037 
.469946 
.4702851 
.470621 1 
.4709538 
. 4712835 
.47161 
.4719334 
. 4  72  2538 
’.472  571  1 
.4728054 
.4731966 
.  4 735049 
.4733102 
.4741 126 
.474412 
.4747085 
.4750021 
.4752929 
.  4755809 
.475866 
. 476 1 483 
.  4764  278 
.4 767046 


i  .7  ' 

l  .  995 

. 4769786 

ii 

.4772499 

2 . 005 

.4775185 

2 . 0  1 

.4777845 

2.015 

.4780477 

2.02 

.4783084 

2 . 025 

.4785664 

2.03 

. 4783218 

2.035 

.4790746 

2.04 

.4793249 

2.045 

.4795726 

2.05 

.4798178 

2 . 055 

.4800606 

2.06 

.4303008 

2 . 065 

.  4805386 

2 . 07 

.4807739 

a 

ERF  (a) 

2 . 0  75 

.4810068 

2.08 

.4812373 

2 . 085 

.4814654 

2.09 

.4816912 

2 . 095 

.4819146 

2. 1 

. 4821356 

2 .  105 

.4823544 

2.11 

.4825709 

2.115 

.4827851 

2.12 

.482997 

2  .  125 

.4832067 

2.13 

.4834142 

2 . 135 

. 4 836 1 95 

2.14 

.4038227 

2. 145 

.4840236 

2.15 

.4842224 

2 . 1 55 

. 48441 9 1 

2 . 1  6 

.4846137 

i  /  rr 

4-  .  1  J 

.4848062 

2.17 

.4849966 

2 . 175 

.435135 

7 . 1  8 

.4853713 

2 . 185 

.4855556 

2 . 1 9 

.4857379 

2 . 195 

.4859182 

i-  m  L. 

. 4860966 

2 . 205 

.486273 

2 . 2  1 

.4864475 

2.21 5 

.48662 

i-  m  i~ 

.4867907 

L  .  'L  v.  b 

. 4869594 

L.  ■  -J 

.4871263 

ii.miL  -J 

.4872914 

L  •  ti  4 

.4874546 

2 . 245 

.487616 

Z  .  Z*5 

.4877756 

2 . 255 

.4879334 

2 . 2  6 

.4880894 

2 . 265 

.4882437 

2.2:7 

.4883962 

2.275 

.4885471 

2.28 

.4886962 

2 . 285 

.4888436 

2.29 

.4889894 

2 . 295 

.4891335 

•7  -‘“i 

.4892759 

2.305 

.4894167 

2.31 

.4895559 

2.315 

.4896935 

v.f  -  « 
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TABLE  VI  ERROR  FUNCTION  VERSUS  a 


a 

ERF  (a) 

2 

.4898296 

2 . 32  S 

.  439964 

i  •  •  '.*  j 

.4900969 

.  4902283 

*"  2 .  34 

.4903581 

2  .  345 

. 4904865 

v- .  35 

.4906133 

2 . 3‘o'5 

. 4907387 

/ 

L  •  O 

.4908625 

fT'  5 

.490985 

. 49 1 1 06 

-  7  cr 

.4912255 

L,  m  OO 

.4913437 

2 . 305 

. 4914604 

2  .09 

.4915758 

•r  “i  “i  i:r 

£  «  -2*  7  •  J 

.4916898 

2.4 

.4918025 

2  .  405 

.4919138 

2.11 

.4920237 

2.4:15 

.4921324 

2.42 

. 4922397 

2 . 425 

.4923458 

l.  .  '4  •  :* 

.4924506 

2.435 

.4925541 

2.44 

■  .4926564 

2.445 

.4927574 

2.45 

.4928572 

2 . 455 

. 4929558 

2.46 

.4930531 

2 . 465 

.4981493 

2.47 

.4932443 

2 . 475 

.4933382 

2.43 

.4934309 

2 . 435 

.4935224 

2.47 

.4936128 

2 . 4  95 

. 493702 1 

2.5 

.4937903 

2 . 505 

.4938774 

2.51 

.4939634 

2.515 

.4940484 

rr  •“ 

£-  a 

.4941322 

-•  c  cr 

u  -  •-><-.  J 

.4942151 

Z  -  5  3 

.4942969 

2 . 535 

.4943776 

2.54 

.4944574 

2.545 

.4945361 

2.55 

.4946138 

2 . 555 

.4946906 

2 . 56 

.4947664 

ERF  (a) 


L.  .  JO  *5 

.4943412 

i.  -  -!j7 

.494915 

■7  cr  “7  rr 

L  .  •  /  J 

.494988 

2.53 

. 49506 

cv.ir 
i-  .  JOJ 

. 495131 

Z  -  5  9 

.4952012 

2 . 5*?b 

.4957704 

2 . 6 

.  4953388 

2 . 605 

.4954062 

2 . 6  1 

.4954728 

2 . 6  1  5 

.4955386 

ii  «  <.J  ii 

.4956035 

2 . 62  5 

.4956675 

Z  -  6  o 

.4957307 

•7*  /  .  rr 

J 

.4957931 

2.64 

.4958547 

2.645 

.4959154 

t  cr 
c  •.  •  J 

.4959754 

2 . 655 

.4960345 

Z  •  /:<6 

.4960929 

2 . 665 

.4961505 

2.67 

.4962074 

2 . 675 

.4962635 

2  . 68 

r  4-  v  £*  1  y  y 

2 . 685 

.  4963735 

2 . 6  V 

.4964273 

2 . 695 

.4964805 

2 . 7 

.496533 

•;  .  705 

.4965847 

2.71 

.  4  9  6  6  3  5  8 

2.715 

.4966362 

*2  m  il. 

.4967358 

~T  •  *.  cr 

<2  a  /  fc-.  J 

.4967849 

2.73 

.4968332 

2 . 735 

.4968809 

2.74 

.496928 

2.745 

.  4969744 

2.75 

. 4970202 

2 . 755 

.4970653 

2 . 76 

.4971099 

2 . 765 

.4971538 

2.77 

.4971971 

2.775 

.4972399 

2.78 

.497282 

2 . 785 

.4973236 

2.79 

.4973645 

2 . 795 

.497405 

O 

C  «  'J 

.4974448 

2.805 

.4974841 
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TABLE  VI  ERROR  FUNCTION  VERSUS  a 


a 

ERF  (a) 

2.81 

.4975229 

2.815 

.4975611 

2.82 

.4975988 

2.825 

.4976359 

2.83 

.4976725 

2 . 835 

.4977086 

2.84 

.4977443 

2.845 

.4977794 

2.85 

.497814 

2.855 

.4978481 

2 .36 

.4978817 

2.865 

.4979149 

2.87 

.4979476 

2.875 

.4979798 

2.83 

.4980116 

2.835 

.4980429 

2.39 

.4980737 

2 . 895 

.4981041 

2.9 

.4981341 

2 . 9(35 

.4981637 

2.91 

.4981928 

2.915 

.4982215 

2.92 

.4982493 

2.925 

.4982776 

2.93 

.4983051 

2.935 

.4983322 

2.94 

.4983589 

2.945 

.4983852 

2.95 

.4984111 

2.955 

.4984366 

2.96 

.4984617 

2 . 965 

.4984865 

2.97 

.4985109 

2.975 

.498535 

2.98 

.4985587 

2.985 

.498582 

2.99 

.498605 

2.995 

.4986277 

3 

.49865 

3.005 

.498672 

3.01 

.4986937 

3.015 

.498715 

3.02 

.4987361 

3.025 

.4987568 

3.03 

.4987772 

3.035 

.4987972 

3.04 

.498817 

3.045 

.4988365 

3.05 

.4988557 

0 

ERF  (a) 

3.055 

.4988746 

3.06 

.4988932 

3.065 

.4989116 

3.07 

.4989296 

3 . 075 

.4989474 

3.08 

.4989649 

3.035 

.4989822 

3.09 

.4989991 

3 . 095 

.4990159 

3.1 

.4990323 

3.105 

.4990485 

3.11 

.4990645 

3.115 

.4990802 

3.12 

.4990957 

3.  125 

.4991109 

3.13 

.4991259 

3.135 

.4991407 

3.14 

.4991552 

3.145 

.4991695 

3.15 

.4991836 

3.155 

.4991974 

3.16 

.4992111 

3.165 

.4992245 

3.17 

.4992377 

3.175 

.4992508 

3.18 

.4992636 

3.135 

.4992762 

3.19 

.4992886 

3.195 

.4993008 

3.2 

.4993128 

3 . 205 

.4993246 

3.21 

.4993363 

3.215 

.4993477 

3.22 

.499359 

3.225 

.4993701 

3.23 

.499381 

3.235 

.4993917 

3.24 

.4994023 

3.245 

.4994127 

3.25 

.4994229 

3.255 

.499433 

3.26 

.4994429 

3.265 

.4994526 

3.27 

.4994622 

3.275 

.4994716 

3.28 

.4994809 

3.285 

.49949 

3.29 

.499499 

3.295 

.4995078 
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TABLE  VI  ERROR  FUNCTION  VERSUS  a 


a 

ERF  (a) 

■Z>  ■  -5 

.4995165 

0  5 

.4995251 

.31 

. 4995385 

313 

.  4995417 

.4995499 

335 

.  4995578 

■  -J 

. 4995657 

•  .  *1  7l 

.4995734 

.34 

. 499581 1 

345 

. 4995835 

>  ."iC 
•  -  -•  -J 

. 4995959 

•~«  i..*  cr 
•»  - 1  *-J 

.. 499603 1 

•  ./ 

. 4996 1 02 

il: 

.  499617.. 

•  ••  • 

>  ,  / 

. 4996241 

,  3  /  j 

.4996309 

.»  .  i“i 

49963  75 

-  •'  i  c 

.  :■  ■  j 

.4996441 

. 4996505 

.  j  '•/ j 

.4996568 

3 . 4 

. 499663 

.  40*3 

.4996691 

.  4 1 

.  4996751 

.4  15 

.4  99';  81 

3 .  43 

. 4996368 

4  .C* 

.  ‘T  i  .  J 

. 49 9 8925 

3.43 

.4996982 

.  43 3 

. 499/037 

3.44 

.4997091 

.44  5 

.  4997144 

3.45 

.4997197 

.  4'V- 

.4997248 

3 .  1 6 

.  4  997299 

.  435 

.4397348 

3 . 4  / 

. 4997397 

.475 

.4997445 

3.48 

.4997492 

•  4  b' 5 

.4997539 

3 . 49 

. 4997584 

.495 

.  4997629 

3.5 

.4997673 

.  505 

.4997717 

3.51 

.4997759 

.515 

.4997801 

3.52 

.4997842 

C*  • cr 
•  Ji_  J 

.4997882 

3.53 

.4997922 

'  •  535 

.4997961 

3.54 

. 4997999 

a_ 

ERF  (a) 

545 

. 4998037 

i  r*  ir 
.  .J  J 

.4998073 

555 

.49981  1 

.56 

.4998145 

565 

.  499818 

;  .57 

.  4  9  V  5  215 

575 

.  4993243 

■ .  58 

.  4  9  v  y  2  y  2 

,  585 

.4998814 

.49  .'8  846 

,  595 

.4998  378 

3 . 6 

..4  993409 

.  605 

.4998487 

3 .  6  1 

.  1  9  8  4  6 

.  6 1 5 

.49 78498 

.4993527 

.  625 

.4998555 

8.63 

/i  -}  o  c.- •  o 

.685 

.  4  V  '7861 

3.64 

. 4 v 986 86 

.  645 

. 4993663 

3.65 

. 4993688 

.  655 

.4998714 

8 . 6  6 

. 4998739 

/. c 
«  L'  u 

. 4998763 

.  /-  ,  7 

.4998787 

.  675 

.  499881 

3.68 

. 4998833 

.  685 

.4998856 

3 . 69 

. 4998878 

.  695 

. 49989 

3 . 7 

.4998922 

.  705 

. 4998943 

3 . 7 1 

. 4993963 

.715 

.4998984 

3.72 

.4999004 

.  725 

.  4999023 

-  r  •- 1 

2*  .  /  O 

. 4999042 

: .  735 

.4999061 

3.74 

.  499908 

1.745 

.4999098 

3.75 

.4999116 

; .  755 

.4999133 

3.76 

.499915 

3 . 765 

.4999167 

3.77 

.4999183 

3.775 

.49992 

3.78 

.4999216 

3 . 785 

.4999231 

3 . 79 

.4999246 

81 
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E.  INVERSE  ERROR  FUNCTION  COMPUTATION 


Def inition 


The  Inverse  Error  Function  is  that  value  of  a  that  yields  a  specified 
value  of  erf  (a).  That  is,  given  the  value  of  erf  (a),  find  a. 

TI-59  METHOD: 

The  following  computation  method  applies  to  positive  values  of  a. 


For  negative  values  of  a  use  the  relationship 
erf  (-a)  =  -  erf  (a) 

That  is,  calculate  erf  (a)  and  change  the  sign  of  erf  (a) 
for  negative  a  values. 


Def ine 


K  =  erf  (a)  ;  0  <  K  <  0.5 


z  *  V-2  Hn(2K) 


^  =  g0  +  giz  +S2 


z  +  . 


Where: 


g o  =  6.55864x10 
gl  =  -0.02069 


=  0.737563 


g3  =  -0.207071 
gA  =  -2.06851x10' 
=  0.03444 


+  g10 


g6  =  -1.17213x10”^ 
g?  =  2. 10941xl0~3 
gg  =  -2.18541x10”^ 
g9  =  1.23163x10”^ 
gin  =  -2.93138xl0-7 


TI-59  METHOD  (Cont'd) 


The  following  listing  PL-12  is  for  the  TI-59  Programmable  Calculator. 

It  utilizes  the  equations  previously  described.  To  use:  enter  the 
value  of  erf  (a) ;  then  press  A  .  The  computed  value  of  a  will 
ultimately  be  displayed.  Table  VI  can  be  used  as  a  check  on  the 
reasonableness  of  the  answer.  Computation  time  is  approximately  10  seconds. 

0  <  erf  (a)  <0.5 

For  negative  values  of  erf  (a),  calculate  a.  for  the  positive  value  of  erf 
(a).  Then  change  the  sign  of  a. 

Example: 

Find  a  for  erf  (a)  =  0.4995 
Enter  0.4995;  Press  j~A  | 

Display:  3.29049 
Thus  a  =  3.29049 


Example: 

Find  a  for  erf  (a)  = 
Enter  0.4975  ;  Press 
Display:  2.8066 
Thus  a  =  -  2.8066 


0.4975 


A 
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PL  -  12 


INVERSE  ERROR  FUNCTION; 

LISTING  FOR  TI-59  PROGRAMMABLE  CALCULATOR 

LOC  CODE  KEY  LOC  CODE  KEY 


044  43  RCL 

000  76  LEL  n45  02  02 

001  11  H  ii46  54  ) 


,  V- 


I 


INVERSE  ERROR  FUNCTION; 

LISTING  FOR  TI-59  PROGRAMMABLE  CALCULATOR 


LOC 

CODE  KEY 

LOC  CODE 

KEY 

OSS 

94  +.■■-- 

1 

-:3 

01 

1 

089 

44  SUM 

1 

34 

Q 

090 

03  03 

1 

35 

01 

1 

091 

53  <; 

1 

36 

07 

*7 

0  9  2 

02  2 

1 

02 

093 

9  3 

1 

”i  o 

01 

1 

094 

00  0 

1 

39 

03 

3 

095 

06  6 

140 

52 

EE 

096 

08  8 

141 

02 

j! 

097 

05  5 

142 

94 

+  /  - 

0  9  8 

0 1  1 

143 

65 

099 

52  EE 

144 

Cj  ■-! 

r 

1 0  0 

02  2 

145 

43 

RCL 

1 0 1 

94  +/- 

146 

02 

02 

102 

65  x 

4 

i 

47 

54 

j 

103 

c;  j  j*‘ 

148 

45 

Y  -:  :; 

104 

43  RCL 

149 

06 

6 

105 

02  02 

1 

50 

54 

j 

106 

54  > 

1 

51 

94 

+  /  — 

107 

45  YX 

1 

52 

44 

SUM 

108 

04  4 

1 

53 

03 

03 

109 

54  ) 

1 

54 

cj  o 

l'* 

110 

94  +/- 

1 

55 

02 

1  1  1 

44  SUM 

1 

56 

93 

H 

112 

03  03 

1 

57 

6i 

1 

113 

53  < 

1 

tr  c> 

00 

0 

114 

03  3 

1 

59 

09 

9 

115 

9  3 

160 

04 

4 

116 

04  4 

161 

01 

1 

117 

04  4 

162 

E7 

■J  C. 

EE 

118 

04  4 

163 

0  3 

3 

119 

52  EE 

164 

94 

+  /  - 

120 

02  2 

165 

65 

121 

94  +/- 

1 66 

53 

( 

122 

65  x 

167 

43 

RCL 

123 

53  < 

168 

02 

02 

124 

43  RCL 

169 

54 

') 

125 

02  02 

170 

45 

Y  :‘:- 

126 

54  > 

171 

07 

«■ 

127 

45  vx 

172 

54 

j 

128 

05  5 

173 

44 

SUM 

129 

54  ;■ 

174 

03 

03 

130 

44  SUM 

175 

53 

131 

03  03 

176 

02 

•L-  -r 

132 

53  <: 

86 


PL-iZ 


INVERSE  ERROR  FUNCTION; 

LISTING  FOR  TI-59  PROGRAMMABLE  CALCULATOR 


LQC 

CODE  KEY 

LQC 

CODE 

KEY 

17? 

Q 

221 

09 

9 

i  O 

U  i  l 

2  2  2 

0  3 

3 

179 

0  8  8 

2  2  3 

01 

1 

1  O  Cl 

ncr  CT 

•~T  .-i 

__ 

i  O  w 

bi  d  4 

IJ  ,~i 

1  o  1 

04  4 

•~i  ■  j  cr 

n  c* 

o 

1  O  bl 

0 1  1 

2  2  6 

52 

EE 

1  83 

52  EE 

71 

07 

1 

1  bH 

U  4  4 

2  2  8 

94 

4-  /  - 

1  ;-!FJ 

Q  4  t/  - 

2  2  9 

65 

1 8  6 

65 

230 

cr  ~i 

i  C'  i 

cr  s-‘ 

2  3  1 

4  3 

RCL 

1  o  o 

4  8  R  U  L 

ci  Z'  2 

02 

02 

1  09 

U  2  U  ui 

2 '  “« 3 

54 

190 

54  > 

2  3  4 

45 

V 

1 9 1 

45  Y>:; 

“I  C. 

01 

1 

192 

OS  8 

2  3  6 

00 

0 

1  9  3; 

54 

2  3  7 

54 

) 

194 

94  +  /  - 

2  3  o 

94 

+  - 

1 95 

44  SUM 

2  3  9 

44 

SUM 

1 96 

03  03 

240 

0  3 

n  3 

1  47 

53 

24  i 

43 

RCL 

i  9 1?. 

0 1  I 

242 

0  3 

0  3 

1 9  9 

9  3  . 

243 

9  1 

R/S 

2  y  u 

02  2 

244 

61 

GT0 

20  i 

0  3  3 

245 

01 

01 

2Q2 

O  l  l 

246 

1 0 

10 

2  U  3 

0  6  6 

247 

O  1 

1 

RST 

4  03 


53  EE 


94  +•■■■'  ~ 

kf;  x 


INVERSE  ERROR  FUNCTION  (HP-67  METHOD) 


Refer  to  PL-10. 


INVERSE  ERROR  FUNCTION  (HP-34C  METHOD) 

The  program  listing  for  the  Inverse  Error  Function  and  other  functions 
for  the  HP-34C  is  given  in  PL-15.  Included  here  is  a  brief  listing  for 
erf  (a)  and  its  inverse  as  well  as  an  example  in  its  use. 


* 

* 

88 
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AD-A103  211  GENERAL  ELECTRIC  CO  UTICA  NY  AEROSPACE  ELECTRONIC  SY--ETC  F/G  20/11 
COMPUTATION  METHODS. (U) 

MAR  80  R  G  LAMBERT  NOOO19-70-C-O4O7 


END 


9 


HP34C 

ERF  (a)  AND  INVERSE  (a) 


erf  (a)  -*■  Rg 


LBL  0 
X2 

2 

CHS 

ex 

RTN 

LBL  I 
GSB  B 
RCL  6 

STO  7 
RTN 

LBL  B 
STO  7 
0 

ENTER  t 
RCL  7 
fO 
2 
v 

X 

r~ 

l/x 

X 

RTN 


89 


HP-34C 

INVERSE  ERROR  FUNCTION 


GIVEN:  ERF(a)  ;  positive  or  negative  value 
FIND:  a 


Put  the  value  of  erf (a)  into  the  X-register 
Key  in  desired  accuracy  * 

Store  the  value  of  erf (a)  into  register  6 


Put  initial  estimate  on  bounds  of  a  .  Upper  and  lower  bounds  must 
bracket  true  value  of  a  .  Lower  bound  is  put  into  X-register.  Then 
press  ENTER  +  .  Then  put  upper  bound  in  X-register.  Press  f 
SOLVE  1. 


Execution  time  is  approximately  as  follows: 

Accuracy 
f  SCI  3 
f  SCI  5 
f  SCI  7 


Time 

7  minutes 
14  minutes 
24  minutes 


EXAMPLE: 

GIVEN:  ERF (a)  =  -.4331927 
FIND:  a 


Put  -. 

4331927  in  X-register 

*  Press 

f  SCI  7  (for  accuracy  of 

7  significant  figures) 

Press 

STO  6 

Press 

2  CHS  (i.e.  -2) 

K  ai  m/i  o  a n  a  fir  mv 

Press 

ENTER  +  y?— 

Dounas  on  answer 

Press 

.1  CHS  (i.e.  -.1 V 

Press 

f  SOLVE  1 

Final  Display  =  -1.500 
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F.  PROBABILITY  OF  FAILURE  COMPUTATION 


TI-59  METHOD 


Many  probability  of  failure  F(a)  expressions  are  of  the  form 
F(«)  =  0.5  +  erf (a) 

This  expression  can  be  conveniently  computed  using  PL-13.  To  use 


enter  the  value  of  a  and  press  A 
ultimately  be  displayed. 


The  value  for  F(rx)  will 


EXAMPLE : 


FIND:  F(a)  for  a  =  -2.305,  -1.365,  -0.68,  -0.365,  -0.125,  0.125, 

0.365,  0.68,  1.365,  2.305 


-2.305 

-1.365 

-0.68 

-0.365 

-0.125 

0.125 

0.365 

0.68 

1.365 

2.305 


0.0105832588 

0.0861265831 

0.24825221581 

0.3575557601 

0.4502617342 

0.5497382658 

0.6424442399 

0.7517478419 

0.9138734169 

0.9894167412 


PL- 13 


PROBABILITY  OF  FAILURE  COMPUTATION 
(TI-59) 


LOC 

CODE 

KEY 

LOC 

CODE 

KEY 

00 

76 

LBL 

042 

0  3 

03 

|-| 

01 

1 1 

H 

043 

35 

i/'N 

j“j 

02 

22 

INV 

044 

53 

j": 

0  3 

86 

STF 

045 

53 

( 

j~i 

04 

01 

Oi 

046 

42 

STD 

i“i 

05 

29 

CP 

047 

02 

02 

j“i 

06 

77 

GE 

048 

45 

i 

i"| 

07 

87 

IFF 

049 

04 

4 

i"| 

OS 

94 

+  /  - 

050 

65 

j~i 

09 

86 

STF 

051 

01 

1 

i“i 

10 

0 1 

0  1 

052 

93 

, 

1 1 

7  6 

LBL 

053 

03 

o 

u 

12 

O  » 

IFF 

054 

0  3 

o 

13 

«=■■". 

•• 

055 

0  0 

0 

14 

42 

STO 

0  5  6 

02 

l“j 

i  CT 

i  J 

0  3 

0  3 

057 

07 

r* 

«"} 

16 

v  2 

053 

04 

4 

u 

17 

22 

INV 

059 

04 

4 

IS 

V ' -i 

LNX 

0  6  0 

02 

i— 

19 

65 

061 

09 

9 

«“l 

20 

02 

062 

75 

- 

(”j 

2  1 

65 

063 

43 

RCL 

l“l 

22 

39 

■ti 

064 

02 

02 

23 

54 

> 

065 

45 

Y 

24 

34 

■i  X 

0  6  6 

0  3 

•»' 

u 

25 

35 

1/X 

067 

65 

u 

2  6 

42 

STO 

068 

01 

i 

t-t 

27 

01 

01 

0  6  9 

93 

• 

93! 

1 J  (  !_i 

0  8 

071 

u 

29 

02 

U2 

i— 

u 

30 

03 

072 

01 

1 

31 

01 

i 

U  <'  3 

02 

d 

l“l 

32 

06 

6 

074 

05 

C* 

•J 

j“j 

■~i 

04 

4 

075 

05 

cr 

l“l 

34 

01 

1 

0  7  6 

09 

9 

u 

35 

09 

9 

U  r  (' 

U  f 

l 

36 

49 

PRD 

078 

08 

r! 

■2*  \ 

0  3 

03 

079 

85 

+ 

l“l 

3  8 

01 

1 

080 

43 

RCL 

39 

44 

SUM 

081 

02 

02 

LI 

40 

03 

03 

082 

45 

Y 

u 

41 

43 

RCL 

083 

02 

92 
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PL-13  (Cont'd) 


PROBABILITY  OF  FAILURE  COMPUTATION 
(TI-59) 


J 


LOC 

CODE 

KEY 

LOC 

CODE 

KEY 

084 

65 

1  26 

01 

01 

085 

01 

t 

127 

54 

'} 

08b 

93 

* 

128 

8 

IFF 

087 

07 

*7 

129 

U  1 

01 

088 

08 

t”s 

1 30 

0 1 

01 

089 

f S 1 

4 

1  3 1 

37 

37 

090 

04 

4 

1  32 

r 

( 

09 1 

07 

I  O  -2‘ 

94 

+  /  - 

092 

07 

“ 

134 

85 

+ 

093 

Li  9 

9 

i  O  cr 
i  J 

0 1 

i 

094 

03 

i  36 

54 

) 

095 

07 

■7 

i  !' 

91 

R/S 

096 

75 

- 

097 

43 

RCL 

098 

02 

02 

099 

65 

1 00 

93 

t 

1 0  1 

0  3 

•*: 

102  05  5 


103 

06 

b 

104 

05 

5 

105 

0  6 

6 

106 

03 

"5 

1 0  { 
i  np 

U  i 

fj  o 

_ 

109 

02 

O 

110 

C;  Cj 

4- 

1 1 1 

9  3 

112 

03 

3 

i  i  ■' 

X  X 

U 1 

i 

i  i  .i 

1 

0  9 

9 

1  15 

03 

■Zl 

1  16 

08 

8 

117 

01 

i 

1  1 8 

05 

c 

.J 

119 

03 

3 

120 

54 

> 

121 

65 

1  22 

‘“I 

“t  j 

RCi 

123 

02 

0 

124 

65 

125 

43 

RCl 

93 
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BASIC  LANGUAGE  METHOD 

PL-14  is  a  modification  of  PL-llto  calculate  F(a). 

F(a)  =0.5+  erf (a) 

TABULAR  METHOD 

Table  VII  tabulates  F(a)  versus  a. 

EXAMPLE: 

FIND:  FC-2.305)  using  Table  VII 
F (-2.3)  =  0.0107241 
F(-2.31)  =  0.010444 
Using  linear  interpolation 
F  (-2 . 305)  =  0.010584 

F(-2'305) ACTUAL  *  °-°10583 

EXAMPLE: 

FIND:  F ( .368)  using  Table  VII 

F( . 36)  =  0.6405764  ;  F(.37)  =  0.6443087 

Using  linear  interpolation:  F(,368)  =  0.6435622 

F ( . 368)  =  0.6435633 

ACTUAL 


i 


V- 


FAILURE  PROBABILITY  PROGRAM  LISTING 
(BASIC  LANGUAGE) 


10  REM  A* ALPHA 
20  DEF  FNE  <  X ) 

30  P*. 2316419 
40  T  =  1  /  { 1 +P*  X ) 

50  El =.3 1938 153 
60  B2  =  -  . 356563782 
70  B3» 1.781477937 
80  B4*-l .821255978 
90  E5*1.330274429 
100  H*. 39894228 
110  Z=H#EXP  < - . 5#  X  1 2 ) 

120  S  =  B  1  ■*T+B2*T»-2+B3*T 1 3 
130  S  =  S  +  B  4  *  T  »•  4 + B  5  #  T 1 5 
140  Q  =  2 * S 
150  FNE= .5  -Q 
160  FNEND 

170  FOR  A*- 2  TO  -1  STEP  .01 
180  IF  A>=0  GO  TO  220 
190  X  =ABS (A) 

200  F  = . 5 -FNE ( X ) 

210  GO  TO  230 
220  F* . 5+FNE (A) 

230  PRINT  A t F 
240  NEXT  A 
250  END 


TABLE  VII  FAILURE  PROBABILITY  VERSUS  a 


a 

FCa) 

-4 

3. 1 6873 e- 05 

-3.99 

3.305452-05 

-3.98 

3.447382-05 

-3.97 

3.595282-05 

-3.96 

3.749132-05 

-3.95 

3.909322-05 

-3.94 

4. 075342 -05 

-3.93 

4.249442-05 

-  3 . 9  2 

4.429372-05 

-3.91 

4.616752-05 

-3.9 

4.811582-05 

-3.89 

5.014242-05 

-3 .83 

5. 22509 e -05 

_  o  o  7 

5.444142-05 

—  6 .36 

5.67175a-05 

-3.85 

5.90831e-05 

-3.84 

6 . 154102-05 

_  “>  O  '"V 

6.409732-05 

—  .  O  i- 

6 . 674902-05 

-  SK.  8 1 

6.951 022 -05 

-  6 . 8 

7.237492-05 

-3.79 

7.535142-05 

-3.78 

7.843972-05 

-3.77 

8.165092-05 

-3.76 

8.498502-05 

-3.75 

8.844582-05 

-3.74 

9.204072-05 

_o  70 

9.576982-05 

_  O  7  O 

w  .  / 

9 . 96441 2-05 

-3.71 

.0001037 

-3.7 

.0001078 

-  o  /.  Q 

.0001122 

-  8  •  6>  8 

.0001167 

-3.67 

. 000 1213 

r_3 .66, 

.0001261 

-3.65 

.0001312 

-3.64 

.0001364 

-3.63 

.0001417 

-  3 . 62 

.0001473 

-3.61 

.0001531 

..-3.,  6  . . 

^3001591 

-3.59 

.0001654 

-3.58 

.0001718 

-3.57 

.0001785 

-3.56 

.0001855 

-3.55 

.0001927 

.-000200 1 

a 

FCa) 

-  O  Ev? 

O  • 

.0002078 

-3.52 

.0002158 

-3.51 

.0002241 

”3-5 

.0002327 

-3.49 

.0002416 

-3.48 

. 0002508 

-3.47 

.0002603 

-3.46 

.0002701 

-3.45 

.0002803 

-3.44 

. 0002909 

-3.43 

.0003018 

-3.42 

.0003132 

-3.41 

.0003249 

-3.4 

.000337 

-3.39 

.0003495 

-3.38 

.0003625 

-3.37 

.0003759 

-3.36 

.0003898 

-3.35 

.0004041 

-3.34 

.0004189 

•~i 

0.0  o 

.0004343 

-  ■  J 

.0004501 

-3.31 

.0004665 

“  O  ■  O 

.0004835 

~3 . 29 

.000501 

-3.28 

.0005191 

*  o  •  7 

.0005378 

-  3 . 26 

.0005571 

-3.25 

.0005771 

-3.24 

.0005977 

_  o  o 

w  «  ii-O 

.000619 

O  ■  <t«  i_ 

.000641 

-3.21 

.0006637 

-3.2 

.0006872 

-3.19 

.0007114 

-3.18 

.0007364 

-3.17 

.0007623 

-3.16 

.0007889 

-3.15 

.0008164 

-3.14 

.0008448 

_-3. 13 

.0008741 

-3.12 

.0009043 

-3.11 

.0009355 

-3.1 

.0009677 

-3.09 

.0010009 

-3.08 

.0010351 

-3.07  . 

.  _ _ J00 10704 

96 
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TABLE  VII  FAILURE  PROBABILITY  VERSUS  a 

4 


FpsO 

.0011068 
.0011443 
.001183 
.0012228 
.001 2639 
.0013063 
. 00135 
. 00 1 395 
.0014413 
.0014891 
.0015383 
.0015889 
.0016411 
.0016949 
.0017502 
.0018072 
.0018659 
,0019263 
.0019884 
.0020524 
.0021133 
.002186 
.0022557 
.  0  0  2 3  2  7 5 
.0024012 
.0024771 
.0025552 
.0026355 
.002718 
. 0028029 
.0028901 
.0029798 
.003072 
.0031668 
.0032641 
.0033642 
.003467 
.0035726 
.0036812 
.0037926 
.0039071 
.0040246 
.0041453 
.0042693 

_  .0043965 

“-2 '61  ~  .0045271 

-2  ,j6  .0046612 


FCfO 

-£759 

.0047988 

-2.58 

.00494 

-2.57 

.005035 

-2.56 

.0052336 

-2.55 

.0053862 

-  2 . 54 

.0055426 

-2.53 

.0057031 

-2.52 

.0058678 

-2.51 

.0060366 

-2.5 

.0062097 

-2.49 

.0063872 

-2.48 

.0065691 

-2.47 

.0067557 

-2.46 

.0069469 

-2.45 

.0071428 

-2.44 

.0073436 

-2.43 

.0075494 

-2.42 

.0077603 

-2.41 

.0079763 

-2.4 

.0081975 

-2.39 

.0084242 

-2.38 

.0086563 

-2.37 

.003894 

-2.36 

.0091375 

-2.35 

.0093867 

-2.34 

.0096418 

***  £  *  3 

.0099031 

~  4.  •  4- 

.0101704 

-  2  1 

.010444 

-2.3 

.0107241 

-2.29 

.0110106 

-2.28 

.0113038 

.  -2.27 

.0116038 

-2.26 

.0119106 

-  2..  25- 

.0122244 

-2.24 

.0125454 

•r-  O 

-  U  • 

.0128737 

~  ’4  >  Z 

.0132093 

-2.21 

.0135525 

-2.2 

.0139034 

-2.19 

.0142621 

-2.18 

.0146287 

-2.17 

.0150034 

-2.16 

•0153S63 

-2.15 

.0157775 

-2.14 

.0161773 

-2.13 

.0165857 

a_ 

-3.06 

-3.05 

-3.04 

-3.03 

-3.02 

-3.01 

-2.99 

-2.98 

-2.97 

-2.96 

-2.95 

-2.94 

-  2 . 93 
-2.92 
-2.91 

-2.9 

-2.39 

-■  o  p 
"  t,  i'ju 

o  7 

-  L.  +'->  J 

-2 . 86 
-2.35 

-  2  •  34 


-  2  •  3 1 

”  L.  • 

-2.79 
-2.78 
-2.77 
-2.76 
-2.75 
-2.74 
-2.73 
-2.72 
-2.71 
-2.7 
-2.69 
-  2 . 68 
-2.67 
-2.66 
-2.65 
-2.64 
-2.63 
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TABLE  VII  FAILURE  PROBABILITY  VERSUS  a 


a 

FCa) 

-2.12 

.0170029 

-2.11 

.0174291 

-2.1 

.0178643 

-2.09 

.0183088 

-2.08 

.0187627 

-2.07 

.0192261 

-2.06 

.0196992 

-2.05 

.0201821 

-2.04 

.0206751 

-2.03 

.0211782 

-2.02 

.0216916 

.-2.01 

.0222155 

.02275 

-  1  .99 

.0232954 

-1.98 

.0238517 

-1.97 

.0244191 

-  1 .96 

.0249978 

-1.95 

.025588 

-  1  .94 

.0261898 

-1.93 

.0268033 

-  1  .92 

.0274289 

-1.91 

.0280665 

-1.9 

.0287165 

-1.89 

.0293789 

-1.88 

.030054 

-1.87 

.0307418 

-  1  . 86 

.0314427 

-1  .85 

.0321567 

-1.84 

.0328841 

-  1  .83 

.0336249 

-1.82 

.0343794 

-1.31 

.0351478 

-1.8 

.0359302 

-  1.79 

.0367269 

-1.78 

.0375379 

-1.77 

.0383635 

-1.76 

.0392038 

-1.75 

.0400591 

-1.74 

.0409294 

-1.73 

.0418151 

-  1.72 

.0427162 

-1.71 

.0436329 

-1.7 

.0445654 

-1.69 

.0455139 

-1.68 

.0464786 

-1.67 

.0474596 

-1.66 

.0484572 

a 

F(a) 

-1.65 

.0494714 

-1.64 

.0505025 

-1.63 

.0515507 

-1.62 

.0526161 

-1.61 

.0536989 

-1.6 

.0547992 

-1.59 

.0559174 

-  1  .58 

.0570534 

-1.57 

.0582075 

-1  .56 

.0593799 

-1.55 

.0605707 

-1.54 

.0617801 

-1.53 

.0630083 

-1.52 

.0642554 

-1.51 

.0655217 

-1.5 

.0668072 

-1.49 

.0681121 

-1.48 

.0694366 

-1.47 

.0707809 

-1.46 

.0721451 

-1.45 

.0735293 

-1.44 

.0749337 

-1.43 

.0763585 

-1.42 

.0778039 

-1.41 

.0792699 

-1.4 

.0807567 

-1.39 

.0822645 

-1.38 

.0837934 

-1.37 

.0853435 

- 1 . 36 

.086915 

-  1 . 35 

.088508 

-1.34 

.0901227 

-  1  . 33 

.0917592 

-  1 . 32 

.0934176 

-  1 . 3 1 

.095098 

-1.3 

.0968005 

-1.29 

.0985254 

-1.28 

. 1002726 

-1.27 

. 1020423 

- 1 .26 

. 1038347 

-1.25 

. 1056498 

-1.24 

. 1074877 

-1.23 

. 1093486 

-1.22 

.1112325 

-1.21 

.1131395 

-1.2 

.1150697 

-1.19 

. 1170232 
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TABLE  VII  FAILURE  PROBABILITY  VERSUS  a 


a  F  (a) 


-1.18 

.1190001 

-1.17 

. 1210005 

-1.16 

.1230244 

-1.15 

.1250719 

-1.14 

. 1271431 

-1.13 

.1292381 

-1  .  12 

. 1313569 

-1.11 

.1334995 

-1.1 

.135666 

-  1 . 09 

. 1378565 

-1.08 

.140071 

-1.07 

. 1423096 

-1.06 

.  1445722 

-1.05 

. 146859 

-1.04 

.  1491699 

-1.03 

. 1515049 

-1.02 

. 1538641 

-1.01 

. 1562476 

- 1 

.  1586553 

-  .99 

.  1610871 

-  .98 

. 1635431 

-  .97 

.  166-02  32 

-  .  96 

. 1685^76 

-  .95 

.1710561 

-  .94 

. 1 736088 

-  .93 

.1761855 

-  .92 

.  1737363 

-.91 

.1814112 

-  .9 

.  1840601 

-  .39 

.  1867329 

- .  yy 

.  1894296 

-  .37 

.  1921501 

-  .36 

. 1943945 

-  .  35 

.1976625 

-  .34 

.2004541 

-  .83 

.2032693 

.206108 

-  .31 

.20897 

-  .  8 

.2118553 

-  .79 

.2147638 

-  .78 

.2176954 

-.77 

.2206499 

-  .76 

.2236272 

-  .75 

.2266273 

-  .74 

.2296499 

-  .73 

.232695 

-  .72 

.2357624 

a_ 

F(ct) 

-  .71 

. 238352 

-.7 

.2419636 

-.69 

.245097 

-  .  68 

.2482521 

-  .67 

.2514288 

-  .66 

.2546268 

-  .65 

.257846 

-  .64 

.2610862 

-.63 

.2643472 

-  .62 

.2676288 

-  .61 

.2709308 

-  .6 

.2742531 

-  .  59 

.2775953 

-.58 

.2809573 

-.57 

.  l.  8  4  3  3  8  y 

-  .56 

.2877397 

- .  5  5 

. 291 1596 

-  .54 

.2945985 

-  .53 

.2980559 

-.52 

.3015318 

-  .51 

.3050257 

-  .5 

.3085375 

-  .49 

.3120669 

-  .43 

.3156137 

-  .47 

.3191775 

-  .46 

.3227581 

-  .45 

.3263552 

-  .44 

.3299686 

-  .43 

.3335979 

-  .42 

. 3372428 

-  .41 

. 340903 

-  .4 

.3445783 

-  .  39 

.3482633 

-  .  33 

.3519727 

-  .  37 

.3556913 

-  .36 

.3594236 

-  .35 

.3631694 

-  .34 

.3669283 

—i 

.3707 

-  .  32 

.3744842 

-  .31 

.3782805 

•  J 

.3820886 

-.29 

.3859082 

-  .  28 

.3897388 

-  .27 

.3935802 

-.26 

.3974319 

-  .  85 

.4012937 

TABLE  VII  FAILURE  PROBABILITY  VERSUS  a 


a 

F(a) 

.24 

.4051652 

.23 

.4090459 

.4129356 

.21 

.4163339 

•1 
■  <c_ 

.4207403 

.  19 

.4246546 

.18 

.4235763 

.17 

. 432505 

.16 

.4364405 

.15 

.4403323 

.  14 

.44433 

.13 

.4432332 

.12 

.4522415 

.  1 1 

.4562046 

-  .  1 

.4601721 

.09 

.4641435 

.08 

.4631135 

.07 

. 4720963 

.06 

.4760773 

.05 

.  43006 1 1 

.04 

.4840465 

.03 

. 4330335 

.02 

.4920216 

.01 

.4960106 

0 

.  5 

.01 

.5039394 

.02 

.5079784 

.03 

.511 9665 

.04 

.5159535 

.05 

.5199339 

.06 

.5239222 

.07 

.5279032 

.08 

.5313314 

.09 

.5353565 

.  1 

.5393279 

.  1  1 

.5437954 

.12 

.5477535 

.13 

.5517163 

.14 

.55567 

.15 

.5596177 

.16 

.5635595 

.17 

.567495 

.  13 

.5714237 

.  19 

.5753454 

•  L. 

.5792597 

.21 

.5331661 

YZ2 

.5370644 

cx 

F(a) 

.5909541 

24 

.5943343 

25 

. 5987063 

26 

.6025681 

27 

.6064198 

i.  o 

.6102612 

29 

.6140918 

■  -J 

.6179114 

31 

.  '6  L- 1  7  *  9  4 

32 

.6255153 

-i  “ 

_•  o 

.6292999 

34 

.6330717 

"  c 

J -J 

. 6363306 

-■6 

.6405764 

37 

.6443087 

I*  o 

.6430272 

39 

.6517317 

.4 

.6554217 

41 

.659097 

42 

.6627572 

43 

.6664021 

44 

.6700314 

45 

.6736447 

46 

.6772419 

47 

.6803225 

43 

.6843363 

49 

.6879331 

.  5 

.6914625 

51 

.6949743 

52 

.6934682 

53 

.7019441 

54 

.7054015 

erer 

,~i 

.7083403 

56 

.7122603 

57 

.7156612 

53 

.7190427 

59 

.7224047 

.  6 

.7257469 

61 

.7290691 

62 

.7323712 

63 

.7356528 

64 

.7339138 

65 

.7421539 

66 

.7453731 

67 

.7485712 

68 

.7517478 

69 

.754903 

100 
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TABLE  VII  FAILURE  PROBABILITY  VERSUS  a 


F(a) 

a 

1  .  17 

F(a) 

. 8789994 

.7 

.7580364 

1  .  18 

.8809998 

.71 

.761143 

1  . 19 

.8829767 

.72 

.7642376 

1.2 

.8849302 

.73 

.767305 

1.21 

.8868605 

.74 

.77035(31 

1  .22 

.8887675 

.75 

.7733727 

1  .23 

.8906513 

.76 

.7763723 

1  .24 

.8925122 

.77 

.7793501 

1  .25 

.8943501 

.73 

.7823046 

1  . 26 

.8961652 

.79 

.7852362 

1 .27 

.8979576 

o 

•  o 

.7831447 

1  . 28 

.8997273 

.31 

.79103 

1  . 29 

.9014746 

o-? 

.793392 

1 .3 

.9(331994 

.83 

.7967307 

1.31 

.904902 

.34 

.7995459 

1  .32 

.9065824 

.35 

.3023375 

1  . 33 

.9082408 

.  S  6 

.3051055 

1  .34 

.9098772 

.37 

.3078498 

1  .35 

.9114919 

.  8  8 

.3105704 

1  .36 

.9130849 

.89 

.3132671 

1.37 

.9146564 

.9 

.8159399 

1  .38 

.9162066 

.91 

.8185338 

1  .39 

.9177354 

.92 

.8212136 

1.4 

.9192432 

.93 

.8238145 

1  .41 

.92073 

.94 

.3263912 

1.42 

.9221961 

.95 

. 8289439 

1  .43 

.9236414 

.96 

.3314724 

1  .44 

.9250662 

.97 

“  .8339763 

1  .45 

.9264706 

.8364569 

1 .46 

.9278549 

.99 

.8339129 

1  .47 

.929219 

1 

.8413447 

.9305633 

1  .01 

.3437523 

1.49 

.9318878 

1  .02 

.8461358 

1.5 

.9331927 

T703 

'  .843495 

~  1 V5 1 

. 93447S2 

1 .04 

.85083 

1 .52 

.9357444 

1 .05 

.8531409 

i  1.53 

.9369916 

1 .06 

.8554277 

j  1.54 

.9382198 

1  .07 

.3576903 

1.55 

.9394292 

1 .08 

.8599289 

- 1.56 

.94062 

-r-.-m  ‘ 

—  .8621434 

1  .57 

.9417924 

1 .  1 

.8643339 

- 1.58 

.9429466 

-1  .11 

.8665004 

1 .59 

.9440826 

1 .12 

.868643 

- 1.6 

.9452007 

“IT  13 

.8707618 

1.61 

.9463011 

1.14 

.8728568 

- T762 

.9473839 

r.  i5 

-  . .874928 

1.63 

.9484492 
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TABLE  VII  FAILURE  PROBABILITY  VERSUS  a 


cx 

FCct) 

1.6  4 

.9494974 

1.65 

.9505285 

1.66 

.9515428 

1  .67 

.9525403 

4  /  «~k 

A  «  O  O 

.9535213 

1.69 

.954486 

1 .7 

.9554346 

1.71 

.9563671 

1  .72 

.9572838 

1 . 73 

.9581849 

1 .74 

.9590705 

1.75 

.9599409 

1.76 

.9607961 

1.77 

.9616365 

1 .78 

.962462 

1 . 79 

.9632731 

1  .8 

.9640697 

1.31 

.9648521 

1  « 

.9656205 

1  Q  O 

.9663751 

1  .84 

.9671159 

1 . 85 

.9678433 

1 . 86 

.9685573 

1  .87 

.9692581 

1  .88 

.969946 

1.89 

.9706211 

1.9 

.9712835 

1.91 

.9719334 

1.92 

.9725711 

1.93 

.9731966 

1  .94 

.9738102 

1.95 

.974412 

1 . 96 

.9750021 

1  .97 

.9755809 

1.98 

.9761483 

1  .99 

.9767046 

2 

.9772499 

2.(31 

.9777845 

2.02 

.9783084 

2.03 

.9788218 

2.04 

.9793249 

2.05 

.9798178 

2.06 

.9803008 

;z .  07“ 

.9807739 

2.08 

.9812373 

Z  .09 

.9816912 

2.1 

.9821356 

a 

F(cO 

2.11 

.9825709 

2.12 

.982997 

r>  i  o 

.9834142 

2.14 

.9838227 

Z .  15 

.9842224 

2.16 

.9846137 

2.17 

.9849966 

2.18 

.9853713 

2.19 

.9857379 

2 . 2 

.9860966 

2.21 

.9864475 

•">  •-> 

iL  ■ 

.9867906 

Z  •  23 

.9871263 

Z  .24 

.9374546 

Z .  25 

.9877756 

2 .26 

.9380894 

2.27 

.9883962 

.  4.3 

.9886962 

2.29 

.9889394 

2.3 

.9892759 

2.31 

.9895559 

.  o 

.9898296 

^  O'*' 

t_  .  o 

.9900969 

2.34 

.9903581 

2 . 35 

.9906133 

2 . 36 

. 9908625 

2.37 

.991106 

.9913437 

2.39 

.9915758 

2.4 

.9918025 

2.41 

.9920237 

2.42 

.9922397 

2.43 

.9924506 

2.44 

.9926564 

2.45 

.9928572 

2 . 46 

.9930531 

2.47 

.9932443 

2.43 

.9934309 

2-.  49 

.9936123 

2.5 

.9937903 

2.51 

.9939634 

2.52 

.9941322 

2.53 

.9942968 

2.54 

.9944573 

Z  .55 

.9946138 

Z  .56 

.9947664 

Z75T 

“  ."794915 

TABLE  VII  FAILURE  PROBABILITY  VERSUS  ct 


rr 


a_ 

F(a) 

a 

F(a) 

2.58 

.99506 

3.05 

.9988557 

2.59 

.9952012 

3.06 

.9988932 

2.6 

.9953338 

3.07 

.9989296 

2.61 

.9954728 

3.08 

.9989649 

L.  •  O  i- 

.9956035 

3.09 

.9939991 

2  •  63 

.9957307 

3.1 

.9990323 

2.64 

.9953547 

3.11 

.9990645 

2 . 65 

.9959754 

3.12 

.9990957 

2 . 6  6 

. 9960929 

3 . 1 3 

.9991259 

2.67 

.9962074 

3.14 

.9991552 

2.63 

.9963188 

3.15 

.9991836 

2.69 

.9964273 

3.16 

.9992111 

2.7 

. 996533 

3.17 

.9992377 

2.71 

.9966353 

•-i  r"v 

•2*  •  X  O 

.9992636 

2.72 

.9967359 

3.19 

.9992336 

■—  -T--5 

i_  *  f 

.9963332 

O  •  *1 

.9993123 

2.74 

.996923 

3.21 

.9993363 

2.75 

.9970202 

■-»  *-> 

■  4_  4_ 

.999359 

2.76 

.9971099 

»  4-  O 

.999381 

2.77 

.9971971 

o  •  4iA 

.9994023 

2 . 78 

.997282 

3  *25 

.9994229 

2.79 

.9973645 

3.26 

.9994429 

Z  •  s 

.9974443 

3.27 

.9994622 

•7  O  1 

i-  ■  1 

.9975229 

3.23 

.9994809 

•  O  4- 

.9975983 

3.29 

.999499 

4_  lO  O 

.9976725 

•-*  •  ..j 

.9995165 

-•  o  4 

4_  • 

.9977443 

3.31 

.9995335 

•7-  OS 

4_  ■  *v-.J 

.997314 

o  o  7 

.9995499 

2.86 

.9973317 

,  .  -2* 

.9995657 

iL  ■  O  / 

.9979476 

3.34 

.999581 

2 .  8  8 

.9930116 

3.35 

.9995959 

2.89 

.9930737 

3.36 

.9996102 

2.9 

.9931341 

3.37 

.9996241 

£.91 

.9931928 

3 . 38 

.9996375 

2 . 9  2 

.9932498 

3.39 

.9996505 

2.93 

.9983051 

3.4 

.999663 

2.94 

.9983589 

3.41 

.9996751 

2.95 

.9934111 

’3.42 

.9996868 

2.96 

.9934617 

3.43 

.9996982 

2.97 

.9985109 

3.44 

.9997091 

2 . 93 

.9985587 

“  3.415 

.9997197 

2.79 

.998605 

3.46 

.9997299 

3 

.99865 

3.47 

.9997397 

3.01 

.9936937 

3.43 

.9997492 

3.02 

.9987361 

k  3.49 

.9997584 

3.03 

.9937772 

3.5 

.9997673 

3.04 

.998817 

‘  3.51  " 

.9997759 
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IABLE  VII  FAILURE  PROBABILITY  VERSUS  a 


F(cQ 


9997842 
9997922 
9997999 
9998073 
9998145 
9998215 
9998282 
9998346 
9998408 
9998469 
9998527 
9998583 
9998636 
9998688 
9998739 
9998787 
. 9998833 
,9998878 
.9998922 
.9998963 
.9999004 
.9999042 
. 999908 
.9999116 
.999915 
.9999183 
.9999216 
.9999246 
. 9999276 
.9999305 
.9999333 
.9999359 
.9999385 
.9999409 
.9999433 
.9999456 
.9999477 
.9999499 
.9999519 
.9999538 
.9999557 
.9999575 
.9999592 
.9999609' 
.9999625 
.999964 
.9999655 


3.99 

4 


9999669 

,9999633 


* 


HV*' 


I 

I 

,  G.  HP-34C  PROGRAMS 

The  following  program  listing  PL-15  is  to  be  used  to  compute  the 
Incomplete  Gamma  Functions,  erf  (a)  and  its  inverse,  and  to  solve  the 
accelerated  test  Transcendental  function  (which  is  discussed  in  a  later 


section) . 


l5 


HP-34C 

25  13  11 

LBL  A 

h 

23  5 

STO  5 

24  3 

RCL  3 

25  3 

y* 

24  0 

RCL  0 

61 

X 

32 

CHS 

rranscendental 

1 

1 

Function 

51 

+ 

f(x) 

24  4 

RCL  4 

25  3 

yx 

24  2 

RCL  2 

61 

X 

32 

CHS 

24  5 

RCL  5 

51 

+ 

1  ) 

25  12 

RTN 

25  13  0 

LBL  0 

15  3 

x^ 

t 

2 

2 

71 

; 

erf (a)  + 

32 

CHS 

-1  . 

15  1 

e* 

erf  ^(a) 

25  12 

RTN 

25  13  1 

LBL  1 

13  12 

GSB  B 

24  6 

RCL  6 

41 

- 

25  12 

RTN 

25  13  12 

LBL  B 

23  7 

STO  7 

0 

0 

31 

ENTER  f 

24  7 

RCL  7 

14  72  0 

SO 

2 

2 

25  73 

IT 

erf (a)  + 

61 

X 

-1 

14  3 

erf  (ot) 

25  2 

1/x 

1 

61 

X 

1 

9C  19 

RTN 

f 

4  J  i-  i- 

25  13  4 

LBL  4  i 

Incomplete  Gamma  Functions 

23  8 

STO  8  1 

y(a,r) 

15  22 

R  + 

23  9 

STO  9 

0 

0 

31 

ENTER  t 

24  8 

RCL  8 

14  72  2 

n 

t'  "te-  t  f  ».ik* 
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H.  ACCELERATED  RANDOM  TEST  LEVEL  COMPUTATION 


The  computation  methods  in  this  section  calculate  the  accelerated  random 
input  vibration  level  that  will  cumulate  the  same  fatigue  damage  to 
a  structural  element  being  stressed  for  T2  hours  as  a  random  input 
vibration  level  x^  for  T1  hours.  Fracture  Mechanics  effects  causes  the 
relationship  between  parameters  to  be  a  transcendental  function  as  follows: 


1/nB 


x„  -  bA  [  1  -  b2  X2  nCe-2)  ]  »  0 


x2  is  computed  by  solving  the  above  function.  These  programs  use  the 
half-interval  technique. 

INPUT  PARAMETERS: 

a A  =  initial  crack  half-length  (inches);  NOTE:  a^  must  be 

>0.004  inches 

x^  =  input  acceleration  rms  level  at  the  service  environment  (g  rms 

x 2  =  input  acceleration  rms  level  at  the  accelerated  test 

environment  (g  rms) 

=  rms  stress  per  g  rms  constant  (ksi/g  rms) 

*  service  environment  duration  (time  units) 

=  accelerated  test  duration  (same  time  units  as  T^) . 

=  damping  linearity  constant 
=  material  crack  growth  rate  constant 
=  geometrical  parameter 

=  material  sinusoidal  fatigue  curve  constant  (ksi) 

=  material  random  fatigue  curve  constant  (ksi) 

=  material  fracture  toughness  (ksi  /In) 
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The  constant  C4  is  the  ratio  of  the  rms  stress  in  the  structural  element 


that  is  cumulating  fatigue  damage  to  the  "black  box"  rms  vibration  input 
acceleration  level.  These  programs  assume  that  the  power  spectral  density 
of  the  input  acceleration  retains  the  same  shape  at  both  x^  and  X£  levels, 
expecially  in  the  vicinity  of  structural  resonant  frequencies.  Otherwise 
C4  will  take  on  a  different  value  at  each  input  acceleration  level.  This 
is  because  the  rms  stress  level  in  the  structural  element  (hence,  fatigue) 
is  related  to  the  square  root  of  the  power  spectral  density  of  the 
acceleration  input  in  the  regions  of  the  resonant  frequencies. 

The  constant  n  denotes  the  linearity  of  the  structural  assembly  in  relating 
the  rms  stress  level  o  of  the  structural  element  to  the  input  acceleration 
level  as  follows: 

=  C4  ;  02  =  C4  X2^ 

n  =  1  represents  the  damping  linearity  system,  n  =  0.833  represents  a 
system  whose  damping  is  controlled  by  internal  stress-strain  hysteresis 
damping,  q  =  1.5  represents  a  system  controlled  by  Coulomb  friction 
damping. 


I 


% 

{■m 

’  V 


Define 


CONSTANTS : 


J_  fAKcl 
ai  Lc5yJ 

(4)- 

9-2 

W 


b3  =  1  -  b2 


n(e  -  2) 


•  ■■  (C  $ 


x2  =  bl  2n 


For  convenience  let  be  also  referred  to  as  x.  The  value  of  x 

(i.e.  x^)  to  be  solved  for  is  that  value  that  sets  the  following  function 
equal  to  zero: 


f(x)  =  x  -  b4  [  1  -  b2  xn(0  "  2)  ] 


>  t+f  ^ 


■fc-  * 


is  the  largest  value  of  x  that  will  not  cause  computational 


x2 

max 

problems  (e.g.  in  of  a  negative  number).  It  represents  the  largest 

value  of  x  that  has  practical  use.  If  >  X2  ,  fatigue  failure 

max 

will  occur  during  the  application  of  the  first  vibration  cycle.  A 

similar  practical  limit  is  imposed  on  the  selection  of  the  value  for 

a..  If  a.  is  chosen  larger  than  ar  (see  below),  x~  will  be  less 
1  1  ^2  ^ 

than  x^  ,  which  is  meaningless  for  an  "accelerated"  test. 


a  =  critical  crack  size  at  x„ 
2  L 


AKc_ 
C5  x2 


n 


inches 


1.  BASIC  LANGUAGE  (PL-16) 


PL-16  solves  for  x. 


INPUT  DATA: 


Al,  C4 ,  N,  0,  K,  Y,  C,  A,  T2,  Tl,  Xl 


represents 


C4 ,  p,  6,  AKC>  Y,  C,  A,  T2»  T^* 


OUTPUT  DATA: 


X2  =  x. 


EXAMPLES : 


WHAT  ARE  Al  •  L  4  »  N  » U  •  r!  t  Y  » L  *  A  »  Tk‘  i  T  t »  X  1 
V  .  007  ili.  8  ••  3 » 4 » 2  0  »  1  .  77  * :  '0 »  1 80  »  1  r  1 000  »  1 


X  Z  =  7.869177 


read  v 
«RUN 


WHAT  ARE  A 1 »  C'  4  »  N  »  0  *  K  >  Y  »C»A»T2»T1*X1 
.  6 1 1 1  •  o  ii  o  *  4  f  2  0  »  1  .  7  7 » 80 »  1 80 »  1  »  1000?  1 

XZ  =  Z. 898453 


y*». 


PL-16 


10  PRINT  "  WHAT  ARE  A1 » C4» N» 0> K»  Y » C» A»  T2*  T 1  *  XI " 

20  INPUT  A1,C4,N,0»K,Y,C. A.T2.T1, XI 

30  E-IE-9 

40  C5=2*A*C4/C 

50  B1=(K/(C5*Y) >**2/Al 

60  B2=( 1/BI >**< <0-2>/2> 

70  B3=Xl**(N*<0-2) )*(-B2)+l 
80  B4=  (1/B3)**<i/<N*0>> 

90  B4=B4*X1*(T1/T2)**<1/(N*0> ) 

100  M=B1**(1/(2*N) ) 

110  DEF  FNX(X)=X-B4*(l-B2*X**(N*(0-2> > >**(1/(N*G> ) 
120  H=M 

130  X4=FNX ( X )*FNX (H ) 

140  IF  X4>0  THEN  230 
150  H=H/2 

160  X3=FNX ( X ) *FNX ( X+H ) 

170  IF  X3<  0  THEN  190 
180  X=X+H 

190  IF  H<E  THEN  210 
200  GO  TO  150 
210  PRINT  ,,X2=";X 
220  GO  TO  240 

230  PR I NT "PHYSICALLY  IMPOSSIBLE  DATA  INPUT  SET" 

240  END 
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2.  TI-59  (PL-17) 


Enter  PL-17  into  the  computer;  then  execute  the  program  as  shown  in 
the  following  example  : 


INPUT 


.RAMETER 

ENTER 

PRESS 

a . 
i 

7  x  10-3 

R/S 

X1 

1 

R/S 

C4 

1 

R/S 

i/t2 

103 

R/S 

n 

.833 

R/S 

6 

4 

R/S 

A 

180 

R/S 

C 

80 

R/S 

Y 

1.77 

R/S 

AKc 

20 

R/S 

PRESS 

R/S 


DISPLAY 


x  2  =  7.8722755 


*1 

fi 


For  =  0.6  inches  ~  3.8934211  g  rms 

Execution  time  -  2  minutes 


W  ’  -tv, .  .  -  ii<K 


TI-59  PL-17  LISTING 


LRN 

LOC  CODE  KEY  LOC  CODE  KEY 


I 


PI, -17 


LOC 

CODE  KEY 

LOC 

CODE  KEY 

i  nr* 

fi4  fi  4 

150 

4  3  ■'  ’=.  L 

1  fi  i 

cr  a. 

i  cr  r 

fi  A  M  - 

4  - 

•  •  ’  :  t  fl 

1  ",  A 

J.  H 

104 

41  41 

i  -• 

4  "  _  L 

1  j  )  . 

i  s" 

3  : 

1  u 

-  1 

1  4:  n 

"4  34 

1  1  1 

-  cr  . 

1  4:  ' 

3  5  i  ; 

J.  ’  '  H 

1  4  *-• 

54 

1  1  - 

•-  t  :  : 

i  4 .  -! 

45 

114 

- 

1  64 

4  3  90  L 

i  cr 

cr.  a 

165 

3  3  4  3 

i  j  j— . 

4  a  r  ■_  ' _ 

1  A  A 

A  4 

1  1  " 

:  p-  :  4. 

167 

54 

1  7  O 

■-.cr 

1  Ah 

1  1  9 

5  a 

1  A  H 

cr.  3:  cr 

H  £ 

1  f  ; 

4  K  L'  L 

1  *7  1 

i  "?  A 

i  :  L_ 

4  3«  K  L  L 

1  A  A 

d 

i  -2= 

16  1 6 

1  A  4 

4 -  t  n 

174 

4  1  ■ 

1  -4 

4  -7=  4  -2= 

1  "7  p 

4  3  RCL 

1  4  A 

cr. A 

1  76 

i  A  i  -A 

1  A*? 

— 

i  7 

5  4  j 

1  -  H 

a 

i  •  A 

54 

129 

r,  L  i_ 

1  79 

42  8  10 

1  f  |  i 

4  b  2t 

1  o  •• 

4  j  i  ': 

i  O  i 

A  A  =! 

1  oo 

4  a  p  f  | 

1  A3 

fi  SJ  ij 

i  As  A 

21  1  21  1 

183 

54 

1  3  4 

54 

■J  l“:  .-! 

i  O  H 

4  2  STD 

i  A  ^ 

6  5 

1  C«  cr 

i 

i  j  7  i  J  9 

i  3  6 

43  RCL 

i  y  6 

43  RCL 

i  y 

A  A  J  -I; 

1 8  7 

3  A  36 

138 

.  J  M-  / 

188 

42  STD 

1 3  3 

'■4  ui  "4-  — 

189 

1  U  1  U 

1 40 

85  + 

1  90 

141 

0 1  1 

191 

43  R  L  L 

142 

54  ) 

192 

1  U  1  U 

143 

4  4  H  i  D 

193 

55 

144 

3‘  4  3  4 

194 

U2  4 

145 

53  <: 

195 

54  ) 

146 

43  RCL 

196 

42  STD 

147 

0 1  0 1 

197 

1 0  1 0 

148 

65 

198 

53  < 

149 

•  J  0 

199 

43  RCL 

lib 


ACCELERATED  x2  (HP-34C) 

The  following  program  listing  solves  for  the  accelerated  random  input 
vibration  level  to  accomplish  the  same  goal  as  PL-16.  PL-16  is 

written  in  BASIC  language.  PL-15  is  written  for  the  HP-34C  programmable 


Calculator.  The  constants  C^  ,  x2 


thru  b^  must  be  calculated 


separately  as  shown  below.  Then  their  values  are  entered  into  PL-15. 


INPUT  PARAMETERS: 


a^  (inches) 

Xj_  (g  rms) 

C^  (ksi/g  rms) 

Vt2 


A  (ksi) 

C  (ksi) 

AKC  (ksi  /In  ) 


CONSTANTS: 


_i__  TaKcI 
’1  -  «t  [c5Yj 

’■  -(4)=. 

0-2 

-dr)  ^ 


n(e  -  2) 


b^  =  1  -  62 


■ (C  $ 


X  =  b  2r* 

x2  1 
max 


'  4  V  r 


I 


,'kf^  I 


HP-34C 


n(0  -  2) 


f(x)  =  x  -  b,  1  -  b-  x 


EXAMPLE: 


HP-34C 


t>2  =  1.1102214  x  10"J  * 
b.  =  9.9888978  x  10-1 


b.  =  7.952523 

4 

n(e  -  2)  =  1.666 
l/ne  =  3.0012005  x  10* 


x2  =  59.35586 

max 


x2  =  x 


f(x)  =  Rc 


-  R2  [l  -  R 


0  X  R5 


R3  R4 


Two  initial  guesses:  x.  ;  x?  (i.e.  1;  59.35) 

1  zmax 

Program  uses  17  lines.  Refer  to  the  following  page  for  execution 
instructions.  The  solution  X2  for  the  above  example  is: 

x  2  =  7.8691  <Jl 

*  =  .007  inches  for  all  of  the  above  parameter  values. 


a^  >  .004  inches  for  these  equations  to  apply. 


TT 


HP-3AC  (PL-15) 


Turn  On 

Enter  values  of  desired  f(x) 
coefficients  into  Rq  -+•  R^ 

Enter  two  initial  guesses  for  x: 
*1 

ENTER  t 

X2  * 
max 

f  SOLVE  A 


*The  value  for  x2  entered  should  be  slightly  less  than  the  actual 
max 

value  of  X2  .  Otherwise  the  program  will  take  the  in  of  a  negative 


max 

number  (which  is  illegal)  and  "ERROR  0"  will  appear.  EXAMPLE: 
3.0199.  For  x2  enter  the  value  of  3.01. 


x2 

max 


max 
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I .  SYMBOLS 


Y(a.  t) 
Q(u,t) 


Q(x  |v) 
P(X2|v) 
Qo (a.t) 


erf p(a) / 

> 

erf (a) 


erf  (a) 


Gamma  Function  with  argument  a 
variable 

dummy  variable,  variable 

variable 

variable 

Incomplete  Gamma  Function  with  argument  parameters  a  and  t. 

Incomplete  Gamma  Function  with  argument  parameters  a  and  t. 

variable 

variable 

variable 

variable 


Incomplete  Gamma  Function  Parameters 


Normalized  Incomplete  Gamma  Function  with 
argument  parameters  a  and  t. 

interpolation  variable 

interpolation  variable 


error  function  defined  by  Papoulis  2 


inverse  error  function 


probability  of  failure  in  N  stress  cycles 


variable 


S1’  •"no 


constants 


variable 


s - 


